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Background:  This  effort  was  a  continuation  of  the  joint  ONR  and  DARPA  programs  of 
the  PI  and  Co-PI’s  initiated  as  part  of  the  DARPA  Friction  Drag  Reduction  Program 
(ATO/TTO).  The  purpose  of  the  investigation  was  to  examine  the  physics  and 
engineering  of  friction  drag  reduction  methods  for  turbulent  boundary  layers  (TBL) 
found  in  hydrodynamic  flows.  Three  methods  of  friction  drag  reduction  (FDR)  were 
examined: 

•  Polymer  Drag  Reduction-  solutions  containing  extensible,  long-chain  molecules 
are  injected  into  a  TBL.  The  polymer  molecules  interact  with  the  underlying 
turbulent  flow  and  lead  to  a  reduction  of  the  correlated  velocity  fluctuations  and, 
hence,  a  reduction  of  the  turbulent  transport  of  momentum  across  the  TBL.  This 
leads  to  local  drag  reduction  of  up  to  ~70%  for  a  TBL  flow  over  smooth  surfaces 
compared  to  a  flow  without  polymer  injection. 

•  Micro-bubble  Drag  Reduction-  air  is  injected  into  the  TBL  with  the  aim  of 
producing  a  large  volume  of  small  bubbles  in  the  near-wall  region  of  the  TBL. 
The  presence  of  the  bubbles  reduces  the  bulk  density  of  the  fluid  near  the  wall  and 
possibly  alters  the  turbulent  velocity  fluctuations.  These  effects  can  combine  to 
locally  reduce  the  friction  drag  over  80%. 

•  Air  Layer  Drag  Reduction-  air  is  injected  into  the  TBL  with  the  aim  of  creating  a 
stable  gas  layer  of  very  high  void  fraction  (>  80%  void  fraction),  separating  the 
liquid  flow  from  the  solid  surface.  This  results  in  friction  drag  reductions  of  over 
80%  compared  to  the  friction  drag  of  an  unmodified  TBL. 
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These  FDR  methods  have  been  studied  for  some  time  on  the  laboratory  scale,  yielding  a 
great  deal  of  information  on  their  engineering  potential  and  underlying  methods. 
However,  the  flow  physics  responsible  for  FDR  does  not  easily  scale  from  the  laboratory 
to  full-scale  hydrodynamic  applications.  Hence,  our  effort  focused  on  conducting  a  series 
of  experiments  at  both  large  size  scales  and  high  Reynolds  numbers  in  order  to  further 
explore  the  flow  physics  and  practicality  of  these  FDR  methods. 

Our  experiments  were  performed  on  a  canonical  TBL  flow  formed  on  a  flat  plate,  zero- 
pressure  gradient  boundary  layer.  The  maximum  length  of  the  test  model  was  ~  10  m, 
and  the  maximum  velocity  of  the  flow  was  ~  20  m/s,  making  the  Reynolds  number  based 
on  maximum  length  ~  200  million.  Our  experiments  were  conducted  in  the  William  B. 
Morgan  Large  Cavitation  Channel  (LCC),  owed  and  operated  by  the  Naval  Surface 
Warfare  Center-Carderock  Division  of  the  U.  S.  Navy. 

During  this  portion  of  our  research  program,  three  rounds  of  testing  in  the  LCC  were 
conducted: 

•  Phase  III  (Polymer  FDR) 

•  Phase  IV  (MBDR/ALDR) 

•  Phase  V  (Polymer  FDR). 

Two  additional  rounds  of  testing  were  conducted  under  the  first  portion  of  ONR  and 
DARPA  support:  Phase  I  was  a  baseline  test  of  the  model,  and  Phase  II  was  on  MBDR. 
Also,  a  portion  of  Phase  V  was  devoted  to  examining  the  influence  of  roughness  on 
Polymer  FDR  and  ALDR.  This  portion  of  the  effort  was  supported  under  DARPA 
contract  “Influence  of  Surface  Roughness  on  Polymer  Drag  Reduction”  (HR001 1-06-1- 
0057),  and  these  results  are  included  in  the  final  report  for  that  effort. 

All  of  the  experimental  efforts  were  conducted  in  coordination  with  modeling  efforts 
conducted  by  DARPA  contractors.  The  planning,  conduct,  and  data  analysis  were  jointly 
coordinated,  and  the  experimental  results  of  this  test  program  were  used  to  verify  and 
validate  the  model  development  for  polymer  FDR  and  MBDR.  The  results  of  these 
comparisons  are  not  presented  here,  but  are  presented  in  the  final  reports  of  those 
contractors.  In  this  report,  the  principal  findings  of  our  effort  are  summarized. 
Moreover,  a  more  detailed  presentation  of  our  methods,  analysis,  and  results  are  included 
as  appendices  (which  represent  manuscripts  intended  for  archival  publication  in  the 
technical  journals). 

Baseline  High  Reynolds  Number  Flow:  The  flow  over  the  test  model  represents  a  high- 
quality  data  set  combining  friction  drag  and  near-wall  velocity  data  for  a  high  Reynolds 
number  TBL.  This  data  and  analysis  are  presented  in  Appendix  A,  and  this  work  has  been 
submitted  for  publication  in  the  Journal  of  Fluid  Mechanics. 

•  Oweis,  G.  F.,  Winkel,  E.  S.,  Cutbirth,  J.  M.,  Ceccio,  S.  L.,  Perlin,  M.,  and 
Dowling,  D.  R.,  “Smooth-Flat-Plate  Turbulent  Boundary  Layer  Measurements  at 
High  Reynolds  Number,”  Journal  of  Fluid  Mechanics,  (submitted)  (2008) 


Polymer  Drag  Reduction:  The  injection  of  solutions  containing  long-chain  polymers 
into  the  TBL  will  result  in  significant  reductions  in  friction  drag,  which  is  a  well- 
documented  phenomenon.  The  use  of  polymer  injection  for  external  flows  presents  a 
singular  challenge.  In  external  flows,  the  mixing  of  the  injected  polymer  in  the  TBL  is 
reduced  at  the  highest  levels  of  drag  reduction.  But,  the  mixing  across  the  TBL  never 
ceases.  Hence,  the  near-wall  concentration  of  the  polymer  will  be  continually  reduced 
downstream  of  the  injection  location  (usually  a  line  source  perpendicular  to  the  stream- 
wise  direction).  As  the  near-wall  concentration  of  polymer  is  reduced,  the  amount  of 
drag  reduction  is  decreased  until  the  effect  is  lost.  In  the  present  experiments,  we 
examined  a  variety  of  parameters  that  influence  the  amount  of  drag  reduction  that  was 
achieved  through  the  injection  of  aqueous  solutions  of  Poly(Ethylene  Oxide).  The 
influence  of  polymer  flux,  concentration,  molecular  weight,  and  injector  geometry  on  the 
developed  friction  drag  reduction  were  examined.  These  data  were  combined  with 
measurements  of  the  local  near-wall  polymer  concentration  and  flow  fields  to  understand 
how  the  mixing  processes  influence  the  reduced  drag.  Finally,  measurements  of  the 
polymer  solution  rheology  prior  to  and  after  injection  were  employed  to  examine  the 
processes  of  flow  induced  polymer  degradation. 

Appendix  B  presents  detailed  results  from  our  polymer  drag  reduction  effort.  This  work 
has  been  submitted  for  publication  in  the  Journal  of  Fluid  Mechanics: 

•  Winkel,  E.  S.,  Oweis,  G.,  Vanapalli,  S.  A.,  Dowling,  D.  R.,  Perlin,  M.,  Solomon, 
M.  J.  and  Ceccio,  S.  L.,  “Friction  Drag  Reduction  of  a  High  Reynolds  Number 
Turbulent  Boundary  Layer  with  Wall  Injected  Polymer  Solutions,”  Journal  of 
Fluid  Mechanics,  (submitted)  (2008) 

Other  details  of  this  work  have  been  published  in  the  following  manuscripts: 

•  Vanapalli,  S.  A.,  Ceccio,  S.  L.,  and  Solomon,  M.  J.,  “Universal  scaling  for 
polymer  chain  scission  in  turbulence,”  Proc.  National  Academy  of  Science,  Vol. 
103,  No.  45,  pp  16660-16665  (2006) 

•  Winkel,  E.  S.,  Elbing,  B.  R.,  Dowling,  D.  R.,  Ceccio,  S.  L.,  and  Perlin,  M.  “High- 
Reynolds-Number  Turbulent-Boundary-Layer  Surface  Pressure  Fluctuations  with 
Bubble  or  Polymer  Additives,”  Proc.  of 2005  ASME  I.M.E.C.E,  Orlando,  (2005). 

•  Garwood,  G.  C.,  Winkel,  E.  S.,  Vanapalli,  S.,  Elbing,  B.,  Walker,  D.  T.,  Ceccio, 
S.  L.,  Perlin,  M.,  Solomon,  M.  J.,  “Drag  Reduction  by  a  Homogenous  Polymer 
Solution  in  Large  Diameter,  High  Shear  Pipe  Flow,”  Proc.  of  2nd  Int.  Symp.  on 
Seawater  Drag  Reduction,  Busan,  Korea  (2005) 

•  Almeida,  T.  G.,  Walker,  D.  T.,  Leighton,  R.  I.,  Alajbegovic,  A.,  Pankajakshan,  R., 
Taylor,  L.  K.,  Whitfield,  D.  L.,  and  Ceccio,  S.  L.,  “A  Reynolds  Averaged  Model 
for  the  Prediction  of  Friction  Drag  Reduction  by  Polymer  Additives,”  Proc.  26'h 
Symposium  on  Naval  Hydrodynamics,  Rome,  (2006) 


Micro-bubble  Drag  Reduction:  The  injection  of  micro-bubbles  into  the  TBL  was 
examined  over  large  size  scales  and  high  Reynolds  numbers.  In  the  first  part  of  the 
DARPA  FDR  program,  we  examined  MBDR  and  found  that  injection  of  gas  would  lead 
to  large,  local  reductions  of  friction  drag.  But,  the  downstream  persistence  of  the  FDR 
was  limited.  It  was  hypothesized  that  the  reason  for  the  loss  of  drag  reduction  was  due  to 
the  stratification  of  the  bubbles  away  from  the  near-wall  region.  This  result  was  reported 
in  the  following  manuscript: 

•  Sanders,  W.  C.,  Winkel,  E.  S.,  Dowling,  D.  R.,  Perlin,  M.,  and  Ceccio,  S.  L. 
“Bubble  Friction  Drag  Reduction  in  a  High  Reynolds  Number  Flat  Plate 
Turbulent  Boundary  Layer,”  Journal  of  Fluid  Mechanics ,  Vol.  552,  pp  353-380 
(2006) 

•  Kunz,  R.  F.,  Gibeling,  H.  J.,  Maxey,  M.  R.,  Tryggvason,  G.,  Fontaine,  A.  A., 
Petrie,  H.  L.,  and  Ceccio,  S.  L.,  “Validation  of  Two-Fluid  Eulerian  CFD 
Modeling  for  Microbubble  Drag  Reduction  across  a  Wide  Range  of  Reynolds 
Numbers,”  Journal  of  Fluids  Engineering,  Vol.  129,  No.  1,  pp  66-79  (2007) 

In  this  phase  of  the  FDR  program,  these  experiments  were  repeated  to  determine  the 
influence  of  surfactants  and  injector  geometry.  Moreover,  optical  and  electrical 
impedance  measurements  were  conducted  to  specifically  examine  the  near-wall  void 
fraction.  These  data  confirmed  that  a  region  largely  devoid  of  bubbles  (a  “liquid  layer”) 
formed  near  the  surface  of  the  model  when  the  friction  drag  reduction  was  lost.  Also,  it 
was  found  that  the  presence  of  surfactants,  while  leading  to  a  reduction  in  the  average 
bubble  size,  did  not  significantly  alter  the  observed  drag  reduction.  In  addition  changes 
in  the  gas  injection  method  had  little  effect. 

Appendix  C  presents  detailed  results  from  our  MBDR  effort.  This  work  has  been 
submitted  for  publication  in  the  Journal  of  Fluid  Mechanics : 

•  Belbing,  B.  R.,  Winkel,  E.  S.,  Lay,  K.  A.,  Ceccio,  S.  L.,  Dowling,  D.  R.,  and 
Perlin,  M.,  “Bubble-Induced  Skin-friction  Drag  Reduction  and  the  Abrupt 
Transition  to  Air-layer  Drag  Reduction,”  Journal  of  Fluid  Mechanics,  (submitted) 
(2008) 

Other  details  of  this  work  have  been  published  in  the  following  manuscripts: 

•  Winkel,  E.  S.,  Elbing,  B.  R.,  Ceccio,  S.  L.,  Perlin,  M.,  Dowling,  D.  R.,  “High- 
Reynolds-number  turbulent-boundary-layer  wall  pressure  fluctuations  with  skin 
friction  reduction  by  gas  injection,”  Journal  of  the  Acoustical  Society  of  America, 
(in  press)  (2008) 


Air  Layer  Drag  Reduction:  If  enough  air  is  injected  into  the  TBL,  the  gas  can  coalesce 
into  a  layer  under  the  influence  of  buoyancy.  This  was  observed  in  the  first  round  of 
MBDR  testing  at  lower  speeds,  where  the  friction  drag  was  reduced  by  over  80%  over  the 
full  stream-wise  extent  of  the  test  model.  In  the  follow-on  experiments,  we  observed  this 
Air  Layer  Drag  Reduction  at  higher  speeds,  with  increased  gas  flux.  The  critical  gas  flux 
and  dynamics  of  the  air  layers  were  examined  for  flow  over  the  smooth  HIPLATE. 

Appendix  C  presents  detailed  results  from  our  ALDR  effort.  This  work  has  been 
submitted  for  publication  in  the  Journal  of  Fluid  Mechanics: 

•  Belbing,  B.  R.,  Winkel,  E.  S.,  Lay,  K.  A.,  Ceccio,  S.  L.,  Dowling,  D.  R.,  and 
Perlin,  M.,  “Bubble-Induced  Skin-friction  Drag  Reduction  and  the  Abrupt 
Transition  to  Air-layer  Drag  Reduction,”  Journal  of  Fluid  Mechanics,  (submitted) 
(2008) 


Conclusions:  The  UM  FDR  effort  helped  reveal  the  underlying  physical  process 
responsible  for  the  development  and  loss  of  friction  drag  reduction  with  the  injection  of 
polymer  solutions  and  gas  into  a  high  Reynolds  number  TBL.  At  the  conclusion  of  our 
study,  we  have  come  to  the  conclusion  that  all  three  drag  reduction  methods  can  lead  to 
significant,  local  reductions  in  friction  drag.  However,  only  ALDR  and  PDR  remained  as 
methods  that  could  lead  to  persistent  reduction  of  friction  drag  over  large  stream-wise 
extents  on  our  test  model.  We  continue  to  perform  cost-benefit  analysis  of  both  ALDR 
and  PDR.  Our  analysis  of  PDR  has  led  us  to  conclude  that  it  is  not  cost  effective,  given 
the  necessary  flux  of  polymer,  the  rate  of  polymer  mixing,  and  the  relative  cost  of 
polymer  compared  to  typical  fuel  costs.  Our  analysis  of  ALDR  has  led  us  to  conclude 
that  ALDR  may  be  cost  effective  if  the  gas  is  not  pumped  into  a  prohibitively  large  back¬ 
pressure,  and  if  the  air  layers  can  made  to  persist  for  longer  distances  than  observed  on 
our  present  test  model. 
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Abstract 

The  characteristics  of  smooth-flat-plate  turbulent  boundary  layers  and  their 
dependence  on  Reynolds  number  have  been  studied  for  nearly  a  century.  However,  there 
is  a  relative  dearth  of  measurements  at  Reynolds  numbers  typical  of  full-scale  marine  and 
aerospace  transportation  systems  (Reg=Ued/v>  105,  where  Ue  is  the  exterior  or  free 
stream  flow  speed,  6  is  the  momentum  thickness  of  the  boundary  layer,  and  v  is  the 
kinematic  viscosity  of  the  fluid).  This  paper  presents  the  results  of  a  new  experimental 
study  of  the  boundary  layer  that  forms  on  a  flat  plate  for  Ree  up  to  and  exceeding 
150,000.  The  experiments  were  conducted  in  the  William  B.  Morgan  Large  Cavitation 
Channel  (LCC)  at  nominal  water  flow  speeds  from  3  to  20  m/s  on  a  flat-plate  test  model 
12.9  m  long  and  3.05  m  wide.  The  model’s  test  surface  was  polished  to  k/lv  =k+  <0.2, 
where  k  is  the  root-mean-square  roughness  height  of  the  surface,  lv  =  v«Jp/ is  the 

usual  viscous  wall  unit,  p  is  the  fluid  density,  and  rw  is  the  wall  shear  stress.  The  flow 
closely  approached  the  zero  pressure  gradient  condition;  side-wall  and  model-surface 
boundary-layer  growth  in  the  LCC’s  test  section  caused  a  streamwise  flow  speed  increase 
of  ~2.5%  over  the  11-m-long  test  surface.  The  overall  data  set  is  uniquely  complete  for 
an  experimental  investigation  in  this  Reynolds  number  range.  Direct  measurements  of 
static  pressure  and  wall  shear  stress  vs.  downstream  distance  were  obtained  with  pressure 
taps  and  floating-plate  skin-friction  force  balances.  Composite  wall-normal  profiles  of  the 
mean  velocity  and  second-order  turbulence  statistics  were  measured  with  conventional 
two-component  laser-Doppler  velocimetry  and  a  unique  near-wall  implementation  of 
particle-tracking  velocimetry.  The  resulting  flow  profiles  span  the  wall-normal 
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coordinate  (y)  range  from  y+  =  y/lv  below  unity  to  y  >  Sgg,  where  is  the  wall-normal 
location  where  99%  of  the  free-stream  flow  speed  is  achieved.  Internal  consistency 
checks  within  the  data  set  show  that  rw  determined  from  the  force  balances,  the  PTV- 
measured  mean  wall-gradient  of  streamwise  velocity,  and  the  sum  of  the  LDV-measured 
mean  viscous  and  Reynolds-shear  stresses,  all  agree  within  experimental  uncertainty. 
This  data  set  provides  a  unique  opportunity  for  evaluating  boundary  layer  scaling  laws. 


1.0  Introduction 


Turbulent  boundary  layers  are  ubiquitous  in  both  natural  and  human-engineered 
fluid  flows.  They  commonly  set  heat  and  mass  transfer  rates  in  the  ocean  and 
atmosphere,  and  are  often  a  critical  factor  in  determining  the  efficiency  and  performance 
of  transportation  systems.  The  simplest  possible  boundary  layer  forms  where  a  uniform 
unidirectional  flow  contacts  a  smooth  flat  surface.  This  boundary  layer  will  be  turbulent 
if  the  Reynolds  number  of  the  flow  is  high  enough  to  sustain  turbulence,  and  an 
intentional  or  naturally  occurring  disturbance  has  caused  transition  from  laminar  to 
turbulent  flow. 

Given  the  importance  of  the  phenomena  and  the  simplicity  of  the  flow  geometry, 
flat-plate  turbulent  boundary  layers  (TBLs)  have  been  investigated  for  the  greater  part  of 
the  last  century.  Much  is  known  about  the  mean  and  fluctuating  components  of  flat  plate 
TBL  flow  from  controlled  laboratory-scale  investigations  conducted  at  momentum- 
thickness-based  Reynolds  numbers  Ree  =  Ue6/v  <  104,  where  Ut  is  the  exterior  or  free 
stream  flow  speed,  G  is  the  momentum  thickness  of  the  boundary  layer,  and  v  is  the 
kinematic  viscosity  of  the  fluid.  Additional  investigations  conducted  in  the  atmospheric 
boundary  layer  and  in  special  testing  facilities  at  higher  Reynolds  numbers  -  where  the 
investigations  cannot  be  as  well  controlled  or  as  complete  as  the  laboratory  scale  studies 
-  have  allowed  scaling  laws  developed  from  the  laboratory  studies,  the  governing 
equations,  and  dimensional  reasoning  to  be  at  least  partially  evaluated  at  much  higher 
Reynolds  numbers.  The  primary  purpose  of  this  manuscript  is  to  report  the 
characteristics  of  the  TBL  that  develops  on  a  smooth  flat  plate  under  well-controlled 
conditions  for  6xl03  <Reg  <1.5xl05  (or,  alternatively  in  terms  of  the  downstream- 
distance-based  Reynolds  number,  3xl06  <Rex=Uex/v  <2.2x\0s).  The  experimental 
results  presented  here  are  unique  in  their  completeness  in  this  Reynolds  number  range. 

Although  the  flat-plate  geometry  is  idealized  compared  to  that  of  practical 
aerodynamic  and  hydrodynamic  devices,  the  structure  and  scaling  of  flat-plate  TBLs 
remain  important  topics  because  moderate-Reynolds-number  results  are  commonly 
scaled-up  when  designing  or  predicting  the  performance  of  fluid  dynamic  devices  that 
operate  at  high  Reynolds  numbers  where  experimental  measurements  are  sparse  or  non 


existent.  Furthermore,  the  numerical  values  and  functional  forms  used  in  computational 
models  of  wall-bounded  turbulence  are  based  on  the  flow  structure  and  empirical 
constants  (k  and  B )  determined  from  TBL  studies  (see  Launder  and  Spalding  1972,  or 
Pope  2000).  Thus,  the  validity  of  such  models  is  enhanced  when  the  empirical  constants 
are  shown  to  be  valid  over  as  large  a  Reynolds  number  range  as  possible. 

In  addition,  over  the  last  two  decades  or  so,  a  variety  of  revisions  to  the  classical 
understanding  of  TBL  scaling  and  flow  structure  have  been  derived  and/or  proposed  for 
the  mean  streamwise  velocity  profile  and  the  Reynolds  stresses.  Given  that  the  intent  of 
this  manuscript  is  to  dispassionately  and  concisely  report  new  high-Reynolds-number 
TBL  results,  these  new  experimental  data  are  tabulated  in  dimensional  form  in  the 
appendix,  and  plotted  with  traditional  inner  (law  of  the  wall)  scaling  along  with  a  few 
simple  scaling-law  tests.  For  conciseness  and  clarity,  detailed  tests  and  comparisons 
involving:  /)  revised  scaling  laws  (George  and  Castillo  1997,  Marusic  et  al.  1997, 
DeGraaff  and  Eaton  2000,  Marusic  and  Kunkel  2003,  Lindgren  et  al.  2004),  if)  recently 
proposed  revisions  to  the  layer  structure  of  TBLs  (Wei  et  al.  2005,  Fife  et  al.  2005),  and 
Hi)  log-law  and  power-law  scaling  of  the  mean  streamwise  velocity  profile  (Osterlund  et 
al.  2000a,  Barenblatt  et  al.  2000,  Osterlund  et  al.  2000b,  Afzal  2001,  Panton  2002, 
Buschmann  and  Gad-el-Hak  2003),  while  clearly  possible,  are  beyond  the  scope  of  this 
manuscript. 

The  experimental  results  presented  here  were  collected  during  three  multi-month 
experimental  campaigns  conducted  between  the  fall  of  2001  and  the  spring  of  2005.  The 
experiments  involved  free  stream  flow  speeds  from  3  to  20  m/s,  a  flat  plate  test  model 
with  a  hydraulically-smooth  test  surface  1 1  m  in  length,  and  the  world’s  largest  low- 
turbulence  water  tunnel,  the  US  Navy’s  William  B.  Morgan  Large  Cavitation  Channel 
(LCC).  Measurements  are  reported  from  static  pressure  taps,  average  skin  friction  force 
balances,  two-component  laser  Doppler  velocimetry  (LDV),  and  near- wall  particle 
tracking  velocimetry  (PTV).  In  particular,  the  static  pressure  measurements  show  a  mild 
favorable  pressure  gradient  that  is  consistent  with  boundary  layer  growth  in  the  test 
section  of  the  LCC.  But  more  importantly,  the  wall  shear  stress  determined  from  the 
average  skin  friction  force  balances,  the  LDV  Reynolds  stress  profiles,  and  measured 


mean  flow  in  the  viscous  sublayer  all  agree,  and  are  consistent  with  the  Von  Karman 
boundary-layer  momentum  integral  equation,  to  within  experimental  error. 

The  published  literature  on  the  topic  of  turbulent  boundary  layers  is  extensive. 
However,  review  articles  (Sreenivasan  1989,  Gad-el-Hak  and  Bandyopadhyay  1994, 
Femholz  and  Findley  1996)  and  recent  texts  (Pope  2000,  White  2006)  provide  a  starting 
place  for  sorting  and  assimilating  this  material.  The  prior  studies  most  relevant  to  the  one 
reported  here  fall  into  two  groups:  smooth-wall  laboratory  investigations  (Fernholz  et  al. 
1995,  Osterlund  et  al.  1999,  DeGraaff  and  Eaton  2000)  where  Ree  reaches  ~104  (or  Rex 

n 

reaches  ~10  )  and  the  rough- wall  atmospheric  and  wind-tunnel-sidewall  boundary  layer 
experiments  (Saddoughi  and  Veeravalli  1994,  Metzger  and  Klewicki  2001,  and  Metzger 
et  al.  2001)  conducted  at  much  higher  Reynolds  numbers.  The  experimental  results  from 
the  present  study  extend  the  Reynolds  number  range  of  the  laboratory  results  to  partially 
close  the  Reynolds  number  gap  between  the  two  groups  of  prior  studies. 

The  remainder  of  this  manuscript  is  organized  into  four  sections.  The  next  section 
describes  the  experimental  setup,  measurement  techniques,  and  data  processing  methods. 
The  third  and  fourth  sections  present  the  average  flow  characteristics  and  turbulence 
statistics,  respectively.  This  research  effort  is  summarized  and  conclusions  are  stated  in 
the  final  section.  The  appendix  provides  tabulations  of  the  experimental  data. 

2. 0  Experimental  Techniques 

The  experiments  were  conducted  in  the  U.  S.  Navy’s  William  B.  Morgan  Large 
Cavitation  Channel  (LCC).  This  closed-circuit  water  tunnel  is  the  largest  facility  of  its 
kind  in  the  world  (Etter  et  al.  2005).  The  LCC’s  test  section  has  a  nominal  length  of  13 
meters  and  a  cross-section  of  3.05  m  x  3.05  m.  The  LCC  may  be  operated  at  steady 
empty-test-section  flow  speeds  from  0.5  m/s  to  more  than  18  m/s  and  at  absolute  test 
section  pressures  from  3.4  kPa  to  414  kPa.  At  the  flow  speeds  of  these  tests,  the  LCC’s 
test-section  centerline  turbulence  level  was  0.2%  to  0.4%. 

The  test  model  was  the  same  one  utilized  by  Sanders  et  al.  (2006),  a  flat  plate  12.9 
m  long  and  18.4  cm  thick  that  horizontally  spanned  the  LCC  test  section  (see  Fig.  1). 
With  these  dimensions,  test  section  blockage  lead  to  flow  speeds  above  the  developing 


boundary  layer  as  high  as  20.2  m/s.  The  17-metric-ton  model  was  composed  of  three 
sections  fitted  together  with  precision-machined  slot-and-key  mechanisms  that  allowed 
repeatable  and  seamless  assembly.  The  edges  of  the  model  were  mated  to  the  test  section 
walls  with  an  inflatable  seal  to  prevent  bypass  flow  between  the  top  and  bottom  sides  of 
the  model.  In  addition,  45°  triangular-wedge  fillet  edge-fairings  (75  mm  maximum 
height  and  width)  were  installed  at  the  model-sidewall  junction  to  control  and  minimize 
cavitation,  and  streamwise  corner  and  junction  vortices.  The  leading  edge  of  the  model 
was  a  4-to-l  ellipse  while  the  trailing  edge  was  a  1 5°-full-angle  triangular  wedge  of  0.6- 
m  length  that  was  terminated  at  25  mm  thickness  with  40°  bevel  angle.  This  asymmetric 
trailing  edge  design  was  intended  to  prevent  near-wake  vortex  shedding  and  its 
propensity  for  flow-induced  model  vibration.  The  Cartesian  coordinate  system  used 
throughout  this  manuscript  is  shown  in  Fig.  1;  x  is  zero  at  the  model’s  leading  edge  and 
runs  parallel  to  the  flow  direction,  y  is  zero  on  the  test  surface  and  increases 
perpendicular  to  it  (and  to  the  x  axis)  while  z  runs  in  the  spanwise  direction  and 
completes  a  right-handed  coordinate  system. 

During  testing,  surface  normal  acceleration  was  measured  with  six  internally 
mounted  Wilcoxin  754-1  accelerometers.  When  integrated,  these  signals  produced  a 
maximum  measured  root-mean-square  (rms)  surface-normal  vibration  velocity  of  ~2 
mm/s  at  the  highest  test  speed.  This  vibration  fluctuation  level  was  far  smaller  than  the 
velocity  resolution  of  the  LDV  system  (see  below)  and  was  also  well  below  the  rms  free 
stream  velocity  fluctuation.  Thus,  model  vibration  was  negligible  in  these  experiments. 

A  distributed-roughness  boundary  layer  trip  was  used.  Otherwise,  the  model  was 
hydraulically  smooth  from  its  leading  edge  back  to  the  beginning  of  the  trailing  edge 
wedge.  The  test  surface  was  304  stainless  steel  polished  to  a  nominal  surface  roughness 
of  0.4  pm  or  less.  The  roughness  parameter,  k+  =  k/lv,  where  k  is  the  average  roughness 

height  and  lv  is  the  viscous  length  ( lv  =  v-yjr ~J~p  for  kinematic  viscosity  v,  wall  shear 

stress  Tw,  and  fluid  density  p),  was  less  than  or  equal  to  0.2  over  the  entire  flat  test  surface 
for  all  flow  conditions.  The  boundary  layer  was  tripped  by  an  intentional  spanwise  strip 
of  distributed  roughness  applied  to  the  elliptical  portion  of  the  model  starting  25  mm 
from  the  model’s  leading  edge  and  extending  250  mm  downstream.  The  distributed 
roughness  was  composed  of  nominally-0.25-mm  (100  grit)  sand  grains  embedded  in  a 


film  of  epoxy  and  spaced  randomly  2  to  5  mm  apart.  The  design  of  this  distributed 
roughness  trip  was  based  on  guidelines  provided  by  Prof.  H.  Nagib  (private 
communication). 

Static  pressure,  skin  friction,  and  fluid  velocity  measurements  were  made  at  three 
primary  free  stream  flow  speeds,  £/»  =  6.65  m/s,  13.2  m/s,  and  19.9  m/s.  These  values 
are  correct  for  the  center  of  the  test  plate;  however  -  as  discussed  in  the  next  section  -  Ue 
varied  from  ~1%  lower  (upstream)  to  ~1%  higher  (downstream)  over  the  plate’s 
stream  wise  extent.  In  addition,  a  subset  of  the  flow  measurements  were  made  at  £/«,  =  3.3 
m/s. 

Static  pressure,  P(x),  was  measured  at  1 1  downstream  (x)  locations  with  a 
Rosemount  305 IP  transducer  and  1 .6-mm-diameter  taps  located  48.3  cm  above  the  plate 
on  the  sidewalls  of  the  LCC  test  section.  Calibration  of  the  Rosemount  transducer  was 
performed  using  a  Druck  DP601  pressure  source  and  a  Paroscientific  DigiQuartz  Model 
740  pressure  transducer.  The  static  pressure  measurements  were  corrected  for  zero-bias 
(or  drift)  error  via  no-flow  pressure  measurements  performed  before  and  after  TBL  data 
acquisition.  In  the  next  section,  these  static  pressure  measurements  are  reported  in  terms 
of  the  pressure  coefficient,  Cp(x)  =  (P(x)- Px)/^pUl  where  Px  in  the  pressure  measured 
at  the  first  tap  at  x  =  1.96  m  and  Ue  is  the  free  stream  flow  speed  over  the  center  of  the 
test  plate.  Test-section  side-wall  and  model-surface  boundary-layer  growth  caused  Ue(x ) 
to  increase  by  ~2.5%  between  the  first  and  last  pressure  taps.  At  the  flow  speeds  of  these 
tests,  hole  error  (see  Benedict  1984)  was  the  dominant  source  of  uncertainty,  leading  to 
an  overall  uncertainty  of  ±0.0005  in  pressure  coefficient  units. 

Average  skin  friction,  Tw(x),  on  the  test  surface  was  measured  with  six  floating- 
element  strain-gage  force  balances  located  at  ~l/3  span  and  at  x  =  1.96  m,  3.41  m,  5.94 
m,  7.43  m,  9.23  m,  and  10.68  m.  The  round  sensitive  surfaces  of  these  balances  had  a 
diameter  of  15.24  cm  (6  inches)  were  adjusted  to  be  flush  with  test  surface  of  the  model. 
The  gap  around  the  circumference  of  the  sensing  surface  had  a  maximum  value  of  ~75 
^m.  This  gap  was  checked  for  uniformity  (and  minimum  clearance)  using  a  51  // m 
(0.002  inch)  stainless  steel  shim.  At  the  highest  flow  speeds,  P  is  approximately  1.7 
microns  in  these  experiments,  so  the  worst-case  average  and  maximum  surface  gaps  were 
30 /+  and  44 /+  respectively. 


Strain  gauge  amplifiers  (Vishay  2310)  were  employed  for  excitation,  nulling,  and 
amplification  of  the  skin  friction  signals.  Data  were  collected  at  a  rate  of  50  Hz  for 
variable  lengths  of  time,  but  measured  skin  friction  forces  were  primarily  determined 
from  10-second  averages.  The  force  balances  were  individually  calibrated  at  least  once 
per  test  day  using  a  suction  cup  attached  via  a  metal  cord  to  a  second  precision  load  cell 
(Omega  LCEB-5)  that  was  mounted  on  a  finely-adjustable  translation  stage.  Calibration 
measurements  were  made  by  temporarily  fixing  the  base  of  the  translation  stage  and  then 
adjusting  the  tension  in  the  metal  cord  with  horizontal  movement  of  the  precision  load 
cell  to  produce  known  and  repeatable  horizontal  loads  on  each  force  balance.  A  linear  fit 
to  the  calibration  data  was  used  to  convert  the  measured  voltage  signals  to  friction  force; 
the  root-mean-square  deviation  of  the  calibration  points  from  a  linear  fit  corresponded  to 
at  most  ±0.02  N.  In  the  next  section,  these  static  pressure  measurements  are  reported  in 
terms  of  the  skin  friction  coefficient,  Cf(x)  =  Tw(x)/^pU*(x). 

For  the  skin  friction  measurements,  the  significant  sources  of  uncertainty  include 
drift  and  calibration  errors.  As  is  common  for  strain  gauge  equipment,  the  force-balance 
outputs  tended  to  drift  mildly  over  periods  of  several  hours.  Thus,  the  strain-gage  bridges 
were  nulled  during  a  no-flow  condition  every  one  or  two  hours  during  a  test  day.  The 
error  associated  with  drift  was  approximately  ±0.05  N.  To  account  for  drift,  the  average 
of  at  least  12  independent  10-second  measurements  is  reported  and  the  drift  uncertainty, 
<ydrift,  is  taken  to  be  twice  the  standard  deviation  of  the  mean  of  these  12+  measurements. 
The  uncertainty  associated  with  calibration,  ocah  was  taken  to  be  twice  the  standard 
deviation  of  the  mean  of  multiple  calibrations  (typically  1  to  2%)  taken  over  a  period  of 
several  days  or  more.  The  combined  experimental  uncertainty,  a  =  -y]<y%rifi  +  o]al ,  varied 

somewhat  between  force  balances,  but  were  generally  less  than  ±3%  of  the  final  skin 
friction  measurement  from  any  force  balance. 

Fluid  velocity  measurements  were  made  with  two  LDV  systems.  The  first,  a 
single  component  system,  was  fixed  in  place  and  used  to  monitor  the  test-section-inlet 
streamwise  velocity  at  x  =  -6.2  cm  and  y  =  14.1  cm.  It  operated  with  a  2.775x  beam 
expander,  1 15  mm  beam  spacing,  and  88  cm  focal  length.  The  second,  a  two  component 
system,  was  mounted  on  a  vertical  traverse  and  used  for  measuring  the  mean  velocity 


(£/,  V)  and  the  second-order  statistics  of  the  velocity  fluctuations  (u,,vA)  in  pairs  of  vertical 
columns  (0.1  mm  <  y  <  300  mm)  near  x  =  5.9  m  and  10.7  m.  This  two-component 
system  consisted  of  Dantec  BSA  57N11  signal  processors,  fiber  optic  probes,  Spectra 
Physics  6  Watt  Argon-Ion  lasers  (model  2017),  Dantec  three-dimensional  traverse,  and 
Dantec  Flow  software.  This  system  was  operated  at  514.5  nm  and  488  nm  wavelengths 
with  a  1.5x  beam  expander,  112  mm  beam  spacing,  and  1,600  mm  focal  length.  The 
burst- processor  parameters  included  transit-time  weighted  averages  for  velocity  bias 
correction,  low  and  high-pass  filtering  dependent  on  free  stream  velocity,  and  particle 
size  rejection  for  the  elimination  of  erroneous  measurements  from  reflections  off  the 
polished  surface.  The  y-coordinate  origin  for  the  second  LDV  was  determined  within 
±0.1  mm  by  measuring  the  near  wall  fluid  velocity  with  the  LDV  beams  above,  touching, 
and  reflecting  from  the  model’s  working  surface,  and  setting  y  =  0  at  the  location  of 
minimum  average  streamwise  velocity. 

For  the  LDV  measurements,  the  water  channel  was  flood-seeded  with  silicon 
carbide  particles  with  a  nominal  diameter  of  3  microns.  The  downstream  LDV  had  a 
prolate  focal  volume  with  a  diameter  of  170  pm  and  a  length  of  6.5  mm  that  was  oriented 
with  its  long  axis  perpendicular  to  the  flow  direction  but  nearly  parallel  to  the  plate 
surface.  To  make  measurements  near  the  test  surface,  the  downstream  LDV  optical  head 
was  tilted  2.2°;  this  increased  the  projected  vertical  dimension  of  the  downstream  LDV 
focal  volume  to  250  pm.  This  increase  in  the  focal-volume  vertical  dimension  becomes  a 
dominant  factor  near  the  wall.  Here,  the  larger  dimension  translates  to  an  increased  range 
of  instantaneous  velocity  measurements  due  to  the  steep  gradient  of  the  velocity  profile 

near  the  wall  and  this  can  lead  to  elevated  u'2  values  (the  overbar  denotes  a  time  or 
ensemble  average).  In  addition,  the  LDV-head  tilt  lead  to  non-normal  optical 
transmission  at  the  vertical  LCC  sidewall  windows.  The  resulting  optical  distortion  and 
reduction  in  coincident  bursts  limited  the  LDV  focal  volume’s  distance  from  the  LCC 
sidewall  to  0.65  m.  However,  this  spanwise  location  is  well  beyond  the  furthest  influence 
of  the  model-side-wall-junction  flow  (~0.20  m  from  the  LCC  sidewall). 

Both  LDV  systems  were  calibrated  with  a  spinning  disk  covered  with  60-grit 
emery  paper  that  was  driven  by  a  CompuMotor  Model  SM233BE-NTQN  and 
CompuMotor  TQ10X  Servo-Controller.  The  major  contributor  to  the  uncertainty  of  the 


calibration  is  the  ±0.4  rps  accuracy  of  the  motor  as  listed  by  its  manufacturer.  The  scatter 
of  this  LDV  calibration  data  about  a  linear  fit  was  less  than  ±0.005  m/s.  The  overall 
velocity  uncertainty  was  ±0.025  m/s  at  the  95%  confidence  level  (Park  et  al.  2003). 

The  measurements  reported  here  are  primarily  drawn  from  the  third  phase  of 
experiments  where  10,000  LDV  bursts  were  collected  at  two  pairs  of  measurement 
locations  at  the  three  primary  flow  speeds  (6.65,  13.2,  and  19.9  m/s).  The  burst 
processors  gain  and  photomultiplier  tube  (PMT)  voltage  were  adjusted  for  each 
measurement  location  so  that  the  burst  rate  remained  relatively  constant  with  respect  to 
the  distance  from  the  plate  for  each  free  stream  velocity.  The  coincident  burst  rates 
varied  from  100  Hz  at  6.65  m/s  to  170  Hz  at  19.9  m/s.  The  coincidence  window  varied 
from  0.02  ms  at  6.65  m/s  to  0.01  ms  at  19.9  m/s,  roughly  translating  to  the  ratio  of  probe 
volume  diameter  to  free  stream  velocity.  Table  1  provides  a  listing  of  flow  parameters  at 
the  downstream  (x)  locations  of  the  LDV  profiles  were  measured.  The  minor  differences 
between  free  stream  speeds  at  the  nearby  LDV  measurement  locations  arise  from  the 
finite  precision  of  the  downstream  LDV  system  and  the  repeatability  of  the  LCC  flow 
speed  on  different  days.  In  the  next  section,  LDV  results  are  provided  for  the  average 

streamwise  velocity,  U,  the  streamwise  and  wall-normal  velocity  variances,  u  2  and  vn , 
and  the  Reynolds  shear  stress,  -mV.  The  water  temperature  varied  between  18  and  24 
°C  during  these  tests;  its  average  was  21.2°C.  The  density  (p),  and  dynamic  (p)  and 
kinematic  (v)  viscosities  of  water  at  this  temperature  are  p  =  998  km  3,  p  =  0.993X10'3 
kgm'V,  and  v=  9.95xlCT7  mV. 

Unfortunately,  the  size  of  the  two-component  LDV  focal  volume  prevented  its  use 
for  near-wall  velocimetry  measurements  (i.e.  for  y+  =  y/lv  less  than  100  or  so).  Thus,  to 
complete  the  flow  profile  measurements,  three  custom-designed  near-wall  particle- 
tracking-velocimetry  (PTV)  systems  were  assembled  in  water-tight  boxes  and  built  into 
the  test  model  at  x  =  1.96  m,  5.94  m,  and  10.68  m.  The  three  PTV  systems  had  the  same 
design,  though  there  were  slight  variations  among  the  three  systems  during  actual 
implementation.  A  design  schematic  for  one  of  these  systems  is  shown  in  Fig.  2 
delineating  the  illuminating  laser  source,  visualization  camera,  and  optical  pathways. 
Here,  the  near-wall  flow  was  imaged  through  a  periscope  prism  that  protruded  into  the 


flow.  The  imaged  region  was  preferentially  (anisotropically)  magnified  in  the  wall- 
normal  direction  a  using  a  cylindrical  lens  combination  that  produced  a  2-mm-vertical 
field  of  view.  The  geometrical  and  hydrodynamic  design  choices  of  the  prism  were 
modeled  after  the  periscope-mirror  LDV  work  of  Compton  and  Eaton  (1996,  1997).  An 
advantage  of  imaging  through  this  prism  design,  in  addition  to  eliminating  the  need  for 
large  focal  length  optics,  was  the  formation  of  a  reflection  of  the  PTV  particle  images 
from  the  polished  test  surface  that  could  be  used  to  accurately  locate  the  wall  in  the  PTV 
images.  The  three  dual-cavity,  double-pulsed  Nd-YAG  lasers  (Big  Sky  120  mJ,  New 
Wave  Solo  200  mJ,  and  New  Wave  Solo  50  mJ)  were  tightly  bolted  to  the  exterior  side  of 
the  LCC  to  minimize  relative  motion  between  the  lasers  and  the  illumination-path  optics. 
Each  laser’s  green  beam  (532  nm)  was  relayed  from  the  exterior  of  the  tunnel  through  a 
straight  40-mm-diameter  water-tight  conduit  to  the  sheet-making  optics  in  the  test-model- 
interior  water-tight  box  and  then  reflected  into  the  flow  via  a  45°  laser  mirror  to  form  a 
light  sheet  lying  in  constant-z  planes  having  a  thickness  of  75  ±  25  |im.  The  sheet 
thickness  was  measured  directly  by  passing  the  beam  at  the  measurement  location 
through  the  objective  lens  of  a  graded  microscope  and  expanded  upon  exiting  from  the 
eyepiece  onto  a  white  screen.  Titanium  dioxide  seed  particles  with  a  nominal  size  of  1 
|i.m  and  specific  gravity  4.26  from  J.T.  Baker  Industries  were  flood  seeded  in  the  LCC 
water  to  scatter  the  laser  light.  Adding  seven  to  ten  kg  of  Ti02  powder  to  unseeded  tunnel 
water  resulted  in  roughly  5  particles  in  a  32  pixel  by  32  pixel  region  of  the  image.  For 
each  PTV  system,  the  receiving  optical  train  began  with  a  periscope  prism  that  collected 
the  particle-scattered  laser  light  and  directed  it  towards  the  recording  camera  with  the  aid 
of  an  additional  75-mm-diameter  mirror.  Converging  (focal  length  =  +75  mm)  and 
diverging  (focal  length  =  -300  mm)  cylindrical  lenses  were  placed  in  front  of  the  cross¬ 
correlation  CCD  camera  to  provide  a  nominal  l-to-3  stretching  of  the  camera’s  field  of 
view  in  the  vertical  direction  to  increase  the  sensitivity  of  the  PTV  systems  to  vertical 
velocity  fluctuations.  The  PTV  camera  was  fitted  with  a  Nikon  105  mm/f2.8  micro- 
Nikkor  lens  attached  to  two  Nikon  TC-201  teleconverters  (2X  each).  Spatial  calibration 
of  the  resulting  images  was  made  possible  through  imaging  of  a  precision  test  target 
placed  in  the  PTV  interrogation  plane  that  was  composed  of  a  square  array  of  20  |im 
round  dots  with  center-to-center  spacing  of  50±0.02  fim.  For  the  systems  at  x  =  1.96  and 


5.94,  the  resolution  of  the  cross-correlation  CCD  cameras  was  1200  by  1600  pixels 
(LaVision  Imager  Pro),  while  that  at  x  =  10.68  was  1024  by  1280  pixels  (LaVision 
FlowMaster  CCD  camera);  in  all  cases  the  cameras  were  oriented  to  utilize  their  higher 
resolution  in  the  wall-normal  direction.  The  final  x,y-frame-sizes  were  5.03  mm  by  2.7 
mm  at  x  =  1.96  m,  4.88  mm  by  1.97  mm  at  x  =  5.94  m,  and  4.78  mm  by  2.01  mm  at  x  = 
10.68  mm.  Thus,  each  camera  pixel  nominally  corresponded  to  a  rectangle  4.4  pm  in  the 
x-direction  and  1.5  pm  in  the  y-direction  so  that  a  single  particle’s  image  was  typically  an 
ellipse  covering  an  area  of  5  pixels  by  13  pixels.  The  lasers,  cameras,  and  data  recording 
were  controlled  by  PC-type  computers  running  DaVis  7.1  software  by  LaVision  Inc.  For 
the  range  of  investigated  free  stream  speeds,  the  time  separation  between  laser  pulses.  At, 
varied  from  5.5  ps  to  35  ps  with  accuracy  better  than  50  ns  so  that  average  particle 
displacements  were  typically  12  to  16  pixels  in  the  mean  flow  direction.  Although  dual¬ 
frame  acquisition  rates  from  4  to  12  per  second  were  achieved,  the  image  sizes  and  flow 
speeds  ensured  that  any  frame  pair  was  temporally  uncorrelated  from  the  preceding  or 
following  pair.  At  each  flow  condition  3500  image  pairs  were  recorded,  and  the  vectors 
extracted. 

The  PTV  processing  was  conducted  as  a  series  of  steps  that  converted  elongated 
(elliptical)  particle  images  to  velocity  field  averages  and  variances.  The  processing 
began  by  using  the  reflection  from  the  test  surface  to  determine  the  wall  location  within 
each  instantaneous  PTV  image  with  an  accuracy  of  better  than  ±  1  pixel.  This  was  done 
via  a  two-dimensional  cross  correlation  computation  between  the  imaged  particles  and 
their  reflection  pattern.  This  step  was  followed  by  identifying  and  accurately  locating  the 
x,y-pixel-location  centroid  of  individual  particle  images  to  better  than  ±0.1  pixel  with  a 
two-dimensional  spatial  matched-filter  search  of  each  image  frame.  Fortunately,  particles 
very  near  the  wall  typically  produced  images  that  merely  overlapped  the  wall  location 
and  particle-image  centers  could  be  reliably  located  within  a  pixel  or  two  of  the  wall. 
Thus,  the  overall  resolution  of  the  PTV  processing  scheme  was  more  than  an  order  of 
magnitude  better  than  the  5-pixel-by- 13-pixel  average  particle-image  size  for  particles 
with  image  centers  more  than  one  pixel  above  the  wall.  The  matched-filter-determined 


particle-image-center  locations  were  converted  back  to  physical  .^-coordinates  using  the 
information  from  the  calibration  image. 

The  primary  challenge  in  PTV  is  to  correctly  identify  the  image  from  the  same 
particle  in  both  frames.  Here,  this  task  was  accomplished  by  selecting  a  particle  image  in 
the  first  frame  and  then  searching  a  spatial  window  in  the  second  illumination  frame 
where  this  particle’s  second  image  was  expected  to  appear.  If  no  particles  were  found  in 
the  second-frame  search  window,  that  particle  was  discarded.  If  one  or  more  particle 
images  were  found  in  the  second-frame  search  window,  then  two  concurrent  tests  were 
conducted  to  determine  the  highest  probability  match.  First,  a  two-dimensional  cross 
correlation  value  was  computed  between  the  particle  intensity  distribution  in  the  first 
frame  and  all  possible  particle  images  in  the  second  frame.  The  correlation  window  was 
24  pixels  x  24  pixels  centered  at  the  particle  images’  centroidal  pixels.  The  second-frame 
image  with  the  highest  cross-correlation  value  was  most  likely  to  be  the  correct  match  to 
the  first-frame  particle.  The  second  test  involved  checking  whether  the  highest  computed 
correlation  value  was  an  actual  correlation  peak,  by  repeating  the  correlation  calculation 
eight  times  corresponding  to  eight  different  shifts  of  the  first  242-pixel  box  around  the 
second-frame  particle  image  centroid.  If  the  correlation  value  at  zero  shift  was  the 
highest,  then  the  image  was  considered  the  correct  match  for  the  particle  in  the  first 
frame.  The  local  particle  velocity  (m,v)  was  then  simply  determined  from  the  particle’s 
centroid  displacement  divided  by  time  between  frames,  At.  The  data  processing  then 
moved  to  the  next  first-frame  particle  image  and  the  above  process  was  repeated  for  all 
possible  first-frame  particles.  At  the  end  of  this  process,  each  frame  pair  produced  an 
instantaneous  velocity  field  that  was  randomly  sampled  in  space. 

Next,  a  vector  post-processing  filter  was  passed  over  each  instantaneous  field  to 
eliminate  spurious  vectors.  The  filter  compared  the  vertical  and  horizontal  components  of 
the  vector  under  consideration  with  the  mean  and  root  mean  square  components  of  all  the 
vectors  within  a  radius  of  80  pixels.  For  the  measurements  at  x  =  1.96  m  and  10.68  m, 
the  spurious  vector  percentage  ranged  between  4%  and  7%.  At  x  =  5.94  m,  the  spurious 
vector  percentage  was  higher,  10%  to  15%,  because  of  non-smooth  particle  images 
possibly  caused  by  optical  distortion.  These  instantaneous  velocity  field  results  were  then 
collapsed  in  the  streamwise  direction  and  across  all  frame  pairs  from  the  same 


experimental  run  by  sorting  the  instantaneous  velocity  vectors  into  vertical  bins.  Small 
bins  were  chosen  close  to  the  wall  (ranging  in  size  from  1.5  to  4  pixels  depending  on  the 
flow  speed)  to  resolve  the  near  wall  flow,  while  far  away  from  the  wall  y*  >  50,  larger 
bins  of  16  pixels  were  used.  Once  compiled,  there  were  typically  10  thousand  vector 
samples  in  the  small  bins  nearest  the  wall,  and  50  thousand  vector  samples  in  the  larger 
bins  farther  from  the  wall.  Averages  and  variances  for  each  bin  were  computed  directly 
from  the  vector  samples  in  that  bin  and  reported  at  the  bin-center  vertical  location. 

For  PTV  to  be  successful,  the  second-frame  search  window  size  and  shape  needs 
to  be  optimized.  When  the  second-frame  search  window  is  too  small  or  poorly  placed, 
correct  second-frame  particle  images  are  less  likely  to  be  found  and  potentially  useful 
flow  information  is  lost.  However,  when  the  second-frame  search  window  is  too  large, 
the  likelihood  of  two  or  more  ambiguous  second  frame  images  increases,  and  the 
probability  of  including  spurious  vectors  increases.  In  both  cases,  averages  and  variances 
may  be  statistically  under-resolved  or  even  biased  by  a  poor  choice  of  the  second-frame 
search  window.  For  the  PTV  results  presented  here,  the  shape  of  second-frame  window 

is  a  based  on  a  power  law  of  the  form  (y-y0)/Y  =  ±((;t-;t0)/X)1/6  where  ( xoty0 )  is  the 
first-frame  particle-image-centroid  coordinate,  distances  are  measured  in  pixel  units,  and 
0  <  x  -  x0  <  X  and  \y  -  ye|  <  Y.  Figure  3a)  shows  this  search-window  shape.  The  size  of 
the  search  envelope  was  allowed  to  vary  with  wall-normal  distance  to  account  for  the 
strong  velocity  gradient  of  the  flow.  Near  the  wall,  y+  <  15,  the  values  of  X  and  Y 
increased  linearly  and  quadratically  with  y,  respectively;  according  to 
X  =  a,  +a2y(rw/n)At  and  Y  =  a3  + a4y2(rw//j)At ,  where  the  a's  are  empirical  constants: 
a,  and  a3  have  units  of  pixels,  is  unitless,  and  a4  has  the  units  of  inverse  pixels.  The 
other  parameters  ( zw ,  pt,  and  At)  have  SI  units  and,  in  the  near-wall  region,  the  grouping 
of  parameters,  y(rw//j.)At  =  y(dU/dy)At  =  UAt,  represents  the  average  downstream 
particle  displacement  in  the  time  interval  At.  A  slight  distance  above  the  wall  (y*  ~  15)  X 
and  Y  reach  constant  values  of  X ^  (=  23  pixels)  and  Y ^  (=  13  pixels),  respectively. 


These  maximum  values  were  determined  to  ensure  that 


and 


Y  >  throughout  the  measurement  domain.  Here,  a,  =  X  and  a3  =  Y  at  the 


wall  (y  =  0)-  The  factor  controls  how  fast  the  envelope’s  streamwise  extent  grows  with 
increasing  y  (optimal:  0.50  to  0.56);  and  a4  controls  how  fast  the  envelope’s  vertical 
extent  grows  with  increasing  y  (optimal:  0.006  to  0.009).  Figure  3b)  shows  how  X  and  Y 
were  varied  for  increasing  y+. 

Figure  4  illustrates  the  PTV  processing  method.  Figure  4a)  shows  first-frame 
particle  images,  with  the  matched-filter-identified  particle  images  marked  at  their 
centroids.  The  diagonal  extent  of  the  particle  images  is  a  remnant  of  minor  misalignment 
in  the  PTV-system  collection  optics.  Figure  4b)  shows  the  identified  first-frame  particles 
centroids  from  Fig.  4a)  along  with  the  search  window  where  second  frame  particles  are 
expected  to  be  found.  The  identified  second-frame  particles  are  also  shown.  The  choice 
of  the  search  window  parameters  for  this  case  (a,  =  1.5  pixel,  Oj  =  0.50,  a3  =  1.25  pixel,  a4 
=  0.0075  pixel'1)  places  second-frame  particles  well  inside  the  search  envelope.  Some  of 
the  first-frame  particles  have  no  matching  second  frame  particles  and  vice  versa  because 
spanwise  particle  motion  may  remove  particles  from  the  laser  sheet  between  laser  pulses, 
and  thresholding  of  the  matched-filter  output  causes  some  weak  particle  images  to  be 
discarded.  For  the  PTV  results  reported  here,  the  typical  selection  of  the  a’s  was  the 
same  as  above  for  a,,  Oj  and  a3,  with  aA  =  0.006  pixel'1. 

To  assess  the  sensitivity  of  the  PTV  measurements  to  the  a’s,  a  case  study  is 
presented  in  Fig.  5,  where  six  different  search  parameter  variations  are  used  determine 

the  profiles  of  u 2  and  v'2 ,  the  turbulence  quantities  found  to  be  most  sensitive  to  choice 
of  the  a’s.  The  factor  that  controls  how  fast  the  window’s  streamwise  extent  grows 
with  vertical  distance  is  varied  from  a  value  of  0.41  (x)  to  0.66  (o);  and  a4  that  controls 
how  fast  the  window’s  vertical  extent  grows  is  varied  from  a  value  of  0.0045  (+)  to  0.012 
(o).  When  the  a’s  are  significantly  smaller  than  the  typical  values  specified  above,  some 
of  the  particle-image  data  can  be  filtered  out  resulting  in  lowered  turbulence  quantities  (+, 
x),  and  the  opposite  is  true  (o)  due  to  an  increased  probability  of  the  occurrence  of 
outliers  with  oversized  envelopes.  For  the  other  three  cases  (*,  V,  •)  [0.50  <  a2<  0.56; 
0.0075  <  a4<  0.012]  the  measured  structure  of  the  boundary  layer  is  nearly  insensitive  to 
the  processing  parameters.  Overall,  for  these  six  cases  that  cover  a  wide  variation  in 
search  parameters,  the  processing-parameter-induced  changes  in  the  measured  profiles 
are  mild  and  are  of  the  same  size  as  the  experimental  uncertainty.  Thus,  the  dependence 


of  the  PTV  measurements  on  the  PTV  processing  parameters  was  not  considered 
significant,  and  all  PTV  results  presented  are  determined  from  the  typical  a-values. 

3.0  Average  Flow  Results  (plot  composite  PTV-LDV  profiles  where  possible) 

3.1  Static  Pressure 

The  dependence  of  the  average  static  pressure  in  the  LCC  test  section  on  the 
streamwise  coordinate  x  is  shown  in  Fig.  5  in  terms  of  the  pressure  coefficient,  CP(x). 
Downstream  of  the  test  model’s  elliptical  leading  edge,  the  tunnel  and  test-model 
hardware  cross  sections  were  geometrically  constant  for  more  11  m.  Thus,  boundary 
layer  growth  on  the  test  model  and  the  LCC-test-section  sidewalls  is  the  primary  cause  of 
the  slightly  favorable  pressure  gradient  that  is  observed  at  all  three  primary  test  speeds. 

Although  easily  measured  and  non-zero,  the  observed  pressure  gradient  is  small 
enough  so  the  measured  boundary  layer  characteristics  should  approach  those  of  a  truly 
zero  pressure  gradient  flow.  At  the  three  primary  test  speeds,  the  acceleration  parameter, 

K  =  (v/U'}(dUe/dx),  is  less  than  10"9;  a  substantial  deviation  from  the  classical  log-law 

only  occurs  for  AT>1.6xlO“6  (Patel  1965).  Furthermore,  when  the  Reynolds-averaged 
streamwise  boundary  layer  equation  is  rendered  dimensionless,  the  near-wall  region 
where  the  momentum  flux  terms  are  small  (y+  ^  30)  will  not  be  influenced  by  a 
streamwise  pressure  gradient  when  Ay+«1  where  A  =  (v/ pu*3}(dP/dx)  and 

u*  =  ^Tlv/p  (see  Afzal,  2001).  For  these  experiments,  A  was  of  order  10~5  or  smaller. 
Similar  considerations  suggest  that  the  outer  region  where  the  viscous  stress  is  small  (y + 
£  100)  will  not  be  influenced  by  the  streamwise  pressure  gradient  when 

Cl=(dg9/Tw)(dP/dx)«\.  However,  in  these  experiments,  Q  ranged  from  0.1  to  0.2,  so 

some  differences  are  expected  between  the  measured  outer  layer  flow  and  a  truly  zero- 
pressure-gradient  TBL  at  the  same  Reynolds  number. 


3.2  Average  Skin  Friction 


The  variation  of  the  time-average  skin  friction  on  the  test  surface  is  shown  in  Fig. 
6  in  terms  of  the  skin  friction  coefficient,  and  the  downstream-distance-based 
Reynolds  number,  Rex  =  Uex/v,  calculated  from  the  local  value  of  Ue.  Results  from  1.96 
10.68  m,  and  Ue  =  3.3,  6.65,  13.2,  and  19.9  m/s  are  plotted  on  this  figure  with 
error  bars  representing  the  95%  confidence  interval.  Three  separate  techniques  were  used 
to  determine  Cf  and  these  are  represented  by  different  symbols  in  Fig.  6.  The  open 
symbols  (O,  □,  A)  signify  direct  measurements  of  rw  made  with  the  force  balances.  The 
closed  symbols  (#,  ■,  ▲)  represent  measurements  of  rw  =  p(dU/dy)y=0  determined  from 

the  near-wall  PTV  measurements  of  U(y )  via  a  linear  least-square  fit  for  0.5  ^  y*  ^  5.5 
using  U  =  0  at  y  =  0,  and  the  line-only  symbols  (+,  x,  V)  symbols  represent  measurements 
of  the  total  stress  rw  ~  p(dU/dy)~  pu'v'  as  determined  from  the  LDV  measurements  in 
the  TBL’s  constant-stress  region,  100  ^  y*  ^  2000.  The  solid  line  is  a  simple  least- 
squares  power-law  fit  to  the  force-balance  measurements. 

As  is  readily  apparent,  the  three  separate  techniques  for  determining  rw  all  agree 
to  within  experimental  accuracy  and  they  are  all  well  represented  by  the  power  law  fit, 

Cf  =  (0.0 1 70  ±  0.0004)  •  Re~x(0 ,237±00067).  (3.1) 

Here  it  must  be  noted  that  the  empirical  constants  in  this  law  may  implicitly  include  the 
influence  of  flow  artifacts  arising  from  the  boundary  layer  trip  and  the  test  model’s 
elliptical  leading  edge.  Yet,  such  influences  should  be  small.  Together  the  distributed 
roughness  trip  and  curved  portion  of  the  test  model’s  leading  edge  occupy  the  first  0.37  m 
of  the  test  model,  only  3.5%  of  the  wetted  length  (10.68  m)  of  the  experiment.  Given  the 
moderate  scatter  in  the  measurements  and  the  relative  success  of  the  power  law  fit,  the 
values  of  Ct  (and  Tw)  used  for  all  subsequent  evaluations  and  normalizations  are 
determined  from  Eq  (3.1). 

3.3  Average  Streamwise  Velocity  Profile 

•  completeness  of  measurements;  from  0. 1  ^  y+  ^  10A5 

linear,  log-law,  and  wake  regions  -  overlap  between  instruments 
good  near-wall  measurements 

as  much  as  2-decades  of  y+  in  the  classical  “log-law”  region 


•  fitted  values  of  “kappa”  and  “B”  across  all  the  current  data  are _ ± _ 

and _ ± _ and  these  are  in  good  agreement  with  prior  values  for  TBLs 

•  wake  component  results  are  consistent  with  that  of  prior  experiments 

values  for  the  wake  parameter  are  “close” 

current  values  are  likely  lower  because  of  mild  acceleration  of  the 

flow 

•  log-law  vs.  power-law  checks  are  uncertain,  but  the  current  data  when 
analyzed  via  fractional  differences  suggests  that  the  log-law  is  a  better  fit  over  the  chosen 
range  of  y. 

3.4  Discrete  Evaluation  of  the  Von  Karman  Boundary  Layer  Integral  Equation 

As  a  final  overall  check  of  the  LDV  measurements,  the  skin  friction  on  the  plate 
midway  between  the  middle  and  downstream  measurement  stations  was  estimated  from  a 
discrete  evaluation  of  the  Von  Karman  boundary  layer  integral  equation, 


|V)  (t/ft-t/fr)  |  (32) 

pL~  **  2  to  ’ 

where  x  =  (x,  +  x2)/2  =  8.3 1  m,  Ax  =  x2  -  jc,  =  4.74  m,  and  the  subscripts  1  and  2  refer  to 
average  conditions  at  xx  =  5.94  m  and  x2  =  10.68  m,  respectively.  The  momentum  (0)  and 
displacement  (S')  thicknesses  were  calculated  from  discrete  evaluations  of  their  usual 

definitions,  6=  i*(U/Ut)(l-U/Ue)dy  and  <5*  =  J*  (1  -U/Ue)dy,  with:  i)  C/  =  0aty  =  0, 

ii)  h  defined  to  be  the  vertical  location  from  the  test  surface  where  the  LDV-measured 
Reynolds  shear  stress  dropped  to  the  free-stream  background  noise  level,  and  Hi)  Ue 
defined  to  be  the  average  of  U  for  first  three  measurement  points  above  y  =  h.  Individual 
values  of  0and  5  are  given  in  Table  1.  For  evaluation  of  Eq.  (3.2),  the  values  of  (/,  were 
drawn  directly  from  the  LDV  measurements,  while  the  values  of  U2  were  consistent  with 
the  LDV  results  but  were  slightly  corrected  (a  0  to  0.8%  velocity  increase)  to  match  the 
measured  CP  drop  between  jc,  and  ;t2.  Here,  the  momentum  thickness  difference  in  the 


first  term  on  the  right  side  of  Eq.  (3.2)  provides  the  dominant  contribution  because  Ux  and 
U2  differ  by  only  a  percent  or  so. 

The  results  of  this  estimation  exercise  showed  that,  at  the  three  primary  test 
speeds,  the  inferred  skin  friction  results  from  Eq.  (3.2)  were  35%,  31%,  and  32%  lower 
than  the  power-law  fit  to  the  direct  skin  friction  measurements  provided  by  Eq.  (3.1). 
The  consistency  of  the  mismatch  across  test  speeds  suggests  that  the  mismatch  was 
caused  by  a  bias  problem  rather  than  random  error.  The  most  likely  source  for  such  a 
bias  error  is  the  unintentional  optical  distortion  in  the  LDV  windows  used  at  x  =  5.94  m 
and  x  =  10.  68  m.  Such  optical  distortion  leads  to  minor  variations  in  the  LDV  beam  foci 
location  and  beam-crossing  angle,  and  can  cause  spatially-varying  bias  error  in  LDV 
average  velocities  of  approximately  ±0.5%  to  ±1%.  While  such  bias  errors  may  be 
negligible  for  assessing  profile  shapes,  they  can  have  an  order  of  magnitude  greater 
influence  (±5%  to  ±  10%)  on  6  (and  <5*).  And,  when  9  (and  (S’)  are  assigned  such 
uncertainties,  the  range  of  possible  skin  friction  values  from  Eq.  (3.2)  overlaps  the  Eq. 
(3.1)  results.  Thus,  the  inferred-vs.-measured  skin-friction  mismatch  was  concluded  to 
result  from  the  optical  imperfections  in  the  LDV  windows,  experimental  artifacts  that 
(regrettably)  were  not  avoided. 

4.0  Second-order  Turbulence  Statistics  (plot  composite  PTV-LDV  profiles  where 
possible) 

inner  variable  collapse  is  imperfect  -  higher  peak  u'2  as  Re  increases 
u 2  at  specific  points  does  follow  prior  proposed  scaling  laws 
D&E’s  mixed  scaling  of  u'2  appears  superior  to  inner  variable  scaling 

•  V2 

inner  variable  collapse  is  OK? 
imperfect  overlap  between 

•  -mV  +  total  TBL  stress 

shear  stress  collapses  well  for  the  inner  layers? 

total  stress  is  consistent  with  the  force-balance  shear  stress 


5.0  Summary  and  Conclusions 

•  Unique  High  Re  Data  Set  -  precisely  the  type  of  info  that  should  be  archived! 

•  (possibly  the  last  exp.  of  this  type  for  sometime  to  come) 

•  largest  possible  flat  plate  model  that  would  fit  in  the  worlds  largest  water  tunnel 

•  Results  fall  between  prior  controlled  lab  tests  and  other  High  RE  experiments 

•  Well-controlled  experiment 

•  low-turbulence  water  tunnel,  small  variation  in  speed  &  temp. 

•  polished  flat  test  surface 

•  four  instruments  providing  measurements  that  are  consistent  within  error 

•  (force  balances,  static  pressure,  LDV,  PTV). 

•  Tabulated  results  for  the  community  to  use 

•  the  present  preliminary  analysis  suggests  the  log-law  to  be  better  than  a  power 
law 

•  u2  increases  according  to  the  proposed  scaling  relationships 

•  other  turbulence  quantities  are  “well  behaved” 
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IV.  TABLES 


Table  1 :  TBL  Characteristics  for  the  LDV  profile  measurement  conditions. 
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Fig.  1.  Schematic  of  the  test  model  showing  the  test  surface  (a)  and  the  instrumentation 
access  side  (b).  (ESWT  Fig.  3) 

Fig.  2.  Schematic  of  the  near  wall  particle  tracking  velocimetry  setup.  The  illumination 
light  pulses  originate  in  a  double-cavity  Nd-YAG  PIV  laser  head  that  is  mounted 
to  the  outside  wall  of  the  water  tunnel.  Sheet-forming  optics  and  a  45°  mirror 
inside  a  water-proof  optics  box  relay  the  fan  the  light  into  a  sheet  that  is  parallel 
the  flow  direction  and  perpendicular  to  the  test  surface.  Side-scattered  light  from 
the  sheet  is  directed  via  another  45°  mirror  to  a  digital  camera  after  passing 
through  two  cylindrical  lenses  (focal  lengths:  +75  mm,  and  -300  mm)  to  provide 
anisotropic,  preferential  magnification  in  the  y-direction. 

Fig.  3.  PTV  second-frame  search  window  shape  (a),  and  size  specifications  as  a  function 
of  y*  (b).  The  search  window  is  X-pixels  (solid  line)  long  with  a  vertical  extent  of 
±Y  pixels  (dash-dotted  line).  Here,  X  and  Y  increase  with  the  distance  from  the 
wall  until  y+  ~  15  where  they  become  constant. 

Fig.  4.  Sample  PTV  processing  images,  (a)  First-frame  particle  images  with  the  matched- 
filter-identified  particle  centroids  marked  with  a  white  cross  in  a  black  circle.  The 
wall  is  located  at  y+  =  0  and  the  flow  is  from  left  to  right.  The  frame  is  125  pixels 
in  the  x-direction  by  100  pixels  in  the  y-direction).  (b)  A  composite  drawing 
superimposing  the  second-frame  matched-filter-identified  particle-image  centers 
(black  dots),  the  particle  image  locations  from  the  first  frame  (white  cross  in  a 
black  circle),  and  search  envelopes,  (c)  The  vectors  resulting  from  matching  the 
particles  in  the  dual-frame  exposures. 

Fig.  5.  PTV-determined  values  of  the  normalized  velocity  variances:  punjxw  =  un  (a), 

and  pv,1lxw  =  v,2+  (b)  for  six  different  combinations  of  the  second-frame  search 
window  size  parameters  for  Ue  =  13.2  m/s  and  x  =  10.68  m:  (x)  ax  =  1.4,  a2  = 

0.41,  a3  =  1.25,  a4  =  0.006;  (+)  ax  =  1.5,  a,  =0.46,  a,  =  1.4,  a4  =  0.0045;  (*)  ax  = 
1.65,  a2  =0.50,  a3  =  1.5,  a4  =  0.0075;  (V)  a,  =  1.9,  a2  =0.53,  a,  =  1.6,  a4  =  0.01 1; 
(•)  a,  =  1.8,  Oj  =0.56,  o3  =  1.6,  a4  =  0.009;  (o)  a,  =  2.2,  a2  =0.66,  a3  =  1.9,  a4  = 
0.012.  The  search  envelope  maximum  sizes  (X^,  Y „,)  in  pixels  are  (21,  1 1)  for 
(x)  and  (+),  (22,  12)  for  (*)  and  (•),  and  (23,  13)  for  (V)  and  (o). 


Fig.  6.  Static  pressure  coefficient,  CP(x)  =  (P(x)~  Px)/^pUl,  vs.  distance  from  the 

leading  edge  of  the  test  model,  x  (in  m),  for  the  three  primary  test  speeds:  U„= 
6.65  m/s  □  [open  square],  13.2  m/s  O  [open  circle],  and  19.9  m/s  A  [open 
triangle].  Here,  Px  is  the  static  pressure  measured  at  the  first  tap  at  x  =  1.96  m. 
Boundary  layer  growth  on  the  model  and  the  LCC  test  section  sidewalls  lead  to  a 
mildly  favorable  pressure  gradient.  (ESWT  Fig.  6) 


Fig.  7.  Friction  Coefficient  vs.  Rex.  Summary  plot  of  all  shear  stress  measurements:  force 
balances,  near-wall  PTV,  LDV  total  stress. 

(ESWT  Fig.  32) 

Fig.  8.  U  vs.  y  -  inner  variables,  0.1  ^  y+  <  105.  Show  linear  fit  for  y+  <  10,  and  log-law 
fit  for  y+  >  10  shown.  (ESWT  Fig.  8)  This  figure  should  be  big  and  plotted  in  the 
landscape  orientation. 

Fig.  9.  Wake  parameter,  A U/ur ,  vs.  Reg  =  Ue6/v,  from  the  LDV  measurements  at  x  = 
5.94  m  A  [open  triangle]  and  x  =  10.68  m  O  [open  diamond].  Also  shown  are 
the  results  of  DeGraaff  and  Eaton  (2000)  □  [open  square]  for  log-law  constants 
of  K=  0.41  and  B  =  5.0.  The  solid  curve  is  the  correlation  from  Coles  (1962) 
which  asymptotes  to  approximately  A U/ux  =2.7  at  high  Reynolds  number.  The 
current  data  likely  fall  below  2.7  because  of  the  mild  favorable  pressure  gradient 
in  the  LCC  test  section.  (ESWT  Fig.  14) 

Fig.  10.  Comparison  of  the  diagnostic  mean-profile  fitting  functions  from  Osterlund 

(2000)  for  the  log-law,  *F  =  (y*  dU+ / dy+}  '  shown  in  open  symbols,  and  power- 
law,  T  =  (y+ / U+MdU+ / dy+}  shown  in  solid  symbols.  Part  a)  shows  *F  and  T  vs. 

y+  for  x  =  5.94  m  and  {/„=  13.2  m/s  (squares),  x  =  5.94  m  and  U„=  19.9  m/s 
(triangles),  x  =  10.68  m  and  (/„=  13.2  m/s  (circles),  and  x  =  10.68  m  and  t/„= 
19.9  m/s  (diamonds).  Part  b)  shows  the  mean  fractional  difference  of  *F  and  T  vs. 
Ree  =  Ue 6/ v  for  y+  >  150  and  y  <  0. 1 5  5*,-  The  lower  mean  fractional  difference 
for  the  log  law  (open  symbols)  suggests  that  it  is  a  better  fit  to  the  experimental 
data  than  a  power  law  (solid  symbols)  in  the  stated  range  for  y.  (ESWT  Fig.  15  a) 
and  b)),  and  (ESWT  Fig.  17) 

Fig.  11.  u’A2  vs.  y,  inner  variables.  (ESWT  Fig.  18) 

Fig.  12.  part  a)  (u’A2)max  vs.  ReT.  (ESWT  Fig.  19) 

part  b)  u’A2  at  y+  =  200  vs.  Rev  (ESWT  Fig.  20) 

Fig.  13.  “Eaton  scaling”  of  u’A2  vs.  y,  inner  variables.  (ESWT  Fig.  21) 

Fig.  14.  v’A2  vs.  y,  inner  variables.  (ESWT  Fig.  22  &  23) 

Fig.  15.  uV  vs.  y,  inner  variables.  (ESWT  Fig.  24  &  25) 

Fig.  16.  Normalized  total  TBL  shear  stress  vs.  y,  inner  variables,  log  horizontal  axis 
part  a)  6.65  m/s  (all  three  measuring  stations  ?) 
part  b)  13.2  m/s  (all  three  measuring  stations  ?) 
part  c)  19.9  m/s  (all  three  measuring  stations  ?) 

(composite  or  summary  of  ESWT  Figs.  26-31) 
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Fig.  5.  Coefficient  of  pressure  vs.  downstream  distance  (ESWT  Fig.  6) 
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figure  should  be  big  and  plotted  in  the  landscape  orientation. 
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Fig.  8.  Wake  parameter,  A U/ut ,  vs.  Reg  =  Ued/v ,  from  the  LDV  measurements  at  x  = 
5.94  m  A  [open  triangle]  and  x  =  10.68  m  O  [open  diamond].  Also  shown  are 
the  results  of  DeGraaff  and  Eaton  (2000)  □  [open  square]  for  log-law  constants 
of  K—  0.41  and  B  =  5.0.  The  solid  curve  is  the  correlation  from  Coles  (1962) 
which  asymptotes  to  approximately  At//wT  =2.7  at  high  Reynolds  number.  The 
current  data  likely  fall  below  2.7  because  of  the  mild  favorable  pressure  gradient 
in  the  LCC  test  section. 
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Fig.  9.  Comparison  of  the  diagnostic  mean-profile  fitting  functions  from  Osterlund 

(2000)  for  the  log-law,  *F  =  (y+  dlf/dy*}  shown  in  open  symbols,  and  power- 
law,  T  =  (y+ / U+}(dU*  /  dy*}  shown  in  solid  symbols.  Part  a)  shows  'P  and  T  vs. 

y*  for  x  =  5.94  m  and  U„=  13.2  m/s  (squares),  x  =  5.94  m  and  (/„=  19.9  m/s 
(triangles),  x  =  10.68  m  and  U„=  13.2  m/s  (circles),  and  x  =  10.68  m  and  U„= 
19.9  m/s  (diamonds).  Part  b)  shows  the  mean  fractional  difference  of 'P  and  Y  vs. 
Reg  =Ue6/v  fory+>  150  and  y  <  0.15<^<>.  The  lower  mean  fractional  difference 
for  the  log  law  (open  symbols)  suggests  that  it  is  a  better  fit  to  the  experimental 
data  than  a  power  law  (solid  symbols)  in  the  stated  range  for  y.  (ESWT  Fig.  15  a) 
and  b)),  and  (ESWT  Fig.  17) 
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Fig.  10.  u’A2  vs.  y,  inner  variables.  (ESWT  Fig.  18) 
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Fig.  1 1.  part  a)  (u’A2)max  vs.  Rer  (ESWT  Fig.  19) 

part  b)  u’A2  at  y+  =  200  vs.  ReT.  (ESWT  Fig.  20) 
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Fig.  12.  “Eaton  scaling”  of  u’A2  vs.  y,  inner  variables.  (ESWT  Fig.  21) 
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Fig.  13.  v’A2  vs.  y,  inner  variables.  (ESWT  Fig.  22  &  23) 
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Fig.  14.  u’v’  vs.  y,  inner  variables.  (ESWT  Fig.  24  &  25) 


Fig.  15.  Normalized  total  TBL  shear  stress  vs.  y,  inner  variables,  (composite  or  summary 
of  ESWT  Figs.  26-31) 
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A  set  of  experiments  has  been  conducted  at  the  William  B.  Morgan  Large  Cavitation  Channel  (LCC)  in 
Memphis,  Tennessee  to  investigate  the  phenomena  of  friction  drag  reduction  at  large  scales  and  high 
Reynolds  numbers  in  a  controlled  test  environment.  Drag  reduction  was  achieved  by  injecting  aqueous 
polyethylene  oxide)  (PEO)  solutions  at  three  different  mean  molecular  weights  into  a  near-zero-pressure- 
gradient  turbulent  boundary  layer  formed  on  a  smooth  flat  test  surface  having  a  length  of  nearly  1 1  m.  The 
test  model  had  an  overall  length  of  12.9  m,  and  spanned  the  3.05  m  width  of  the  LCC  test  section.  Friction 
drag  was  measured  with  six  floating-plate  force  balances  at  downstream-distance-based  Reynolds  numbers 
up  to  220  million  and  free  stream  speeds  up  to  20  ms 1 .  For  a  given  polymer  type,  the  level  of  drag 
reduction  was  measured  for  a  range  of  free  stream  speeds,  polymer  injection  rates,  and  concentrations  of 
the  injected  solution.  Polymer  concentration  fields  in  the  near-wall  region  (0  <  y  <  2  mm)  were  examined 
at  three  locations  downstream  of  the  injector  using  near-wall  pulsed  laser-induced-fluorescence  imaging. 
Comparisons  of  the  development  and  extent  of  friction  drag  reduction  are  drawn  between  the  various  test 
results  from  the  current  study,  and  between  the  current  and  prior  PEO-solution-based  friction  drag 
reduction  experiments.  The  evolution  of  near  wall  polymer  concentration  is  divided  into  three  regimes:  the 
development  region  near  the  injector,  where  the  injected  polymer  is  highly  stratified,  forming  filaments  in 


the  near-wall  flow;  the  transitional  mixing  region,  where  the  polymer  begins  to  mix  across  the  boundary 
layer;  and  the  final  region,  where  the  polymer  mixing  occurs  at  the  rate  of  boundary  layer  growth.  The 
observed  drag  reduction  and  mixing  processes  are  compared  to  previously  reported  results,  using  the 
standard  K-Factor  scaling,  the  physical  basis  and  usefulness  of  this  scaling  is  discussed.  Finally,  the  near 
wall  polymer  concentration  is  correlated  to  the  drag  reduction  achieved  for  the  three  polymer  molecular 
weights  employed.  These  data  are  related  to  the  polymer  solution  rheology  the  local  Weissenberg  number, 
and  the  Weissenberg  numbers  for  drag-reduction  onset  and  shear  based  degradation.  The  data  are  shown  to 
correlate  to  the  universal  drag  reduction  curves  first  proposed  by  Virk  et  al.  (1967)  for  pipe  flows. 
However,  a  scaling  which  would  relate  the  differences  in  the  intrinsic  drag  reduction  and  concentration 
with  changes  in  the  molecular  weight  was  not  found,  possibly  de  to  the  presence  of  significant  flow 
induced  polymer  degradation. 

1.  Introduction 

Reduction  of  turbulent  boundary  layer  skin  friction  in  external  flows  is  an  enduring  priority  for  fuel  savings 
in  marine  transportation  systems.  Researchers  have  known  for  more  than  fifty  years  that  the  presence  of 
high  molecular  weight  polymers  in  the  near-wall  region  of  a  turbulent  boundary  layer  (TBL)  flow  can  lead 
to  reduction  of  the  friction  drag  by  as  much  as  70%  compared  to  the  friction  produced  by  a  corresponding 
flow  of  pure  water,  even  when  the  polymer  concentration  is  on  the  order  of  a  few  tens  of  weight  parts  per 
million  ( wppm ).  The  mechanism  underlying  this  polymer  drag  reduction  has  been  under  investigation  for 
several  decades,  and  an  overall  picture  of  the  interaction  of  the  long-chain  polymers  with  the  underlying 
turbulent  flow  has  emerged  through  analytical,  experimental,  and  numerical  studies.  Dubief  et  al.  (2004) 
discusses  in  detail  how  the  stretching  and  recoiling  of  polymer  molecules  in  the  near  wall  flow  can  modify 
the  self-sustained  near-wall  turbulence  regeneration  cycle.  Near  wall  turbulence  in  Newtonian  flow  is 
sustained  through  the  production  of  near  wall  vortices  in  the  buffer  region  and  streaks  in  the  viscous  sub¬ 
layer  (Jimenez  and  Pinelli,  1999).  The  presence  of  the  polymer  molecules  modifies  the  non-linear 
interactions  by  extracting  energy  from  the  buffer  layer  vortices  and  releasing  it  in  the  near  wall  streaks. 
This  leads  to  enhanced  stream-wise  velocity  fluctuations,  reduced  wall-normal  velocity  fluctuations,  and 
increased  spacing  between  stream-wise  vortices  in  the  near-wall  region.  The  balance  of  stresses  in  the 


flow  near  the  wall  are  then  composed  of  a  combination  of  viscous,  Reynolds,  and  (non-trivial)  polymer 
stresses.  Finally,  the  state  of  maximum  drag  reduction  (MDR)  is  achieved  when  there  is  a  sufficient 
concentration  of  polymer  molecules  in  the  near-wall  region  to  produce  a  new,  self-sustained  regeneration 
cycle.  Once  MDR  is  achieved,  further  increases  in  polymer  concentration  will  no  longer  directly  result  in 
turbulence  modification  and  further  reduction  in  friction  drag.  Instead,  increasing  polymer  concentration  in 
an  MDR  flow  may  increase  skin  friction  through  an  increase  in  the  bulk  shear  viscosity.  The  detailed 
interactions  of  the  polymer  with  the  underlying  turbulent  flow  is  still  the  object  of  study,  but  this 
phenomenological  description  is  consistent  with  experimental  observations  of  drag  reduction,  including  the 
recent  work  of  Warholic  et  al.  (1999),  Ptasinski  et  al.  (2003),  and  White  et  al.  (2004). 

Polymer  additives  are  used  commercially  to  reduce  the  friction  developed  by  flows  in  pipes  and  other 
conduits  (Sellin,  1982).  In  most  internal  flow  applications,  the  polymer  is  introduced  and  mixed  into  the 
flow  to  produce  a  relatively  uniform  and  sustained  polymer  distribution  in  the  flow  (in  the  absence  of 
polymer  degradation).  However,  application  of  polymer  drag-reduction  to  external  flows  is  complicated  by 
the  need  to  inject  the  polymer  solution  and  allow  it  to  be  mixed  with  free-stream  fluid  in  the  TBL  flow. 
Even  at  MDR,  mixing  in  the  TBL  continues.  Without  continual  addition  of  polymer,  the  near-wall  polymer 
concentration  decreases  as  the  TBL  entrains  free-stream  fluid.  Once  the  near-wall  polymer  concentration  is 
sufficiently  reduced,  skin-friction  drag  reduction  is  diminished  and  then  lost.  Therefore,  practical 
application  of  polymer  drag  reduction  to  external  flows  requires  a  careful  balance  between  the  cost  of 
carrying,  mixing,  and  injecting  the  polymer,  and  the  benefit  of  the  resulting  drag  reduction. 

Experimental  investigations  of  external  TBL  drag-reduction  by  polymer  additives  have  been  conducted  at  a 
number  of  scales  and  flow  speeds.  A  variety  of  polymers  having  different  mean  molecular  weighs,  Mw, 
have  been  used  for  TBL  drag  reduction  studies,  with  poly(ethylene  oxide)  (PEO)  and  polyacrylamide 
(PAM)  being  the  two  most  commonly  reported.  Polymer  skin-friction  drag  reduction  has  been  examined 
on  flat  plate  TBL  at  downstream-distance-based  Reynolds  numbers  up  to  45  million  (Vdovin  and 
Smol’yakov,  1981)  and  speeds  up  to  16.8  ms1  (Petrie,  et  al .,  1996a).  In  general,  a  polymer  solution  is 
injected  into  a  developed  TBL  via  a  span-wise  slot  or  line  source  on  the  wall.  The  solution  concentration, 


C„  and  flux  per  unit  span,  Qh  are  varied  to  determine  their  influence  on  the  magnitude  and  extent  of  friction 
drag  reduction  that  is  achieved.  With  the  injection  of  high  concentration  polymer  solutions  (>  500  wppm ), 
the  drag-reduction  peaks  down-stream  of  injection  and  then  decreases  with  distance  from  the  injector.  The 
decrease  in  drag  reduction  with  downstream  distance  results  from  dilution  in  the  TBL  and  possibly  from 
flow-induced  polymer  degradation. 

Transport  and  dilution  of  the  injected  polymer  solution  in  a  TBL  have  been  measured  experimentally  by 
adding  a  dye  tracer  to  the  injected  polymer  and  extracting  samples  of  near-wall  fluid  downstream  from  the 
injection  location  (Fruman  and  Tulin,  1976;  Vdovin  and  Smol’yakov,  1978  and  1981)  or  by  using  laser- 
induced  fluorescence  of  a  dye  mixed  into  the  polymer  solution  (LIF)  (Fontaine  et  al. ,  1992).  Turbulent 
mixing  of  polymers  injected  into  the  near- wall  region  of  a  TBL  (Poreh  and  Hsu,  1972)  is  reduced  when 
compared  to  mixing  of  a  passive  scalar  (Poreh  and  Cermak,  1964),  and  this  is  consistent  with  the  observed 
reduction  in  skin  friction  in  the  polymer-injected  flow. 

Degradation  (chain-scission)  of  long-chain  polymer  molecules  in  the  turbulent  flow  is  also  of  practical 
concern.  Investigations  of  polymer  degradation  have  been  conducted  primarily  in  pipe  flows  and  rotational 
Taylor-Couette  flow  devices.  Currently  no  measurements  of  polymer  degradation  in  TBL  flows  have  been 
reported.  In  pipe  flows,  Horn  and  Merrill  (1984)  have  shown  that  the  polymer  chains  break  near  the  mid¬ 
point,  effectively  reducing  the  polymer  molecular  weight  by  50%.  Vanapalli  et  al.  (2005)  and  Vanapalli  et 
al.  (2006)  developed  relationships  relating  the  near-wall  shear  to  the  degradation  of  PEO  with  various 
molecular  weights,  and  their  results  indicate  that  chain-scission  occurs  in  PEO  solutions  at  shear  rates  much 
lower  than  those  of  previous  TBL  studies  and  the  present  experiments. 

Application  of  polymer  drag  reduction  to  large-scale  external  flows  necessitates  an  understanding  of  how 
the  injected  polymer  mixes  and  possibly  degrades  as  it  is  transported  downstream  of  the  injection  location. 
Moreover,  it  is  useful  to  understand  how  experimental  results  at  lower  Reynolds  numbers  compare  to 
similar  results  at  higher  speeds  and  larger  length  scales.  In  particular,  it  is  important  to  determine:  (1)  the 
downstream  persistence  of  TBL  drag  reduction  given  the  aforementioned  effects  of  mixing  and  polymer 


degradation;  (2)  the  local  near-wall  concentration  with  downstream  distance;  and  (3)  the  influence  of  the 
molecular  and  rheological  properties  of  the  injected  polymer  solutions  on  the  resulting  friction  drag 
reduction  and  near-wall  concentrations. 

Presented  here  are  the  results  of  large-scale  experiments  designed  to  examine  TBL  polymer  friction  drag 
reduction  on  a  smooth  surface.  A  comprehensive  data  set  was  gathered  during  the  experiments,  including 
the  local  average  skin  friction,  and  near-wall  velocity  and  concentration  fields  measured  with  micro- 
Particle  Image  Velocimetry  (PIV)  and  micro-Planar  Laser  Induced  Fluorescence  (PLIF)  systems.  Here,  the 
relationship  between  polymer  mean  molecular  weight,  injection  concentration,  injection  flux,  free-stream 
speed,  and  the  resulting  drag  reduction  and  turbulent  polymer  diffusion  is  reported.  The  relationship 
between  near- wall  polymer  concentration  and  the  drag  reduction  is  also  examined,  along  with  the 
possibility  of  TBL-flow-induced  polymer  degradation. 

2.  Experimental  Methods 

2.1.  Test  facility  and  model 

Experiments  were  conducted  in  the  William  B.  Morgan  Large  Cavitation  Channel  (LCC),  the  world’s 
largest  low-turbulence  (free  stream  turbulence  <  0.5%),  re-circulating  water  tunnel.  The  dimensions  of  the. 
test  section  are  16  m  (length)  and  3.05  m  x  3.05  m  (width  and  height).  Details  of  the  facility  and  its  use  can 
be  found  in  Etter  et  al.  (2005).  The  test  model  was  a  12.9  m  long,  3.05  m  wide,  and  0.18  m  thick  rigid  flat 
plate,  spanning  the  entire  width  of  the  LCC  test  section.  The  model  and  its  installation  in  the  LCC  is 
described  in  Sanders  et  al.  (2006),  so  only  a  brief  summary  is  provided  here.  The  model’s  leading  edge 
was  a  4:1  ellipse  and  its  trailing  edge  was  a  15°  full-angle  truncated  wedge.  The  model  was  positioned  near 
the  vertical  centerline  of  the  LCC  test  section  with  the  test  surface  facing  downward.  The  boundary  layer 
was  tripped  via  distributed  roughness  applied  to  the  first  0.30  m  of  the  model.  Otherwise,  the  test  surface 
was  polished  304  stainless  steel  with  a  root-mean-square  surface  roughness  of  k  <  0.4  nm,  small  enough  for 
the  surface  to  be  considered  hydro-dynamically  smooth.  The  test  surface  was  instrumented  with  a  suite  of 
flow  measurement  devices,  as  shown  in  Figure  1. 


2.2.  Polymer  mixing  and  injection 

The  polymers  used  in  this  experiment  were  polyethylene  oxide)  (PEO  or  POLYOX)  water-soluble  resins 
manufactured  by  Dow  Chemical  Company:  WSR-N-60K,  WSR-301,  and  WSR-308.  These  polymers  have 
manufacturer-supplied  mean  molecular  weights,  Mw ,  of  approximately  2, 4,  and  8  million,  respectively.  In 
the  case  of  WSR-N60K,  the  manufacturer’s  mean  molecular  weight  was  checked  by  means  of  gel 
permeation  chromatography  (GPC)  with  light  scattering  detection  (Vanapalli  et  a/.,  2005,  provides  details 
of  the  methods  used).  The  WSR-N60K  measured  Mw  and  polydispersity  MJMn  were  2.3  x  106  g  mole1 
and  2.44  respectively.  Attempts  to  assess  the  manufacturer-supplied  molar  mass  of  WSR-301,  and  WSR- 
308  by  the  same  technique  failed,  the  most  likely  reason  being  molecular  aggregation  (Vanapalli,  2006). 

The  raw  polymers  were  provided  as  dry  powders  that  must  be  hydrated  before  injection  into  the  TBL.  Such 
polymer  solutions  were  prepared  by  sifting  the  powder  into  a  jet  of  water  filling  a  reservoir.  An  activated 
carbon  filter  (Model  RT-2260-4,  Aquapure  Technologies)  was  used  to  remove  chlorine  from  the  water 
prior  to  mixing  because  of  chlorine’s  potential  for  chemical  degradation  of  PEO’s  drag  reduction  capability 
(Petrie  et  al. ,  2003).  After  the  initial  wetting,  the  polymer  solution  was  stirred  in  the  4.55  m3  reservoir 
(Chem  Systems,  Inc.)  with  a  4-blade  0.75  m  diameter  impeller  at  30  RPM  until  homogeneous;  this  process 
usually  took  1  to  3  days.  In  addition,  two  small  trolling  motors  were  used  near  the  reservoir  surface  to 
break-up  and  homogenize  any  remaining  agglomerations.  The  effectiveness  and  repeatability  of  the  mixing 
process  and  the  stability  of  the  polymer  solutions  were  verified  using  a  4.57  mm  diameter  pipe-flow 
apparatus.  The  pressure  drop  along  the  pipe  for  a  range  of  flow  rates  was  used  to  characterize  the  behavior 
of  the  polymer  solution.  Such  pressure-drop  test  results  spanning  115  hours  showed  that  the  polymer 
remained  stable  over  a  period  at  least  that  long,  and  this  is  consistent  with  the  reported  observations  of  Ellis 
etal  (1972). 

After  hydration  and  mixing,  the  polymer  was  transferred  to  a  smaller  1.140  m3  tank  which  fed  a  low  shear- 
rate,  stainless  steel  progressive  cavity  pump  (Moyno  model  2E012G1SSQAAA- 158587).  This  pump 
delivered  the  polymer  to  a  manifold  with  40  evenly  spaced  ports  along  the  injector  span.  Inside  the  injector 
manifold,  three  layers  of  perforated  plates  and  brass  screens  generated  a  pressure  drop  to  ensure  span-wise 


uniformity  of  the  injection  flow.  The  actual  injector  was  a  simple  tapered  slot  -  wider  at  its  inlet,  narrower 
at  its  outlet  -  oriented  with  mean  angle  of  25°  \vith  respect  to  the  working  surface  and  having  an  outlet  gap 
of  1  mm.  This  injector  is  similar  to  one  used  in  an  optimization  study  by  Walker  etal.  (1986). 

The  mixing,  pumping,  and  flow  upstream  of  the  injector  exit  of  the  polymer  solutions  may  all  degrade  the 
polymer  unintentionally  before  it  flows  into  the  TBL.  Vanapalli  et  al.  (2005)  conducted  a  series  of 
experiments  to  determine  the  critical  shear  rate  for  degradation  required  to  initiate  degradation  of  the  PEO 
polymer  solutions.  This  relationship,  presented  in  Section  4,  is  used  to  determine  the  critical  shear  rate  for 
degradation.  This  degradation  strain  rate  was  then  used  as  a  constraint  for  the  injection  manifold  design. 
In  particular,  care  was  taken  to  design  the  pumping,  piping,  and  distribution  systems  such  that  the 
maximum  shear  rates  were  at  least  a  factor  of  three  below  the  critical  value  for  degradation  initiation. 
Moreover,  samples  of  the  polymers  solutions  that  passed  through  the  pumping,  piping,  distribution,  and 
injection  systems  were  collected  at  the  injector  exit,  and  their  drag  reducing  potential  was  characterized 
from  pressure  drop  measurements  using  the  4.57-mra-diameter  pipe-flow  apparatus.  In  all  cases,  the  drag 
reducing  properties  of  the  delivered  polymer  were  found  to  be  unchanged  from  the  polymer  drawn  from  the 
upstream  reservoir. 

2.3.  Experimental  test  matrix 

A  total  of  four  independent  parameters  were  varied  throughout  the  polymer  tests.  There  were  three  free- 
stream  velocities,  6.65,  13.2,  and  19.9  ms1.  Three  different  PEO  polymers  having  different  values  of  Mw 
were  used,  and  for  each  polymer,  three  solution  concentrations  were  injected:  1000,  2000,  and  4000  wppm. 
Lastly,  polymer  solutions  were  injected  at  three  flow-rates  corresponding  to  Q{  =  0.14,  0.28,  and  0.71  ls~xm~ 
\  Following  established  convention,  the  polymer  solution  flow-rate  is  normalized  with  a  standard  volume 
flux  per  unit  span,  Qs=  67.3  v,  where  v  is  the  kinematic  viscosity  of  the  free  stream  liquid.  This  value, 
established  by  Wu  and  Tulin  (1972),  is  obtained  by  assuming  a  linear  velocity  profile  and  integrating  the 
streamwise  velocity  over  0  <  y+  <  1 1.6,  where  y+  is  the  wall  normal  distance  normalized  with  the  viscous 
length  ly  =  V'yjp/ Tw  ,  p  is  the  freestream  liquid’s  kinematic  viscosity,  and  rw  is  the  TBL’s  wall  shear  stress. 

The  three  polymer  flow  rates  in  this  study  represent  QJQS  =  2, 4,  10. 


2.4.  Test  procedure 

The  friction  drag  reduction  resulting  from  polymer  injection  was  measured  using  the  following  procedure. 
First,  a  volume  of  polymer  was  prepared  and  staged  for  injection.  The  free-stream  velocity  of  the  water 
tunnel  was  brought  to  zero,  and  the  skin-friction-drag  sensors  were  electronically  balanced  and  trimmed  to 
zero  output  at  the  desired  test  section  static  pressure.  Next,  flow  was  initiated  in  the  water  channel,  and 
when  the  desired  free-stream  velocity  and  static  pressure  were  achieved,  the  baseline  skin  friction,  r*,,  was 
measured.  At  this  point  polymer  injection  was  initiated,  and,  after  a  transient  period,  the  new  reduced 
friction,  was  measured.  When  polymer  injection  was  stopped,  the  baseline  friction  level  returned  after 
a  period  of  time  during  which  the  highly  viscoelastic  polymer  solution  continued  to  be  drawn  from  the 
injector  in  strands.  In  general,  rw  was  recovered  within  a  few  tens  of  seconds  or  a  few  minutes  after 
injection  of  polymer  solution  ceased. 

The  potential  buildup  of  active  polymer  in  the  tunnel  free-stream  flow  was  monitored  carefully.  If  the 
baseline  shear  stress  after  injection  was  reduced  by  more  than  2%,  a  controlled  amount  of  chlorine  was 
added  to  the  tunnel  water  to  degrade  the  background  polymer.  The  free  chlorine  is  consumed  as  it  reacts 
with  the  polymer,  so  the  amount  of  chlorine  added  between  runs  was  adjusted  to  leave  only  trace  amounts 
of  much  less  than  1  ppm.  The  background  chlorine  level  was  continually  monitored  to  insure  that  an  excess 
of  chlorine  did  not  occur.  As  a  second  measurement,  samples  of  the  tunnel  water  were  removed 
periodically  after  large  quantities  of  polymer  were  injected  and  tested  in  the  aforementioned  4.57  mm 
diameter  pipe  flow  apparatus  to  ensure  that  tunnel  water  mildly  contaminationed  with  degraded,  remnant 
polymer  had  the  same  pipe-flow  pressure-drop  characteristics  as  uncontaminated  water.  In  addition, 
throughout  the  course  of  these  multi-month  tests,  the  entire  LCC  (5,300  m3)  was  flushed  and  refilled  on 
weekends  when  polymer  had  been  used  during  the  prior  week. 

2.5.  Skin  Friction  Drag  Balances 

Direct  skin-friction  measurements  were  made  at  six  stream-wise  locations  using  floating-plate  strain-gauge 
type  friction  drag  sensors  at  the  locations  shown  on  Figure  1.  The  floating  plate  elements  of  the  sensors 
were  polished  17-4PH  stainless  steel,  15.2  cm  in  diameter,  0.79  cm  thick.  In  each  sensor,  the  floating  plate 


was  attached  to  the  load-sensing  element,  a  beryllium  copper  flexure  instrumented  with  a  full  Wheatstone 
bridge  of  semiconductor  strain  gages.  The  plate  and  flexure  were  mounted  within  a  bowl-shaped  housing 
so  that  the  plate  and  housing  edge  were  flush  with  the  model’s  test  surface.  The  gap  between  the  floating 
plate  and  housing  was  measured  with  feeler  gauges  to  be  60  ±  20  pirn,  and  varied  slightly  from  sensor  to 
sensor.  The  strain  gages  were  excited  at  7  VDC  using  signal-conditioning  amplifiers  (Model  #2310, 
Vishay  Measurements  Group).  The  sensor  outputs  were  amplified  500  fold  and  low-pass  filtered  at  100  Hz 
with  the  same  Vishay  units.  The  conditioned  signals  were  recorded  simultaneously  at  50  Hz  with  a 
National  Instruments  NI-DAQ  data  acquisition  board  and  a  LabView  virtual  instrument.  The  typical  data 
collection  intervals  were  approximately  one  to  three  minutes. 

The  skin-friction  sensors  were  calibrated  using  a  four-point  calibration  from  zero  to  two  pounds  with  both 
increasing  and  decreasing  loads  and  a  second  load  cell  (Model  LCEB-5,  Omega  Engineering)  mounted 
horizontally  downstream  of  the  friction  balance.  This  second  load  cell  was  calibrated  from  zero  to  two 
pounds  using  weights.  A  thin  cable  tied  between  the  second  load  cell  and  a  suction  cup  on  the  floating 
plate  element  of  the  sensors  was  used  to  apply  the  load.  The  second  load  cell  was  mounted  on  a  linear 
traverse  used  to  increase  or  decrease  the  load  on  the  sensor.  This  procedure  eliminated  friction  problems 
associated  with  pulleys.  The  error  associated  with  the  sensor  slope  was  taken  to  be  twice  the  standard . 
deviation  of  multiple  calibration  slopes  (typically  five  or  more)  and  varied  between  0.8%  and  2.9%.  The 
zero-level  of  the  load  cells  in  the  friction  balances  drifted  with  changes  in  free-stream  temperature  and 
pressure,  so  corrections  were  built  into  the  calibration  and  data  reduction  process.  The  overall  uncertainty 
in  any  particular  time-averaged  measurement  of  wall  shear  stress  from  one  sensor  was  typically  ±  5%. 

For  the  present  study,  skin  friction  was  measured  to  determine  the  drag  reduction,  %DR ,  obtained  in 
polymer-laden  flow: 


%DR  = 


(2.1) 


Here,  sensor  drift  errors  are  largely  eliminated  because  ^  and  the  steady  shear  stress  measured  during 
polymer  injection,  twJ,  were  measured  within  a  few  minutes  of  each  other.  Thus,  the  uncertainties  for  the 
reported  %DR  values  are  also  ±  5%  of  the  reported  value. 

2.6.  Polymer  concentration  measurements 

Planar  Laser  Induced  Flourescence  (PLIF)  was  used  to  measure  the  polymer  concentration  in  the  near-wall 
region  of  the  TBL.  Measurements  were  made  at  three  stream-wise  locations  X-X,  =  0.64,  4.62,  and  9.36  m, 
where  X  is  the  stream-wise  distance  downstream  of  the  leading  edge  of  the  test  model,  and  X,is  the  stream- 
wise  location  of  the  injector  slot.  A  known  concentration  of  fluorescent  dye  (Rhodamine  6G  supplied  by 
Sigma  Chemical)  was  thoroughly  mixed  into  the  injected  polymer  solutions.  The  dye-laden  polymer 
solution  was  illuminated  with  a  light  sheet  normal  to  the  plate  surface  and  aligned  with  the  mean  flow 
direction.  Light  sheets  were  formed  from  the  beams  of  pulsed  Nd-YAG  lasers  operating  at  the  532  nm 
wavelength  and  were  measured  to  be  ~125  fim  thick  at  the  plate  surface.  The  plane  of  illuminated  polymer 
and  dye  was  imaged  with  a  high-resolution  camera  (LaVision  Imager  Pro),  where  the  field  of  view  was 
approximately  5  mm  (stream-wise)  by  2  mm  (wall-normal).  Periscope  prisms  protruding  5  mm  from  the 
test  surface  were  used  to  redirect  the  line  of  sight  of  the  cameras  housed  in  waterproof  boxes  inside  the  test 
model.  The  prism  was  located  4  mm  from  the  light  sheet.  A  long-pass  optical  filter  was  used  to  attenuate 
the  original  532  nm  light  from  the  laser  and  return  only  the  Stokes-shifted  light  from  the  fluorescent  dye. 
Finally,  a  pair  of  cylindrical  lenses  was  used  to  stretch  the  wall-normal  component  of  the  images  by 
approximately  a  factor  of  2.5.  Further  details  are  provided  in  Winkel  etal.  (2006). 

The  intensity  of  the  light  collected  by  the  cameras  was  used  to  determine  the  polymer  concentration  via 
separate  calibration  measurements  completed  as  follows.  A  calibration  cell  was  filled  with  -  2500  cm 3  of 
dyed  polymer  solution  and  fitted  onto  the  plate  surface,  immersing  the  light  sheet  and  collection  optics.  A 
pump  was  used  to  slowly  circulate  the  solution  to  prevent  local  heating  of  the  liquid  due  to  repeated  laser 
pulses.  The  fluorescent  return  from  these  solutions  was  used  to  relate  the  intensity  at  each  imaging  pixel 
with  the  dye  concentration,  accounting  for  variations  in  light  sheet  intensity.  As  discussed  below,  the 
dilution  level  of  the  injected  polymer  solutions  varied  between  factors  of  10  and  1000,  depending  upon  the 


experimental  conditions  and  the  measurement  location.  Hence,  care  was  taken  to  set  the  initial  dye 
concentration  and  calibration  range  such  that  the  PLIF  measurement  at  a  particular  location  produced 
sufficient  signal-to-noise  ratio.  However,  it  was  usually  not  possible  to  select  a  single  dye  concentration 
that  would  result  in  acceptable  measurements  at  all  three  PLIF  measurement  location,  necessitating 
repeated  experimental  runs  at  multiple  injected  dye  concentrations.  The  uncertainty  of  the  PLIF 
measurements  varied  with  the  amount  of  dilution,  but  is  typically  ±20%  for  dilution  of  less  than  a  factor  of 
100  of  the  initial  concentration,  and  ±30%  for  dilutions  greater  than  a  factor  of  100.  Background  light  and 
electronic  noise  were  the  limiting  factors  for  accurate  large-dilution  concentration  measurements. 

3.  Baseline  Flow  Results 

The  baseline  flow  measurements  from  this  experiment  are  summarized  in  Table  1  for  the  six  skin-friction 
measurement  locations  ( X  -  X,  =  0.64,  2.09,  4.62,  6.11,  7.91,  9.36  m),  and  the  three  PLIF  measurement 
stations.  Boundary  layer  growth  on  the  model  and  test  section  side-walls  caused  a  slight  favorable  pressure 
gradient  along  the  plate;  the  three  free  stream  test  speeds  -  6.65,  13.2,  and  19.9  ms~]  -  are  valid  for  the 
nominal  streamwise  center  of  the  plate.  Full  baseline  boundary  layer  profiles  were  measured  at  the  same 
locations  as  the  third  and  sixth  shear  stress  sensors  with  an  LDV  system.  The  local  free  stream  velocity  U , 
the  99%  boundary  layer  thickness  (5^,  and  the  momentum  thickness  <9,  inferred  from  these  measurements, 
are  listed  in  bold  in  Table  1.  The  remaining  entries  for  U  in  Table  1  were  inferred  from  static  pressure 
measurements  made  along  the  plate  and  the  steady  Bernoulli  equation.  The  remaining  entries  for  5*?  and  0 
in  Table  1  were  obtained  from  simple  power-law  fits  to  the  LDV-profile-determined  values  of  <5^,  and  0. 
In  all  cases,  the  values  obtained  from  the  fits  were  within  a  few  percent  of  what  would  be  obtained  from 
standard  flat-plate  TBL  formulae.  In  addition,  Table  1  provides  values  for  the  downstream-distance-based 
Reynolds  number,  Rejr  =  f/x/v,  the  skin  friction  coefficient,  Cf=Twl\pU2,  the  acceleration 

parameter  KH=  (v/  U2)(dU  /dx),  the  friction  velocity  ux  =  (Tw  /  p)1/2  ,  the  viscous  length  and  the  ratio, 

of  the  test-surface  roughness  height  divided  by  ly.  When  K'  <  10"*,  the  flow  acceleration  should  be 
sufficiently  mild  so  that  there  will  be  little  or  no  deviation  in  the  boundary  layer’s  profile  from  the  typical 


log-law  (Patel  1965);  for  the  current  experiments,  K '  <  10'9.  And,  when  k+  ^  1  or  so,  the  plate  surface  can 
be  considered  hydro-dynamical ly  smooth. 

Baseline  skin-friction  coefficients  measured  in  this  experiment  are  shown  in  Figure  2  with  the  standard 
friction  correlations  of  White  (1991),  Eq.  (3.1),  and  Schultz-Grunow  (1941),  Eq.  (3.2), 

Cf  =  0.455  In'2  (0.06  Re ,)  (3.1) 

C7  =  0.370  log’2  584  (Re,)  (3.2) 

Skin-friction  results  from  two  separate  tests  conducted  in  2004  and  2005,  termed  Phases  III  and  IV,  are 
provided  to  demonstrate  the  repeatability  of  the  measurements.  The  measured  skin-friction  is 
approximately  5%  below  that  predicted  by  White’s  correlation,  but  shows  good  agreement  with  the 
Schultz-Grunow  correlation.  A  power-law  curve  fit  (solid  line)  to  the  experimental  data  is  also  provided, 
given  by 


Cf  =  (0.0170  ±0.0004)  •  Re^01237±00067)  (3.3) 

The  Reynolds  numbers  given  here  employ  the  kinematic  viscosity,  v  =  9.95 xlO-7  m2s  \  based  on  the 
average  temperature  of  the  water  (20.4  °C).  A  complete  presentation  of  the  baseline  flow  measurements 
obtained  during  these  experiments,  including  turbulence  quantities  and  near-wall  velocity  measurements,  is 
provided  in  Oweis  et  al.  (2007). 

4.  Polymer  Rheology 

The  rheology  of  high  molar  mass  aqueous  PEO  solutions  is  complex.  The  relationship  between  the 
mechanical  behavior  of  PEO  solutions  and  the  detailed  characteristics  and  interactions  of  individual  PEO 


molecules  continues  to  be  an  active  area  of  investigation  (Vlassopoulous  and  Schovvalter,  1994; 


Kalashnikov,  1994;  Tirtaatmadja  et  al. ,  2006).  For  drag  reduction,  the  most  important  characteristic  of  the 
high  molecular  weight  polymer  is  its  intermittent  extension  and  relaxation  in  the  turbulent  flow. 
Treatments  of  polymer-laden  turbulent  flow  have  captured  this  physical  process  by  means  of  various 
constitutive  models  of  flow-polymer  interactions.  For  example,  the  FENE-P  (Finite  Elastic  Non-Linear 
Extensible-Peterlin)  dumbbell  constitutive  model  has  been  widely  used  because  it  incorporates  the 
rheological  effects  of  polymer  finite  extensibility  in  a  computationally  efficient  way  (Beris  and 
Dimitropoulos,  1999;  Housiadas  and  Beris,  2003;  Dubief  et  al.,  2004;  and  Almeida  et  al.,  2006).  The 
FENE-P  dumbbell  model  relates  the  polymer  stress  tensor  to  the  local,  pre-averaged  extension  of  the  dilute 
polymer  molecules.  The  distribution  of  length  and  orientation  of  the  polymer  molecules  determines  the 
polymer  conformation  tensor.  The  polymer  stresses  in  the  flow  are  then  a  function  of  the  polymer 
conformation  tensor.  Molecular  parameters  in  the  FENE-P  dumbbell  model  are  the  equilibrium  polymer 
coil  size,  the  maximum  extensibility  of  the  polymer  molecule  and  the  polymer  concentration.  These 
experimental  parameters  implicitly  determine  rheological  quantities  such  as  the  zero-shear  viscosity,  the 
extensional  viscosity  and  the  polymer  relaxation  time  (Dealy  and  Larson,  2006).  Formally,  the  FENE 
dumbbell  models  are  appropriate  as  a  molecular  description  of  dilute  polymer  solutions  below  the  overlap 
concentration,  c*. 

The  solutions  of  high  molecular  PEO  injected  into  the  turbulent  boundary  layer  in  the  present  study  were  of 
molecular  weight,  Mw ,  2.3,  4,  and  8  million  g  mole1 ,  and  of  concentrations  1000,  2000,  or  4000  wppm.  The 
intrinsic  viscosity,  [rj\0,  of  PEO  solutions,  in  units  of  cm3  g1,  can  be  determined  from  Mw  using  the  Mark- 
Houwink  relationship  provided  by  Bailey  and  Callard  (1959): 

[tj]0  =  0.0125  AO78  (4.1) 

The  dilute  limit  for  polymer  solutions  is  then  determined  by  a  common  definition  of  the  overlap 
concentration: 

C*  =  l/[71]0  (4.2) 

Table  2  presents  the  results  for  the  three  PEO  molecular  weights  used  here.  Note  that  the  injected 
concentrations  exceed  c*  for  all  the  conditions  of  this  study,  and  there  is  a  substantial  portion  of  the  flow 


where  the  concentration  exceeds  c*.  As  the  polymer  mixes  with  the  free-stream  fluid,  the  concentration 
will  be  reduced,  passing  below  c*  at  some  distance  downstream  of  the  injector. 


The  shear  viscosity  of  the  polymer  solution  is  a  function  of  polymer  molecular  weight,  concentration,  and 
shear  rate.  Kalashnikov  (1998)  provides  a  relationship  for  A  =  (7]0  —TJ^)  /  ,  where  7]0  and  7]^ are  the 

limiting  shear  viscosities  at  zero  and  infinite  shear  rates: 


(4.3) 


where  C  is  the  concentration  (mg  mg1)  and  [rj)0  is  the  intrinsic  viscosity  (cm3g '),  and  the  constants  in  (4.3) 


are  dimensional.  Table  2  lists  the  value  AC  for  injected  solutions.  In  the  low  shear-rate  regime,  the  shear 


viscosity  of  the  highly  concentrated  injected  solutions  will  be  much  higher  than  the  solvent  viscosity,  but 


the  solutions  will  undergo  substantial  shear  thinning. 

Kalashnikov  (1998)  also  examined  the  polymer  relaxation  time  for  PEO  solutions,  and  developed  a 
relationship  for  0*for  PEO  solutions  with  [tj]0>  280  cm3g! : 


(4.4) 


where  Tc  is  the  mixture  temperature  in  °C.  Kalashnikov  (1998)  reports  that  this  relationship  is  generally 
valid  over  a  range  of  concentrations  of  order  5  to  2000  wppm  for  a  PEO  solution  with  [t]]0  =  2650  cm 3g  I. 
Note,  however,  that  the  relaxation  time  goes  to  zero  at  zero  concentration,  and  an  alternate  expression  for 
the  relaxation  time,  such  as  the  Zimm  time,  0Z,  should  be  used  for  very  low  concentrations.  The  Zimm 
time,  appropriate  for  single  molecules  immersed  in  a  good  solvent,  is: 


6Z  =  0.422 


[^7  Ip  Mwvjs 
RT 


(4.5) 


where  R  is  the  ideal  gas  constant  (Dealy  and  Larson,  2006).  Shin  and  Shaqfeh  (2005)  report  relaxation 
times  that  are  somewhat  higher  that  those  produced  by  the  Kalashnikov  relationship  for  WSR-301,  but  both 


groups  report  that  the  relaxation  times  of  the  concentrated  polymer  solutions  are  much  larger  that  the  Zimm 
time. 

In  a  strong  flow,  flexible  PEO  molecules  can  experience  a  sufficient  deformation  rate  to  be  oriented  and 
extended.  In  turbulent  flows,  this  condition  is  the  onset  deformation  rate  for  friction  drag  reduction  (Virk  et 
al.,  1967).  Although  the  detailed  validity  of  the  molecular  interpretation  as  a  coil-stretch  transition  is 
unclear  above  the  overlap  concentration,  the  onset  condition  is  a  useful  quantifier  in  both  conduit  and 
boundary  layer  flows.  This  quantity  has  been  characterized  for  varying  aqueous  concentrations  of  PEO;  the 
data  provided  by  Virk  (1975)  are  used  here  to  determine  a  dimensional  relationship  for  the  minimum  shear 
rate  for  drag  reduction  onset,  y  *  (Vanapalli  et  al. ,  2005): 

y*  =  (3.35*109)AV'  (4.6) 

Another  possible  determinant  of  drag  reduction  behavior  is  polymer  scission  and  degradation.  Vanapalli  et 
al.  (2005)  describe  a  series  of  experiments  to  determine  the  critical  shear  rate  for  degradation,  yD ,  of  PEO 
solutions. 

yD  =  3.23  *  1018  *  M~220  (4.7) 

where  Mw  is  the  molecular  weight  in  g  mole1,  and  yD  is  in  s'1.  Table  2  presents  these  rheological 
parameters  for  the  three  polymer  types  used  in  this  study. 

A  potential  complication  of  assessing  the  polymer  rheology  is  that  the  treatments  discussed  above  are 
applicable  to  polymer  chains  described  by  a  well-characterized  molar  mass.  However,  as  discussed  in  the 
literature,  high  molar  mass  PEO  is  susceptible  to  intermolecular  aggregation  in  aqueous  solution  (Dunlop, 
and  Cox,  1977;  Polverari  and  van  de  Ven,  1996;  Faraone  et  al.,  1999;  and  Ho  et  al.,  2003).  For  the 
solutions  considered  here,  the  polymer  molar  mass  is  ill  defined  in  the  context  of  the  above  relationships. 
Moreover,  the  molar  mass  is  difficult  to  characterize  and  is  potentially  a  sensitive  function  of  sample 


preparation  and  deformation  history.  If  aggregation  occurs  as  the  polymer  molar  mass  is  increased  above  a 
certain  value,  then  equations  based  on  scaling  (such  as  4.6  and  4.7)  will  perform  poorly.  Vanapalli  (2006) 
discusses  light  scattering  observations  that  are  suggestive  of  aggregation  behavior  in  WSR-301.  This 
observation  is  likewise  consistent  with  the  finding  that  WSR-301  relaxation  times  greatly  exceed  Zimm 
model  predictions,  even  below  the  overlap  concentration  (Shin  and  Shaqfeh,  2005;  Tirtaatmadja  et  al„ 
2006).  Keeping  this  important  factor  in  mind,  these  relationships  for  the  rheological  properties  of  PEO  will 
be  used  to  help  interpret  results. 

5.  Drag  Reduction  Resulting  From  Polymer  Injection 

5.7.  Stages  of  polymer  mixing 

Injection  of  polymer  solution  into  the  TBL  results  in  a  distribution  of  polymer  across  the  boundary  layer 
that  evolves  in  the  stream- wise  direction.  Poreh  and  Hsu  (1972)  describes  four  mixing  stages:  initial, 
intermediate,  transition,  and  final.  In  the  initial  stage,  the  injected  polymer  is  largely  near  the  wall  and 
within  the  sub-layer.  In  the  intermediate  stage,  the  polymer  lies  within  a  layer  thinner  than  the  boundary 
layer  thickness  but  thicker  that  the  viscous  sub-layer  and  buffer  layer.  In  the  transition  stage,  the  polymer 
is  diffused  throughout  most  of  the  TBL.  In  the  final  stage,  the  growth  rate  and  thickness  of  the  polymer¬ 
laden  portion  of  the  TBL  is  the  same  as  that  of  the  TBL  itself.  This  process  can  be  compared  to  the 
similarly  defined  mixing  stages  of  a  passive  scalar,  discussed  by  Poreh  and  Cermack  (1964).  Since  the 
injected  polymer  flux  is  usually  larger  than  Qs,  the  layer  of  injected  polymer  is  almost  always  thicker  than 
the  viscous  sub-layer,  and  hence  it  is  artificial  to  distinguish  between  the  initial  and  intermediate  stages  as 
defined  by  Poreh  and  Cermack  (1964).  Therefore,  these  stages  are  combined  herein  and  referred  to  this  as 
the  development  stage.  Although  these  stages  are  defined  for  the  mixing  process,  the  skin-friction  drag 
reduction  has  a  corresponding  character  in  each  of  these  regimes. 


5.2.  Drag  reduction  in  the  development  stage 

The  concentration  needed  to  achieve  MDR  varies  with  the  molecular  weight  and  other  physical  properties 
of  the  polymer  and  the  topology  and  Reynolds  number  of  the  flow,  but  typical  values  for  PEO  at  moderate 


to  high  Reynolds  numbers  range  from  10  to  100  wppm.  Consequently,  injection  of  PEO  solutions  with 
concentrations  on  the  order  of  1000  wppm  in  the  near  wall  region  should  result  in  MDR.  Moreover,  MDR 
should  persist  for  a  significant  downstream  extent  after  injection,  possibly  even  into  the  transition  stage  of 
mixing,  creating  an  extended  zone  or  plateau  of  MDR  some  distance  downstream  of  the  injector.  However, 
such  an  MDR  plateau  was  not  observed  in  these  experiments.  Instead,  a  region  of  increasing  drag 
reduction  immediately  downstream  of  the  injection  location  was  observed,  a  phenomenon  also  noted  in 
Petrie  etal.  (1996)  and  Petrie  etal.  (2005). 

The  length  of  the  development  region  where  the  drag  reduction  is  growing  was  observed  to  increase  with 
increasing  polymer  ^(concentration),  and  Q,  (injection  rate),  and  with  decreasing  free-stream  speed. 
In  the  present  study,  this  development  region  was  observed  at  distances  exceeding  5  m  for  the  lowest  free- 
stream  speed  (6.65  ms1)  at  the  highest  Mw  and  injection  rates.  The  development  region  was  never  observed 
at  19.9  ms 1  and  only  in  two  cases  at  13.2  ms 1 .  Two  examples  of  the  development  region  for  the  WSR-301 
polymer  are  shown  in  Figure  3  which  shows  %DR  defined  by  (2.1)  as  a  percentage  vs.  distance  downstream 
from  the  injector,  X  -  X(  (m). 

Polymer  concentration-field  measurements  show  that  the  polymer  is  highly  non-uniform  in  the , 
development  region  and  in  many  cases  the  polymer  is  observed  in  long  filaments.  In  some  cases,  the 
polymer  concentration  field  is  nearly  binary,  essentially  either  the  injected  concentration  C,  or  nearly  zero 
polymer  concentration.  Polymer  solutions  at  high  concentration  have  a  very  large  extensional  viscosity, 
and  the  polymer  does  not  exit  the  injection  slot  in  a  two-dimensional  sheet.  Instead,  the  polymer  is  drawn 
from  the  injector  by  the  high-speed  water  flow.  This  results  in  a  region  of  highly  intermittent  polymer- 
water  mixture.  Also,  the  shear  viscosity  of  the  polymer  solutions  can  be  several  times  that  of  the  pure 
water,  and  this  can  locally  increase  the  skin  friction  near  the  injector,  even  though  there  may  be  a  reduction 
in  turbulent  momentum  transport.  Lastly,  the  injected  polymer  solutions  have  relaxation  times  on  the  order 
of  0.1  seconds  at  concentrations  near  C,.  A  finite  extent  of  stream-wise  flow  is  needed  for  the  polymers  to 
stretch,  mix,  and  interact  with  the  TBL  turbulence.  Moreover,  once  it  has  passed  the  injector,  the  boundary 
layer  turbulence  requires  a  finite  time  (or  distance)  to  adjust  to  the  new  polymer-modified  regenerative 


cycle  that  leads  to  high  levels  of  drag  reduction.  Consequently,  the  maximum  drag  reducing  potential  of 
the  polymer  solution  is  reached  some  distance,  perhaps  several  meters,  downstream  of  the  injector. 


5.3.  Maximum  drag  reduction 

The  maximum  levels  of  drag  reduction  observed  in  the  current  tests  are  70,  67,  and  60  %DR  at  6.65,  13.2, 
and  19.9  ms  ',  respectively.  These  levels  of  drag  reduction  are  observed  with  the  highest  flux  (C,  =  4000 
wppm,  Qt  =  10 Q,)  of  the  highest  Mw  polymer  (8  million),  WSR-308,  and  are  shown  in  Figure  4.  The  peak 
drag  reduction  was  observed  at  the  first  measurement  station  (X  -  X,  =  0.64  m)  only  at  19.9  ms1.  As 
discussed  in  the  previous  section,  peak  drag  reductions  occurred  approximately  6  and  2  m  downstream  of 
the  injector  at  6.65  and  13.2  ms  ',  respectively.  The  maximum  level  of  drag  reduction  observed  in  these 
experiments  agrees  with  prior  work.  Vdovin  and  Smol’yakov  (1981)  observed  maximum  drag  reduction 
between  65  and  70  percent  for  injection  of  WSR-301  solutions.  Petrie  et  al.  (1996)  report  maximum  drag 
reduction  between  64  and  67  percent  for  a  maximum  injection  concentration  of  2000  wppm  WSR-301.  As 
Mw  increases,  the  amount  of  drag  reduction  that  is  achieved  for  a  given  value  of  Q,C,  also  increases.  WSR- 
308  has  an  increase  of  around  1.1  times  %DR  compared  to  WSR-301,  and  both  are  higher  than  WSR-N60K 
by  10 %DR  to  20 %DR.  These  differences  are  less  than  what  is  expected  based  solely  on  molecular  weight 
or  intrinsic  viscosity,  and  this  will  be  discussed  further  below. 

5.4.  Drag  reduction  in  the  transitional  and  final  mixing  regions 
The  transition  and  final  mixing  regions  occur  once  the  level  of  %DR  peaks  and  begins  to  fall  with 
increasing  downstream  distance  (Poreh  and  Hsu  1972).  Figures  4  and  5  illustrate  how  drag  reduction 
decays  with  downstream  distance  from  the  injector.  Warholic  et  al.  (1999)  suggested  that  aggregates  play 
an  important  role  in  polymeric  drag  reduction.  In  their  experiments  they  injected  equivalent  fluxes  of 
polymer  while  varying  the  injected  polymer  concentration  and  injection  rates.  For  equivalent  fluxes  of 
polymer,  drastically  larger  drag  reductions  were  observed  with  500  wppm  solution  than  100  wppm,  45%  vs. 
19%,  respectively.  This  indicates  that  the  higher  concentration  polymer,  where  aggregates  are  more  likely 


to  occur,  will  be  more  effective  at  reducing  drag.  In  contrast,  results  from  the  work  of  Vdovin  and 
Smol’yakov  (1978,  1981)  indicate  that  the  product  of  injection  rate  and  concentration,  Q,Ch  will  scale  the 
drag  reducing  capacity  of  the  injected  solution.  Results  from  the  present  experiment  are  shown  in  Figure  5 
for  two  pairs  of  equivalent  fluxes  of  WSR-301  at  two  concentrations.  The  present  results  indicate  that  the 
total  polymer  flux,  Q,Cit  is  the  primary  factor.  Since  the  experiments  of  Warholic  et  al.  (1999)  were 
performed  in  a  relatively  slow  channel  flow  with  speeds  up  to  2  ms1,  the  role  of  aggregation  may  have  been 
more  pronounced. 


5.5.  K- Factor  scaling  of  the  drag  reduction 

A  simple  scaling  relationship  for  the  observed  drag-reduction  based  on  the  polymer  injection,  downstream 
distance,  and  flow-speed  was  developed  by  Vdovin  and  Smol’yakov  (1981).  They  observed  that  the  near- 
wall  concentration  decreased  exponentially  in  the  transition  region,  and  proposed  a  scaling  parameter  that 
led  to  a  reasonable  collapse  of  their  drag-reduction  data.  Their  parameter,  K,  is  given  by 


K  = 


(5.1) 


where  X,  is  the  distance  from  the  line  source  ejecting  the  polymer.  When  the  percentage  of  drag  reduction, 
%DR,  is  plotted  against  K,  data  from  the  intermediate  region  show  a  logarithmic  decrease  in  drag  reduction 
with  decreasing  K  (i.e.  increasing  X,  after  fixing  the  other  variables).  Poreh  and  Hsu  (1972)  demonstrate 
that  such  behavior  can  be  explained  using  a  Lagrangian  similarity  hypothesis  for  mixing  of  the  polymer  in 
the  logarithmic  layer  of  the  TBL,  as  first  proposed  by  Batchelor  (1957).  The  rate  of  mixing  is  related  to 
the  local  friction  velocity  that,  in  turn,  is  a  function  of  the  near-wall  polymer  concentration.  Consequently, 
the  mixing  rate  is  not  expected  to  be  a  simple  logarithmic  relationship.  Nevertheless,  the  K  scaling  has 
proven  useful  as  a  way  of  organizing  and  comparing  %DR  data,  and  it  can  help  reveal  some  underlying 
flow  processes. 


Figure  6  presents  %DR  vs.  K  for  the  three  polymers  used.  The  development  region  is  evident  for  higher 
values  of  K  ( e.g .  the  measurements  closest  to  the  injector).  %DR  grows  and  reaches  a  maximum  near  K  ~ 
2x1 0~7,  then  decreases  with  decreasing  K  value  ( e.g .  increasing  distance  from  the  injector).  The  product 


QjCj  adequately  collapses  the  data  for  a  given  speed,  suggesting  that  it  is  the  total  flux  of  polymer  that  is  the 
relevant  quantity,  whether  it  is  slowly  injected  at  high  concentration  or  rapidly  injected  at  low 
concentration,  over  the  range  of  variables  of  these  tests.  However,  the  length  of  the  initial  region  varies 
with  free-stream  speed,  with  the  location  of  the  peak  drag  reduction  decreasing  with  increasing  U.  In  turn, 
increasing  Mw  extends  the  development  region  and  level  of  maximum  drag  reduction,  with  the  differences 
between  the  low  molecular  weight  polymer,  WSR-N60K  much  more  pronounced  compared  to  the  two 
higher  molecular  weight  polymers,  WSR-301  and  WSR-308.  Figure  7  illustrates  the  speed  dependence  of 
the  K  plots  for  WSR-301.  The  slopes  of  %DR  vs.  K  for  the  transition  region  are  similar,  but  the  start  of  the 
transitional  region  varies  with  speed.  Also  plotted  is  a  line  representing  the  best-fit  slope  for  the 
transitional  regime  reported  by  Vdovin  and  Smol’yakov  (1981)  for  WSR-301. 

The  slopes  of  the  present  data  are  similar  to  this,  although  the  length  of  the  development  region  varies  with 
speed.  Table  3  describes  the  parameters  and  conditions  under  which  the  data  were  collected.  Vdovin  and 
Smol’yakov  (1981)  only  report  an  overall  range  of  speeds  and  do  not  distinguish  their  %DR  data  with  their 
specific  injection  conditions  and  free-stream  speed.  Moreover,  their  stated  aim  was  to  examine  the 
transitional  region,  and  they  avoided  injection  conditions  that  would  lead  to  measurement  of  an  initial  and 
intermediate  region,  and  do  not  report  any  data  illustrating  this  regime.  But,  their  data  are  similar  in  slope 
and  maximum  %DR  to  the  low  and  medium  speed  data  from  the  present  study,  along  with  the  approximate 
K  value  for  the  beginning  of  the  transition  region. 

Examination  of  other  published  results  for  WSR-301  also  shows  considerable  variability.  Figure  8  presents 
%DR  vs.  K  for  several  experiments,  sorted  into  low,  medium,  and  high-speed  investigations.  Table  3 
presents  details  on  the  experiments  and  how  these  data  were  collected.  There  is  considerable  spread  in  the 
plot,  and  the  maximum  drag  reduction  level  varies  considerably.  Recall  that  in  Fontaine  et  al.  (1992)  and 
Petrie  et  al.  (2003),  a  relatively  large  drag  balance  was  employed  that  could  mask  the  initial  region  and 
peak  level  of  drag  reduction.  Moreover,  in  Fontaine  et  al.  (1992)  and  Petrie  et  al.  (1996),  the  slot 
comprised  only  the  middle  50%  of  the  span.  Span-wise  diffusion  of  the  polymer  is  likely  limited  near  the 
injector  and  at  higher  concentration.  However,  as  the  polymer  dilutes,  spanwise  mixing  may  take  place, 


and  this  may  be  responsible  for  the  steeper  decrease  in  %DR  with  decreasing  K  for  these  data  sets.  Such 
change  in  slope  between  30  %DR  and  40  %DR  can  be  discerned  in  the  data  from  these  experiments, 
although  the  scatter  in  the  data  makes  definitive  conclusions  difficult  The  low  speed  data  are  most 
scattered,  and  the  location  of  peak  drag  reduction  varies  most.  As  the  free-stream  speed  increases,  K  at 
maximum  %DR  increases,  e.g.  maximum  %DR  moves  upstream.  In  the  present  experiments,  the  peak  in 
%DR  was  not  observed  at  the  highest  speed,  since  it  occurred  before  the  first  shear-stress  measurement 
station.  The  data  of  Petrie  et  al.  (1996)  and  Petrie  et  al.  (2005)  appear  to  show  a  maximum  in  %DR  near 
the  injector,  although  it  may  also  be  a  plateau.  This  variability  in  the  %DR  vs.  K  plots  is  not  surprising, 
considering  that  the  K  parameter  is  essentially  a  similarity  variable  based  on  outer  scales  that  does  not 
incorporate  the  detailed  physical  processes  involved  in  the  different  mixing  zones.  It  is  interesting  to  note 
that  the  logarithmic  rate  of  decrease  in  %DR  is  similar  for  these  cases,  at  least  in  the  region  of  drag 
reduction  greater  than  -30  %DR ,  when  the  flow  is  in  this  transitional  mixing  region,  and  polymer  is 
moving  into  the  log-layer. 

One  should  also  consider  the  physical  dimension  of  the  initial  mixing  region  in  the  present  and  previous 
experiments.  In  the  present  case,  the  peak  of  the  drag  reduction  usually  lays  between  the  injector  and  the 
second  measurement  station,  which  was  2.09  meters  from  the  injector.  Table  3  shows  the  range  of  the. 
development  length,  XD  defined  as  the  distance  between  the  injection  and  the  location  where  K  =  2x1  O'7,  for 
the  range  of  Q,C,  in  each  experiment  at  U  =  -7  and  -20  ms1.  For  the  smaller  scale  experiments,  it  is  often 
not  possible  to  resolve  the  longest  initial  regions  at  the  lower  speeds  since  XD  can  extend  beyond  the 
farthest  measurement  location.  In  the  present  experiment,  XD  is  the  largest  at  the  lowest  speed  and  highest 
QjCh  and  this  distance  is  comparable  to  the  experiments  of  Vdovin  and  Smofyakov  (1981)  and  about  a 
meter  longer  than  that  reported  by  Petrie  et  al.  (1996).  At  the  highest  free-stream  speeds,  XD  is  reduced, 
and  there  is  less  spread  in  the  %DR  vs.  K  data  from  the  different  experiments.  The  %DR  vs.  K  plots  also 
reveal  the  end  of  the  transition  and  final  mixing  zones,  where  the  polymer  has  filled  the  boundary  layer  and 
mixes  at  a  rate  similar  to  that  of  a  passive  scalar.  Here,  the  %DR  has  decreased  to  less  than  10%.  The  K 
value  at  which  transition  to  the  final  zone  mixing  occurs  varies  with  Af*,,  but  is  generally  between  2  x  109  < 


K<6x  10 9. 


Variation  in  the  geometry  of  the  injection  slots  used  in  these  studies  as  well  as  variation  in  the  exit  velocity 
of  the  injected  polymer  (typically  1%  to  10%  of  the  free-stream  speed)  did  not  appear  to  strongly  influence 
the  length  of  the  initial  region.  Vdovin  and  Smol’yakov  (1981)  compared  sampled  polymer  using  a  near 
tangential  injector  on  a  small  (-1  meter  long)  flat  plate  (Vdovin  and  Smol’yakov,  1978)  with  the  similar 
measurements  on  their  larger-scale  test.  From  these  concentration  measurements,  they  suggested  that  near- 
tangential  injection  of  the  polymer  solution  would  be  advantageous.  However,  Petrie  et  al.  (2003)  and 
Petrie  et  al.  (2005)  employed  an  injector  with  a  near  tangential  exit,  and  the  resulting  drag  reduction  was 
not  substantially  different  when  compared  to  data  obtained  with  finite  angle  slots,  as  evidenced  by  the 
comparative  %DR  vs.  K  plots.  This  observation  is  consistent  with  those  of  Walker  et  al.  (1986)  who 
reported  that  the  developed  drag  reduction  was  not  sensitive  to  changes  in  a  slot  injector  angle  between  15 
and  25  degrees. 

Finally,  while  %DR  vs.  K  plots  are  helpful  in  organizing  and  studying  the  data  from  comparable 
experiments,  it  may  be  problematic  to  use  these  data  to  predict  far  beyond  the  parameters  of  the  original 
experiments  from  which  they  are  derived.  Injection  of  very  dilute  or  highly  concentrated  solutions,  or 
injection  at  very  high  fluxes  in  excess  of  \0QS,  say,  can  produce  different  flows  near  the  injector  and 
mixing  rates  in  the  initial  region,  especially  since  very  high  injection  fluxes  may  result  in  polymer  solution 
being  driven  from  the  surface.  Findings  for  the  detailed  mixing  processes  associated  with  the  initial  and 
transition  regions  are  discussed  next. 

6.  Evolution  of  the  near-wall  polymer  concentration 

6.1 .  Near-wall  polymer  concentration  profiles 

Measurements  of  the  polymer  concentration  were  made  in  the  near-wall  region  at  three  stream-wise 
locations  using  pulsed  laser  induced  fluorescence  (PL1F).  Examples  of  wall-normal  profiles  of  the  polymer 
at  the  first  two  PLIF  measurement  stations  (X  -  X,  =  0.64  m  and  4.62  m)  are  shown  in  Figure  9  for  WSR- 
301,  6.65  ms'1,  and  Q,  =  lOQ^.  The  sharpest  concentration  gradients  near  the  wall  occur  at  the  first 


measurement  station.  Note  that  the  peak  polymer  concentration  does  not  occur  at  the  wall,  but  a  few 
hundred  microns  away.  This  phenomenon  was  only  observed  at  the  Station  1  at  the  lowest  test  speed. 
While  an  unintended  artifact  of  the  optical  train  may  be  suspected  for  this  anomaly,  this  was  not  observed 
in  the  calibrations  performed  with  uniform  polymer/dye  concentrations.  Moreover,  this  phenomenon  has 
also  been  observed  by  Brungart  et  al.  (1991).  They  found  peak  concentrations  occurred  at  10  <  y+  <  50 
from  the  wall  0.12  m  from  the  injector  at  4.6  ms'1.  However,  they  ascribe  this  result  to  the  varying  index  of 
refraction  with  polymer  concentration.  At  the  second  and  third  measurement  stations,  the  concentration 
profiles  are  much  smoother  and  nearly  uniform  over  the  measurement  volume  0  <  y  <  2  mm.  For  the 
remaining  diffusion  analysis,  the  maximum  concentration  will  be  used  which  is  often,  but  not  necessarily 
the  wall  concentration. 

Vdovin  and  Smofyakov  (1978  and  1981),  Fontaine  et  al.  (1996),  Petrie  et  al.  (1996),  and  Petrie  et  al. 
(2005)  measured  concentration  across  the  boundary  layer  and  discuss  the  evolution  of  the  diffusion  zone 
50%  thickness  (A,)  relative  to  the  boundary  layer  thickness,  A/^9,  as  a  function  of  downstream  distance 
from  the  point  of  injection.  Sommer  and  Petrie  (1992)  also  discuss  the  evolution  of  the  concentration 
profile,  C(y)/cM,  within  the  different  mixing  zones,  where  cM  is  the  maximum  concentration  in  the  near-wall 
region.  The  concentration  measurements  of  the  present  experiments  do  not  extend  over  the  entire  thickness 
of  the  polymer  diffusion  layer,  but  are  limited  to  between  1%  and  10%  of  the  Newtonian  boundary  layer 
thickness,  making  the  typical  measured  thickness  much  less  thanA/<^9. 

6.2.  Near-wall  concentration  vs.  K 

The  scaling  of  polymer  concentration  as  a  function  K  is  provided  in  Figure  10.  As  expected,  the  near  wall 
concentration  is  reduced  as  K  is  reduced  (e.g.  as  the  distance  from  the  injection  location  increases).  This 
scaling  produces  a  reasonable  collapse  of  cM,  but  the  scatter  increases  if  the  data  are  normalized  with  the 
injection  concentration,  C,.  Most  of  the  data  collected  at  the  second  and  third  measurement  stations  are  in 
the  transitional  and  final  mixing  zones,  while  a  few  data  points  are  collected  in  the  initial  region.  (The  data 
for  WSR-N60K  also  includes  data  presented  by  Petrie  et  al .,  2005).  There  is  considerable  scatter  in  the 
data,  but  the  data  consistently  suggest  that,  with  increasing  speed,  cM  is  reduced  for  a  given  value  of  K. 


The  experimental  conditions  of  the  present  study  are  comparable  to  those  of  previously  reported 
experimental  studies  of  Fruman  and  Tulin  (1976),  Vdovin  and  Smol'yakov  (1978)  and  (1981),  Fontaine  et 
al.  (1992),  and  Petrie  et  al.  (2005),  presented  in  Table  3.  As  in  the  present  work,  these  studies  injected 
polymer  solutions  mixed  with  dye,  and  the  polymer  concentration  is  inferred  from  the  dye  concentration. 
However,  in  the  first  three  references,  the  near  wall  flow  was  physically  drawn  from  the  test  model  and 
piped  to  a  concentration  measurement  system.  Fruman  and  Tulin  (1976)  discuss  the  proper  selection  of 
sampling  flux,  since  it  is  important  that  only  the  near  wall  fluid  be  sampled.  However,  even  if  the  proper 
boundary  layer  flux  is  prescribed,  it  is  possible  to  erroneously  sample  the  polymer  solution,  since  the 
solution  is  highly  non-Newtonian  (Lotto  and  El  Riedy,  1981).  Hence,  these  indirect  measurements  of  wall 
concentration  are  potentially  less  certain  than  those  resulting  from  direct  in  situ  PLIF  interrogation. 
Moreover,  it  is  unclear  how  Fruman  and  Tulin  (1976)  managed  the  build-up  of  background  polymer  and 
dye  in  their  experiment.  Nevertheless,  we  will  include  these  measurements  in  our  comparison  with  the 
present  data  as  they  illustrate  the  overall  mixing  process. 

Figures  1 1  and  12  present  the  dimensional  and  normalized  development  of  the  near- wall  concentration  as  a 
function  of  K.  Three  regimes  of  polymer  mixing  can  be  identified.  The  development  region  is, 
characterized  by  a  wall  concentration  that  is  decreasing  as  approximately  K02.  At  a  value  of  1  x  10 7  <  K  < 
2  x  1  O'7,  the  wall  concentration  of  the  polymer  begins  to  decrease  at  a  must  faster  rate  of  K21.  This  is  the 
transitional  mixing  regime  where  the  polymer  is  being  mixed  across  the  log-layer  of  the  TBL.  Once  fully 
mixed  across  the  TBL,  the  polymer  is  much  more  diluted  and  the  level  of  drag  reduction  is  much  reduced. 
In  the  final  mixing  region,  the  near  wall  concentration  reduces  at  the  rate  of  TBL  growth,  The  data  in 
Figure  12  are  scattered,  but  the  overall  trends  and  dilution  rates  are  remarkably  similar.  The  data  of 
Fruman  and  Tulin  (1976)  lie  outside  the  general  envelope  of  the  data  from  more  recent  experiments,  but  it 
should  be  noted  that  these  experiments  employed  relatively  high  polymer  fluxes  (10  <  Q/Qs  <  30)  and  high 
injection  velocities  (10%  to  40%  of  the  free-stream  speed)  compared  to  the  conditions  of  the  other 
experiments  shown  in  Table  4. 


6.3.  Mixing  in  the  development  region 

The  development  region  was  not  commonly  measured  in  the  present  experiments,  since  its  length  was  often 
less  than  the  -0.6  meters  from  the  injector  to  the  location  of  the  first  measurement  station.  However,  we 
can  infer  the  development  length  through  examination  of  both  the  %DR  vs.  K  plots  and  the  cM  vs.  K  plots. 
Both  show  the  start  of  the  precipitous  drop  in  %DR  and  cM  with  decreasing  K  ( e.g .  increasing  distance  from 
the  injector  location)  between  1x10 1  <  K  <  2xl07,  with  the  length  being  somewhat  longer  at  lower  speeds. 
This  is  consistent  with  the  development  of  the  50%  thickness  of  the  concentration  layer  described  by  Petrie 
et  al.  (2003)  and  Petrie  et  al.  (2005).  They  reported  that  the  measured  thickness  of  the  concentration  layer 
was  typically  much  less  than  X/^  =  0.1,  corresponding  to  the  development  region  for  1x10 7  <  K  <  2x10 7. 
For  K  <  lxl O'7  the  concentration  thickness  rapidly  increases,  vVith  X/^  >  0.5  for  K  <  10*.  The  rate  of 
polymer  mixing  is  variable,  but  roughly  corresponds  to  a  rate  of  K02,  although  the  slope  varies 
considerably. 

The  length  of  the  development  region  depends  on  the  mixing  rate  of  the  near- wall  polymer  flow  formed 
downstream  of  the  injector.  If  we  consider  the  flow  through  the  Newtonian  boundary  layer,  a  flux  of  Q/Qs 
=  1  to  10  exists  in  a  region  of  the  fluid  bounded  by  y+  =  10  to  60,  and  the  average  velocity  in  the  layer 
would  be  5 ux  to  12wr  This  would  make  the  initial  layer  have  thicknesses  on  the  order  of  100  microns  for 
the  current  experiment,  and  this  is  consistent  with  the  near  wall  profiles  of  such  as  that  shown  in  Figure  10. 
In  this  region  the  rate  at  which  the  near  wall  polymers  mix  into  the  TBL  is  the  slowest  immediately 
downstream  of  the  injection  location.  Poreh  and  Zsu  (1972)  scale  the  mixing  process  as  a  diffusion  layer 
with  a  characteristic  length,  LD  =  Uh2IDDy  where  U  is  a  velocity  scale,  h  is  the  layer  thickness,  and  Dd  is  a 
diffusivity.  They  show  that  mixing  due  to  molecular  diffusivity  alone  is  much  too  low  to  be  consistent  with 
the  observed  mixing  rates,  and  they  suggest  that  the  actual  diffusivity  is  that  of  the  near-wall  turbulent  flow. 
This  assumption,  however,  implies  a  much  higher  rate  of  mixing  and  shorter  development  region  than  is 
typically  observed. 

The  injected  polymer  forms  a  near-wall  sheet  in  only  the  most  general  sense.  Recall  that  the  injected 
polymer  solutions  have  relaxation  times  on  the  order  of  10  milliseconds.  The  polymer  need  not  exit  the 


injector  as  a  planar  sheet,  but  is  drawn  from  the  slit  as  a  result  of  the  very  high  elasticity  of  the  solution. 
The  flow  then  strains  the  polymer,  forming  filaments,  which  will  be  discussed  in  Section  6.5.  The 
relaxation  times  of  the  highly  concentrated  injected  solutions  are  larger  than  the  integral  timescales  of  the 
flow,  Sgg/U  ~  10'3  5  near  the  injector.  Hence,  while  the  mean  shear  level  near  the  wall  is  quite  large,  the 
high-concentration  polymer  solutions  will  not  be  easily  deformed  by  the  relatively  high-frequency  velocity 
fluctuations  in  the  incoming  TBL.  Instead,  the  injected  polymer  will  be  strained  by  the  mean  near-wall 
shear  flow. 

Examination  of  the  K  vs.  %DR  plots  indicates  that  the  length  of  the  development  region,  LD,  is 
approximately  scaled  with  a  K  value  of  KD  ~  10'7,  and  the  present  data  suggest  that  increasing  speed 
decreases  LD  and  KD.  Nevertheless,  the  observation  that  the  /^-scaling  is  somewhat  successful  at  scaling  the 
length  of  the  development  region  is  revealing,  since  this  supports  the  notion  that  the  mixing  processes  in 
the  initial  region  is  grossly  scaled  with  the  outer  variables  of  the  TBL.  We  will  return  to  the  discussion  of 
the  mixing  in  the  development  region  during  the  discussion  of  polymer  filamentation. 

6.4.  Mixing  in  the  transitional  and  final  regions 

The  transitional  region  extends  from  roughly  from  2xl0‘7  >  K  >  2x10^.  The  mixing  rate  is  not  constant 
over  this  region,  but  roughly  scales  as  K27  (compared  with  the  exponent  of  2.5  reported  by  Petrie  et  al. 
(2005),  for  WSR-N60K.)  Again,  it  is  helpful  to  compare  using  K  vs.  %DR  plots.  Over  this  range  of  AT,  the 
drag  reduction  drops  from  its  peak  to  around  10  to  20  %DR. 

Vdovin  and  Smol'yakov  (1978  and  1981),  Fontaine  et  al.  (1992)  suggest  that  the  polymer  mixing  process 
can  be  scaled  for  a  given  polymer  type  through  determination  of  the  “concentration  decay  length”,  L,  such 


where  L  is  the  e-folding  distance  ( i.e .  the  length  from  the  injection  location  to  where  the  near-wall 
concentration  falls  to  ex  =  37%  of  the  injected  solution  concentration).  It  has  been  supposed  that  L  is  a 
function  of  the  injection  parameters  and  the  free-stream  speed: 


l=KtQ£ \ 

pu 


(6.2) 


where  KT  is  a  constant  that  is  fixed  for  a  given  polymer-solvent  combination  and  molecular  weight.  Data 
presented  by  Vdovin  and  Smol'yakov  (1978  and  1981),  Fontaine  et  al.  (1992),  and  Petrie  et  al.  (1996) 
indicate  that  4.9x1 06  <  KT<  6.25x1 06  for  WSR-301  in  water.  If  KT  is  solely  a  function  of  polymer  type, 
then  K  will  be  sufficient  to  scale  the  mixing  processes.  For  this  to  be  the  case,  the  length  of  the 
development  region  must  be  small  compared  to  L,  since  we  are  taking  as  our  length  the  distance  from  the 
injector  and  not  from  the  end  of  the  development  region.  Also,  we  must  assume  that  the  mixing  rate  in 
linear  with  free-stream  speed.  The  plot  of  cJCj  vs.  K  in  Figure  12  shows  an  exponential  decay  in 
concentration  after  the  development  region,  but  the  scatter  in  the  data  suggests  that  the  decay  length  is  not 
completely  scaled  with  K.  (Note  that  the  data  from  Vdovin  and  Smol'yakov  (1981)  is  re-plotted  using  a 
constant  value  of  L  derived  from  their  best  fit  for  KT.  We  believe  that  the  cJCj  data  presented  in  their 
manuscript  have  been  normalized  by  the  specific  value  of  L  derived  from  each  individual  experiment. 
However,  we  are  forced  to  use  their  best-fit  value  of  L  to  re-plot  their  data.  Therefore,  the  scatter  in  their 
data  due  to  variation  in  L  is  not  represented.  Figure  1  of  their  manuscript  suggests  a  variation  in  L  of 
±20%). 


In  the  transitional  mixing  region,  the  wall  concentrations  decrease  from  100’s  of  wppm  to  around  10  wppm , 
but  significant  turbulence  modification  is  still  possible  with  concentrations  as  low  as  a  few  wppm.  Finite 
levels  of  drag  reduction  are  still  observed.  The  mixing  rate  decreases  in  the  final  region,  where  the  polymer 
has  been  diluted  to  less  than  10  wppm  and  is  largely  mixed  throughout  the  TBL,  at  K  ~  10"*.  The  rate  of 
concentration  reduction  scales  with  the  growth  of  the  TBL  due  to  entrainment,  with  a  growth  rate  of  ~K718, 
which  corresponds  to  the  entrainment  rate  of  a  high  Reynolds  number  TBL  with  little  or  no  friction  drag 
reduction.  By  K  ~  10"9,  the  drag  reduction  has  ceased,  and  the  concentration  has  fallen  to  less  than  1  wppm. 


6.5.  Filamentation  of  high  concentration  polymer  solutions 
The  PLIF  measurements  of  the  near-wall  concentration  field  often  exhibited  significant  stratification.  A 
layer  of  injected  polymer  was  often  observed  near  the  wall,  but  the  distribution  of  polymer  was  not 
necessarily  uniform  and  continuous.  Particular  images  would  reveal  the  near  absence  of  polymer  near  the 
wall  with  local,  high  concentration  strands  of  polymer  away  from  the  wall.  Vlachogiannis  and  Hanratty 
(2004)  observed  similar  processes  of  polymer  filamentation  in  their  study  of  wall-injected  polymer  in  a 
channel  flow.  They  concluded  that  the  presence  of  polymer  filamentation  resulted  in  enhanced  drag 
reduction,  and  this,  in  turn  may  be  due  to  the  formation  of  polymer  aggregates,  as  previously  suggested  by 
Cox  et  al.  (1974).  As  discussed  above,  it  is  likely  that  there  is  substantial  aggregation  of  the  injected 
polymer  solutions.  The  high  elasticity  of  the  injected  solutions  leads  them  to  be  drawn  from  the  injector, 
which  will  introduce  three-dimensionality  to  the  near-wall  polymer  concentration  distribution. 

Examination  of  the  PLIF  images  shows  that  for  K  >  2.4  x  10'7,  the  concentration  field  is  very  stratified, 
with  large  regions  of  nearly  pure  water  and  high  concentrations  of  polymer  filaments.  Figure  13  shows 
three  sample  images.  As  the  polymer  diffuses  from  the  filaments  into  the  surrounding  solvent,  the 
background  concentration  of  the  polymer  increases.  Between  5  X  10 *<  K  <  2.4  X  10‘7,  the  near  wall  region 
consists  of  polymer  filaments  with  a  background  of  somewhat  homogeneous  polymer  solution  (Figure  14). 
In  this  region  the  peak  drag  reduction  has  occurred  and  the  mixing  into  the  log-layer  has  commenced. 
Finally,  with  K  <  5  x  1  O'8,  the  filaments  have  mixed  into  the  flow  and  the  concentration  fields  are  nearly 
uniform  throughout  the  field  of  view.  The  rheology  of  these  highly  stratified  mixtures  is  potentially  much 
more  complicated  than  that  of  the  uniform  mixtures,  but  the  relaxation  time  of  the  filaments  will  be  on  the 
order  of  that  for  the  injected  polymer  solutions.  As  discussed  in  Section  6.3,  the  stretching  and  mixing  of 
these  filaments  will  then  take  place  over  outer  flow  length  and  time-scales  that  result  in  the  development 
region,  until  the  flow  is  sufficiently  diluted  and  mixed  to  near-wall  concentrations  less  than  100  wppm. 


7.  Relationship  between  polymer  type,  concentration,  and  the  resulting  drag  reduction 


7.7.  Development  of  the  Weissenberg  number  with  flow  conditions 
We  next  consider  how  the  observed  drag  reduction  can  be  related  to  the  polymer  solution  rheology  and  the 
near-wall  flow  of  the  TBL.  The  onset  of  drag  reduction  occurs  when  the  near-wall  shear  rate  exceeds  y  *  . 

We  can  define  a  Weissenberg  number  for  onset  of  drag  reduction, 

2 


We*  =  -^~ 
vy  * 


(7.1) 


Note  that  We*  >  1  because  onset  scales  on  the  magnitude  of  the  fluctuating  velocity  gradient.  These 
fluctuations  are  about  an  order  of  magnitude  smaller  than  the  turbulent  scales  used  in  (7.1)  (Dimitropoulos 
ei  al.,  1998).  The  We*  onset  condition  can  be  compared  to  the  Weissenberg  number  based  on  the 
relaxation  time  of  the  polymer  at  the  particular  conditions  of  the  experiment, 


(7.2) 


where  6  is  the  value  provided  by  Kalashnikov  (1998)  for  values  above  the  Zimm  time,  which  is  the 
minimum  value.  Experiments  in  dilute  polymer  solutions  indicate  that  onset  typically  occurs  when  the 
deformation  rate  in  the  flow  exceeds  a  value  necessary  to  stretch  the  long-chain  polymers,  and  this 
corresponds  to  the  onset  value  of  We\  resulting  in  a  We*  that  is  typically  much  greater  than  unity 
(Sreenivasan  and  White,  2000;  Housiadis  and  Beris,  2003).  Table  5  shows  the  spatial  evolution  of  We*  for 
WSR-301  for  ranges  of  K  that  correspond  to  the  development,  transitional,  and  final  mixing  regions  for  the 
intermediate  speed  of  13.2  ms'.  (Note  that  We*  at  any  downstream  location  on  the  plate  is  a  function  of  the 
local  polymer  concentration.)  At  these  high  values  of  Weissenberg  number,  the  individual  polymer 
molecules  can  achieve  significant  extension.  For  We*  »  We*,  polymer  modification  of  the  turbulence 
will  become  increasingly  insensitive  to  further  increases  in  We*.  This  inequality  appears  satisfied  for  most 
experiments  except  the  dilute  flows  of  WSR-N60K  at  6.65  ms'.  These  results  are  consistent  with  the 
findings  of  Housiadis  and  Beris  (2003)  who  suggest  that  at  these  values  of  We*,  the  flow  will  be  in  the  low 
to  moderate  drag  reduction  regime.  The  flows  with  WSR-301  and  WSR-308  both  have  sufficiently  high 


We\  and  these  flows  are  in  the  high  drag  reduction  regime  for  all  but  the  most  dilute  ( i.e .  downstream) 
portions  of  the  flow. 


Finally,  the  potential  for  polymer  degradation  can  be  evaluated  using  (4.7)  with  the  definition  of  the 
degradation  Weissenberg  number,  WeD 


WeD  = 


(7.3) 


Significant  levels  of  polymer  degradation  are  expected  when  WeD  »  1,  which  is  likely  to  occur  with  the 
present  flow. 


7.2.  Comparison  of  observed  and  expected  levels  of  drag  reduction  with  varying  Mw 
With  the  data  presented  above,  we  can  construct  the  universal  drag  reduction  curves  for  the  flow  over  the 
plate.  These  curves  can  be  used  to  determine  the  intrinsic  drag  reduction  and  concentration  of  a  particular 
polymer,  quantities  that  were  introduced  for  pipe-flows  by  Virk  et  al.  (1967).  Figure  15  shows  %DRIcM  as 
a  function  of  cM  for  the  three  polymer  molecular  weights.  The  intrinsic  drag  reduction  [%DR/cM]  and 
intrinsic  concentration  [C]  was  determined  as  a  result  of  a  best  fit  of  the  data  to  a  function  of  the  form 


%DR  /  cM  = 


[%£/?/ Cj 

1  +  CM/[C] 


(7.4) 


Table  6  presents  the  values  for  each  molecular  weight  at  each  test  speed  and  an  average  curve  fitted  to  the 
data  from  of  all  of  tested  speeds.  The  data  are  scattered  at  the  lowest  concentrations,  but  some  trends  are 
apparent.  First,  the  lowest  molecular  weight  polymer,  WSR-N60K,  has  much  lower  intrinsic  drag 
reduction  values  than  WSR-301.  The  average  intrinsic  drag  reduction  of  WSR-N60K  is  -5  times  less  than 
that  of  WSR-301.  The  average  intrinsic  drag  reduction  of  WSR-308  is  -1.2  times  that  of  WSR-301. 


The  relationship  between  the  intrinsic  drag  reduction  and  the  underlying  polymer  rheology  is  not  readily 
apparent;  however,  qualitatively  we  see  that  this  quantity  is  more  influenced  by  the  drag-reducing  polymer 


than  the  flow  speed.  The  existence  of  the  intrinsic  drag  reduction  from  the  functional  form  of  7.4  implies 
dilute  solution  behavior,  as  evidence  by  the  low  concentration  plateaus  in  Figure  15.  Thus,  we  might 
expect  that  [%DR/cM]  would  scale  with  the  Zimm  relaxation  time  (or  some  multiple  thereof),  but 
examination  of  Table  2  shows  that  the  6Z  for  WSR-301  is  ~2.7  times  that  of  WSR-N60K,  and  Bz  for  WSR- 
308  is  ~3.4  times  that  of  WSR-301.  Thus,  consideration  of  the  Zimm  time  alone  is  not  sufficient  to  scale 
these  results.  Use  of  the  Kalashnakov  relaxation  time  introduces  two  further  complications.  First,  this 
relationship  is  valid  for  relatively  high  concentrations.  Tirtaatmadja  et  al.  (2006)  discuss  the  influence  of 
molecular  weight  and  concentration  on  the  relaxation  time  for  PEO  solutions  in  good  solvents.  They 
showed  that,  for  solutions  with  C/c*  <  10'2,  the  relaxation  time  approaches  the  Zimm  time.  In  regimes 
where  the  concentration  is  C/c*  >  10*2,  we  can  use  the  Kalashnikov  relaxation  time  to  compare  the 
performance  of  the  polymer  solutions.  This  corresponds  to  concentrations  of  approximately  5  wppm  and 
above.  Second,  dK  ~  Vc  which  would  not  lead  to  a  constant  value  of  the  intrinsic  drag  reduction  in  the 
limit  as  C  — >  0,  a  result  that  appears  inconsistent  with  the  data  from  Figure  15. 

Virk  et  al.  (1967)  reported  a  phenomenological  relationship  between  the  intrinsic  drag  reduction  and 
intrinsic  viscosity,  where 

[%DR  I cM\~ [77]/  ln(We*)  (7.5) 

The  changing  free-stream  speed  will  lead  to  slight  variation  in  We*.  For  example,  changing  the  speed  from 
6.65  ms'1  to  13.2  ms'1  increases  ln(We*)  by  ~1.3  for  all  three  polymers.  Similarly,  further  changing  the 
speed  from  13.2  ms'1  to  19.9  ms'1  increases  ln(We*)  by  ~1.2.  Combining  these  results  with  Eq.  4.1,  the 
scaling  of  Eq.  7.5  implies  that  WSR-301  would  be  ~4.6  as  effective  at  WSR-N60K,  and  WSR-308  would 
be  ~6.3  times  as  effective  as  WSR-301.  This  empirical  scaling  appears  to  capture  the  difference  in  the 
intrinsic  drag  reduction  between  WSR-301  and  WSR-N60K.  But,  the  intrinsic  performance  of  the  dilute 
WSR-308  solutions  is  much  less  than  predicted.  We  argue  that  an  important  contributing  factor  to  this 
discrepancy  is  the  process  of  flow-induced  polymer  degradation. 


The  intrinsic  drag  reduction  tends  to  decrease  with  increasing  speed,  and  these  changes  are  more 
pronounced  with  increasing  molecular  weight.  Examination  of  WeD  for  the  13.2  ms1  conditions  in  Table  5 
indicates  that  all  the  polymers  are  potentially  susceptible  to  degradation,  with  WSR-308  being  the  most 
likely  to  degrade.  Vanapalli  etal.  (2005)  has  shown  that  high  molecular  weight  PEO  solutions  like  WSR- 
308  undergo  significant  degradation  in  turbulent  shear  flow  for  WeD  »  1.  Continued  exposure  to  flows 
with  sufficient  stretch  for  scission  results  in  continual  reduction  in  Mw,  until  the  polymer  molecules  are 
sufficiently  shortened  that  they  become  scission  resistant  at  the  given  deformation  rate  of  the  flow. 
Degradation  of  WSR-308  will  result  in  a  polymer  solution  that  behaves  more  like  WSR-301,  as  seen  in  the 
present  study. 

The  observation  that  the  intrinsic  drag  reduction  decreases  with  speed  is  consistent  with  the  presence  of 
increased  degradation  at  high  speeds  and  higher  Mw  for  a  given  speed.  Thus,  under  the  conditions  of  the 
present  study,  the  flow  in  the  TBL  will  continually  degrade  polymer  molecules,  reducing  the  mean 
molecular  weight  of  the  polymers  in  solution.  (Note  that  scission  of  the  molecules  will  increase  the  number 
of  lower  Mw  molecules,  but  maintain  a  constant  mass-fraction  concentration).  Vanapalli  et  al.  (2006) 
analyzed  the  scaling  of  polymer  scission  in  wall-bounded  turbulent  shear  flow  and  show  that  the  bulk  of  the 
polymer  degradation  takes  place  in  the  log-layer  of  the  TBL,  since  this  is  where  the  appropriate, 
combination  of  volume  flux  and  shear  rate  is  highest.  Consequently,  the  rate  at  which  the  polymers 
degrade  is  expected  to  be  a  function  of  numerous  flow  parameters,  including  both  the  inner  and  outer  scales 
of  the  TBL,  and  the  distribution  of  polymer  within  the  TBL.  The  data  presented  here  suggest  that  the 
stream-wise  extent  of  the  test  model  used  in  the  current  study  was  sufficiently  long  to  reveal  the  presence 
of  significant  polymer  degradation  for  the  higher  molecular  weight  polymers.  Note  that  the  most  dilute 
polymer  solutions  occurred  farthest  downstream  from  the  injector,  and  the  polymers  solutions  there  have 
had  the  longest  exposure  to  the  turbulent  flow. 

Examining  the  intrinsic  drag  reduction  of  WSR-308  reveals  that  [%DR/cM]  is  much  higher  at  the  lowest 
speed  and  is  about  twice  that  of  WSR-301  under  equivalent  conditions  (although  the  data  are  sparse). 
These  low-speed  results  illustrate  that  WSR-308  can  achieve  a  higher  intrinsic  drag  reduction  if  the  level  of 


degradation  is  reduced.  However,  as  the  flow  speed  is  increased,  the  rate  of  degradation  becomes 
substantial,  and  the  drag  reduction  achieved  by  the  degraded  WSR-308  approaches  that  of  WSR-301. 

8.  Conclusions 

A  set  of  experiments  has  been  conducted  at  the  Large  Cavitation  Channel  (LCC)  in  Memphis,  Tennessee  to 
investigate  the  phenomena  of  friction  drag  reduction  at  large  scales  and  high  Reynolds  numbers  in  a 
controlled  test  environment.  Drag  reduction  was  achieved  by  injecting  aqueous  solutions  of  polyethylene 
oxide)  (PEO)  of  three  different  molecular  weights  into  the  nearly  zero-pressure-gradient  turbulent  boundary 
layer  on  a  smooth  flat-plate  test  model.  The  measurements  covered  more  than  nine  meters  of  flow 
development  length  downstream  of  the  injector  at  free  stream  flow  speeds  up  to  20  ms"'. 

With  sufficiently  large  mass  flux  of  injected  polymer,  maximum  drag  reduction  (MDR)  was  reached  over  a 
limited  extent  of  the  test  model  surface.  However,  a  plateau  at  MDR  was  not  observed.  Instead,  a 
development  region  near  the  injector  produced  a  broad  peak  in  drag  reduction  where  MDR  was  nearly 
achieved,  then  the  continued  process  of  mixing  led  to  a  steady  decrease  in  the  level  of  drag  reduction. 
Increasing  polymer  molecular  weight  resulted  in  higher  levels  of  drag  reduction,  as  did  increasing  the  net 
polymer  flux  injected  into  the  TBL.  These  results  were  collected  at  high  Reynolds  numbers  and  large 
physical  scales,  and  they  are  broadly  consistent  with  previously  reported  results  acquired  on  smaller  test 
models  and  at  lower  speeds. 

The  process  of  polymer  mixing  in  the  TBL  can  be  broadly  divided  into  three  regimes:  the  development 
region,  the  transitional  region,  and  the  final  region.  The  development  region  occurs  just  downstream  of  the 
injection  location.  The  highly  viscoelastic  polymer  solutions  are  stretched  and  mixed  into  the  TBL, 
forming  a  stratified  mixture  of  polymer  and  solvent.  As  this  mixing  process  occurs,  the  friction  drag 
reduction  increases  to  its  peak  value.  The  rate  of  polymer  mixing  in  the  development  region  is  relatively 
slow.  However,  as  the  polymer  begins  to  mix  across  the  TBL  in  the  transitional  regime,  the  mixing  rate 
dramatically  increases.  Mixing-induced  dilution  of  the  polymer  reduced  the  level  of  friction  drag 
reduction.  And,  in  the  final  region,  polymer  dilution  is  due  to  entrainment  of  free-stream  liquid  into  the 


TBL.  At  best  these  mixing  processes  are  adequately  scaled  with  the  /^-Factor,  a  parameter  that  includes  the 
injected  polymer  flux,  free-stream  speed,  and  distance  from  the  injection  location.  By  injecting  high- 
concentration  polymer  solution,  a  relatively  long  development  region  was  produced.  However,  the  length 
of  this  high  drag  reduction  never  extended  over  the  entire  test  model. 

The  drag  reducing  performance  of  the  polymers  increased  with  increasing  weight  concentration  and 
molecular  weight  of  the  polymer  solution.  However,  the  level  of  drag  reduction  for  a  given  near  wall 
polymer  concentration  was  reduced  with  increasing  speed.  And,  the  highest  molecular  weight  solutions 
produced,  on  average,  levels  of  drag  reduction  that  were  only  modestly  higher  than  solutions  with  one-half 
the  molecular  weight.  These  observations  are  consistent  with  the  presence  of  flow-induced  degradation,  a 
process  that  was  expected  to  occur  based  on  the  known  molecular  scission  potential  of  the  PEO  solutions 
used.  Attempts  to  scale  the  observed  intrinsic  drag  reduction  with  changes  in  PEO  molecular  weight  were 
not  successful,  suggesting  that  the  solution  intrinsic  viscosity  alone  is  not  sufficient  to  scale  the  observed 
levels  of  drag  reduction  in  the  diluted  polymer  solutions.  We  hypothesize  that  the  presence  of  significant 
polymer  degradation  will  lead  to  a  failure  of  such  simple  scalings.  Moreover,  the  relationship  between  the 
solution  relaxation  time,  concentration,  and  resulting  turbulence  modification  may  be  too  complex  for  such 
simple  scaling  ideas,  especially  in  regions  of  high  drag  reduction.  Future  study  is  needed  to  resolve  this 
issue. 
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Table  1:  Single  phase  free- stream  velocity  U ,  Reynolds  number  Re*,  99%  boundary  layer  thickness  5*,,  momentum 
thickness  0,  acceleration  parameter  A",  skin  friction  coefficient  Cp  the  friction  velocity  uv  the  viscous  length  and  the 
surface  roughness  parameter  k \  at  the  downstream  locations  of  the  injector,  the  six  skin-friction  measurement  balances, 

and  the  three  PLIF  measurement  stations. 
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86 

7.8 

3.6 

0.0017 

0.39 

2.6 

0.2 

20.22 

217 

88 

7.5 

2.4 

0.0016 

0.57 

1.8 

0.2 

Table  2:  Polymer  Solution  Properties. 


PEO 

Polymer 

Mw 

(g  mole1) 

[Til 

(cm1  g') 

c* 

(wppm) 

4/C 

(wppm1) 

&Z 

(s) 

ejc"2 

(s) 

Y  ♦ 

(s') 

Yd 

(s') 

WSR-N60K 

2.3  x  106 

1.1  x  103 

8.7  x  102 

4.4  x  102 

4.6  x  10" 

1.30 

1.5  x  103 

2.4  x  104 

WSR-301 

4x  106 

1.8  x  103 

5.7  x  102 

8.1  x  102 

1.2  x  10 3 

2.17 

8.4  x  102 

7.1  x  103 

WSR-308 

8  x  106 

3.0  x  103 

3.3  x  10* 

1.7  x  103 

4.2  x  IQ'3 

3.58 

4.2  x  102 

1.5  x  103 

Table  3:  Experimental  setup  and  conditions  of  previous  experiments  reporting  drag  reduction  vs.  K-Factor  for  line- 
source  injection  into  a  turbulent  boundary  layer  using  WSR-301. 


TBL 

Thickness 

at 

Injector 

(mm) 

Maximum 
Measurement 
Distance 
from  the 
Injector 
(m) 

Speed 

Range 

(ms1) 

Cone. 

Range 

(wppm) 

Q/Q. 

Balance 

Type 

Injector 

Type 

Injector 

Span 

(%of 

total) 

Xo 

K=2xl07 
U*  6.7 
ms-1 
(m) 

XD 

Jfs2xl0’ 
U  =  20 
ms-1 
(m) 

Vdovin  and 
Smol'yakov 
(1981) 

8  to  50 

4.375 

2.5  to 
10 

250  to 
4000 

1.1  to 
23 

Local 
(20  mm  x 
20  mm) 

20  degree  slot 
(0.4  x  140 
mm  or  0.8  x 
140  mm) 

87% 

0.02  to 
5.16 

0.01  to 
1.72 

Fontaine  et  al. 
(1992) 

-4* 

0.36  to 
balance 
trailing  edge 

4.5  to 
18.3 

500  and 
1000 

1 .6  to  4 

Integrated 
(318  x 
152  mm) 

25  degree  slot 
(1  mm  x  152 
mm) 

50% 

0.05  to 
0.22 

0.02  to 
0.08 

Petrie  et  al. 
(1996) 

-10  to 
-20** 

1.87 

7.6  to 
16.8 

500  to 
2000 

2.5  to 
10 

Local 
(38 x  127 
mm) 

25  degree  slot 
(2.5  mm  x 
610  mm) 

50% 

0.07  to 
1.12 

0.02  to 
0.37 

Petrie  et  al. 
(2003) 

-4* 

0.38  to 
balance 
trailing  edge 

4.57  to 
13.7 

200  to 
2000 

2  to  20 

Integrated 
(318  x 
152  mm) 

Contoured 
Plenum  Slot 

81% 

0.02  to 
2.24 

0.01  to 
0.75 

Petrie  et  al. 
(2005) 

-4* 

0.33 

6  to  18 

100  to 
500 

2  to  20 

Local 
(38  x 
127) 

Contoured 
Plenum  Slot 

81% 

0.01  to 
0.56 

0.00  to 
0.19 

Present 

Results 

16  to  20 

9.36 

6.6  to 
20 

1000  to 
4000 

2  to  10 

Local 

(152.4 

mm 

diameter) 

25  degree  slot 
(1  mm  exit 
cross  section 
x  2650  mm) 

88% 

0.11  to 
2.04 

0.04  to 
0.75 

*  reported  by  Madavan  et  al.  (1985);  **  estimated  assuming  the  origin  of  theTBL  is  the  model  leading  edge 


Table  4:  Experimental  setup  and  conditions  of  previous  experiments  reporting  near-wall  concentration  measurements 
for  line-source  injection  into  a  turbulent  boundary  layer  using  WSR-301  and  WSR-N60K  (present  results  and  Petrie  et 

a/.,  2005). 


Maximum 
Measurement 
Distance 
From  the 
Injector 
(m) 

Speed 

Range 

(ms1) 

Cone. 

Range 

(wppm) 

e/e. 

Measurement 

Method 

Injector 

Type 

Injector 
Span 
(%  of 
total) 

Fruman  and  Tulin  (1976) 

2.337 

10.65 

100  to  1000 

11.2  to 
27.3 

Sampling  Slits 

38  mm  x  0.25 

mm 

0  degree 

slot  (0.5  x 
300  mm) 

-100% 

Vdovin  and  Smol'yakov 
(1978) 

0.557 

2  to  12 

560  to  5000 

1.3  to  12.6 

Sampling 

Ports  0.5  mm 
diameter  and 
Traversed 
Sampling 
Tubes  with 
0.15  x  1.5  mm 
openings 

7  degree 
slot  (0.7  x 
120) 

80% 

Vdovin  and  Smol'yakov 
(1981) 

4.5 

2.5  to  10 

250  to  4000 

1.1  to  23 

Sampling 

Ports 

20  degree 

slot  (0.4  x 
140  mm 
or  0.8  x 
140  mm) 

87% 

Fontaine  etal.  (1992) 

-0.3 

4.5  to  18.3 

500  and 
1000 

1 .6  to  4 

PLIF 

25  degree 

slot 

(1  mm  x 
152  mm) 

50% 

Petrie  et  al.  (2005) 

-0.3 

6  to  18 

100  to  500 

2  to  20 

PLIF 

Contoured 

Plenum 

Slot 

81% 

Present  Results 

9.36 

6.6  to  20 

1000  to 
4000 

2  to  10 

PLIF 

25  degree 

slot  ( r 
mm  exit 

cross 

section  x 
2650  mm) 

88% 

Table  5:  Evolution  of  the  Weissenberg  Numbers  with  decreasing  K-Factor  at  U  =  13.2  ms1.  The  intrinsic  drag 
reduction  and  intrinsic  concentration  were  derived  from  all  the  near-wall  concentration  measurements  at  all  speeds  for  a 


given  polymer. 


WSR-N60K 

K 

CM 

(wppm) 

e 

(i) 

(ms1) 

We* 

We* 

WeD 

Development  Region 

>2x  10- 

7 

>200 

>0.02 

-0.3 

>1800 

-60 

-4 

Middle  of  the  Transitional  Region 

~6x  10- 
8 

~20 

-0.01 

-0.3 

-720 

-60 

-4 

End  of  the  Transition  Region 

~2x  10 

8 

~5 

-0.004 

-0.4 

-640 

-107 

-7 

Final  Region 

~2x  10- 

9 

<1 

-0.0005 

-0.4 

-160 

-107 

-7 

WSR-301 

Development  Region 

>2x  10' 

7 

>500 

>0.05 

-0.3 

>4000 

-107 

-12 

Middle  of  the  Transitional  Region 

-  6  x  1  O' 
8 

~20 

-0.01 

-0.3 

-1000 

-107 

-12 

End  of  the  Transition  Region 

~2x  10 

8 

-2 

-0.004 

-0.4 

-640 

-190 

-23 

Final  Region 

-  2  x  1  O' 

9 

<  1 

-0.001 

-0.4 

-480 

-190 

-23 

WSR-308 

Development  Region 

>2x  10 

7 

>500 

>0.05 

-0.3 

>9900 

-214 

-59 

Middle  of  the  Transitional  Region 

-  6  x  10* 
8 

~50 

-0.035 

-0.3 

-3150 

-214 

-59 

End  of  the  Transition  Region 

~2x  10 

8 

~5 

-0.011 

-0.4 

-1760 

-381 

-105 

Final  Region 

~2x  10 

9 

<  1 

-0.004 

-0.4 

-1600 

-381 

-105 

Table  6:  The  intrinsic  drag  reduction,  [ %DR/cM ],  and  intrinsic  concentration,  [C],  for  the  near-wall  polymer  solutions. 


The  values  for  each  speed  and  for  all  speeds  combined  are  presented.  The  fitted  curves  are  shown  in  Figure  15. 


U 

(ms1) 

[%DR/cm] 

(wppm') 

si 

WSR-N60K 

6.65 

6 

7 

13.2 

6 

8 

19.9 

4 

12 

1 

All 

5 

10 

WSR-301 

6.65 

32 

2 

13.2 

25 

2 

19.9 

17 

3 

All 

26 

2 

WSR-308 

6.65 

73 

1 

13.2 

27 

3 

19.9 

14 

5 

AD 

32 

2 

Figure  1:  Schematic  drawing  of  the  large  flat-plate  test  model  with  injection  and  measurement  locations. 
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Figure  2:  Baseline  skin-friction  coefficient,  C/t  vs.  Reynolds  Number  from  two  separate  phases  of  the  experiments 
(open  symbols  are  from  Phase  Ill  and  closed  symbols  are  from  Phase  IV)  for  U  =  6.65  (A),  13.2  (O),  and  (□)  19.9  ms  '. 
Also  shown  are  the  power-law  fit  to  the  experimental  data  (Equation  3.3,  solid  line)  and  two  standard  skin-friction 
correlations  from  White  (1991)  (Equation  3.1,  light  dashed  line)  and  Schultz-Grunow  (1941)  (Equation  3.2,  heavy 

dashed  line). 
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Figure  3:  The  %DR  vs.  downstream  distance  from  the  injector,  X-X„  illustrating  the  development  region.  The  injected 
polymer  is  WSR-301  with  C,  =  4000  wppm,  and  Q/Qs=  10  for  U  =  6.65  (O)  and  13.2  ms '(•). 
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Figure  4:  The  %DR  vs.  downstream  distance  from  the  injector,  X-X„  for  WSR-308  for  C(  =  4000  wppm  polymer 
injected  at  Q/Qs=  10  for  U  =  6.65  (O),  13.2  (•),  and  19.9  ms1  (□). 


10 


Figure  5:  The  %DR  vs.  downstream  distance  from  the  injector,  X-Xit  for  WSR-301.  (O)  U  =  6.65  ms1,  C,  =1000  wppm 
and  QJQS~  4;  (•)  U  =  6.65  ms1,  C,  =2000  wppm  and  Q,/Qs=  2;  (□)  U  =  13.2  ms1,  C,  =1000  wppm  and  Q/Qs=  4;  (■) 

U  =  13.2  ms1,  C,  =2000  wppm  and  Q/Qs  =  2. 
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Figure  6:  %DR  vs.  K  for  all  three  Mw  polymers  at  U  =  6.65  (O),  13.2  (•),  and  19.9  (□)  ms'. 
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Figure  7:  %DR  vs.  K  for  WSR-301  at  U  =  6.65  (O),  13.2  (•),  and  19.9  (□)  ms1.  Also  shown  is  the  trend  line  reported 
by  Vdovin  and  S mol* yakov  (1981)  for  flow  in  the  transitional  region  (solid  line).  Similar  slopes  are  found  at  all  three 
speeds,  although  the  K  value  at  which  the  transition  region  starts  changes  with  speed  (dashed  lines). 


%DR 


80% 


70% 

60% 

50% 

40% 

30% 

20% 

10% 


•  "6.65  m/s  WSR-301" 

•  6.65  ms1  (present  results) 

O  Vdovin  and  Smol’yakov  (1981) 
A  Petrie  etal.  (1996)  7.6  ms1 
□  Petrie  et  al.  (2005)  6  ms 1 

♦  Petrie  et  al.  (2003)  4.57  ms1 

❖  Petrie  etal.  (2003)  7.62  ms'1 

*  Fontaine  etal.  (1992)  4.5  ms 1 


A 


□ 

□ 


A 

♦  □ 


—  rf  ♦  a 


0%  I - - - - - 1 

1.0E-10  1.0E-09  1.0E-08  1.0E-07  1.0E-06 

K 


(a) 


%DR 


80% 


70% 


60% 


50% 


40% 


30% 


20% 


10% 


0% 


•  13.2  ms1  (present  results) 

O  Vdovin  and  Smol’yakov  (1981) 
▲  Petrie  etal.  (1996)  10.7  ms1 

*  Petrie  etal.  (1996)  13.7  ms 1 
□  Petrie  etal.  (2005)  12  ms1 

♦  Petrie  etal.  (2003)  10.7  ms1 

❖  Petrie  era/.  (2003)  13.7  ms1 

-  Fontaine  et  al.  ( 1 992)  9. 1 4  ms1 
4  Fontaine  etal.  (1992)  13.7  ms1 


1.0E-10 


1.0E-09 


1 .0E-08 
K 


1.0E-07 


£ 


1.0E-06 


(b) 


%DR 


80% 


70% 


60% 


50% 


40% 


30% 


20% 


10% 


0% 


•  19.9  ms1  (present  results) 

A  Petrie  etal.  (1996)  16.8  ms1 
□  Petrie  etal.  (2005)  18  ms  ' 

*  Fontaine  etal.  (1992)  18.3  ms1 


□  1 

,*a  » 

B.  * 

3  CP 


*  S 

□ 


•  r 


.  .  •»»**  */ 

— ® - * — +-*-4 — *-*A- 


▲  n*  v 


# 


1.0E-10 


1.0E-09 


1.0E-08 

K 


1.0E-07 


10E-06 


(C) 

Figure  8:  vs.  %D/?  from  the  present  experiments  and  those  reported  previously  for  WSR-301.  Data  is  segregated 
into  (a)  low,  (b)  medium,  and  (c)  high  free-stream  speeds.  The  experimental  descriptions  and  conditions  are  provided  ’ 

in  Table  3. 
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Figure  9:  Near-wall  concentration  profiles  of  WSR-301  at  X-X,  =  0.64  m  (a)  and  4.62  m  (b)  for  U  -  6.65  ms  1  and 

Q/Q*=  io. 
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Figure  10:  The  maximum  near  wall  concentration,  cM,  vs.  K  for  (a)  WSR-N60K,  (b)  WSR-301,  and  (c)  WSR-308  at  U 
=  6.65  (O),  13.2  (•),  and  19.9  (□)  ms1.  The  data  for  WSR-N60K  from  Petrie  etal.  (2005)  are  also  included  (X). 
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Figure  11:  cM  vs.  K  for  WSR-301  of  the  present  experiment  (U  =  6.65  (•),  13.2  (■),  and  19.9  (A)  ms'  with  previously 
reported  results;  Fruman  and  Tulin  (1976)  with  C,  =  100  wppm  (O),  500  wppm  (A),  and  1000  wppm  (□);  Vdovin  and 
Smol’yakov  (1978),  assuming  C,  =  1000  wppm  (+);  Vdovin  and  Smol'yakov  (1981),  assuming  C,  =  1000  wppm-,  (-); 
Fontaine  etal.  (1992)  (X).  The  three  lines  correspond  to  exponents  of  0.2, 2.7,  and  0.857.  Details  of  the  previous 

experiments  are  provided  in  Table  4. 
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Figure  12:  cJC{  vs.  K  for  WSR-301  of  the  present  experiment  (U  =  6.65  (•),  13.2  (■),  and  19.9  (A)  ms  '  with 
previously  reported  results;  Fruman  and  Tulin  (1976)  with  C,  =  100  wppm  (O),  500  wppm  (A),  and  1000  wppm  (□); 
Vdovin  and  Smol’yakov  (1978),  assuming  C,  =  1000  wppm  (+);  Vdovin  and  Smol’yakov  (1981),  assuming  C,  =  1000 
wppm\  (-);  Fontaine  etal.  (1992)  (X).  The  three  lines  correspond  to  exponents  of  0.2,  2.7,  and  0.857.  Details  of  the 

previous  experiments  are  provided  in  Table  4. 


Figure  13:  Three  images  of  the  near-wall  polymer  concentration  fields  taken  in  the  fi lamented  region  where 
K>2A\  10'7  (WSR-301 ,  QJQ ,  =  10,  C,  =  4000  wppm ,  U  =  6.65  ms1).  Three  types  of  images  are  common:  1)  Images 
with  minimal  polymer  concentration  (left  image),  2)  Images  that  appear  to  have  a  thin  sheet  of  polymer  flowing  along 
the  model  surface  (middle  image),  and  3)  Images  with  large  amounts  of  Filaments  both  along  the  model  surface  and 
lofted  away  from  the  surface  (right  image).  The  flow  is  from  left  to  right,  and  the  dashed  line  denotes  the  location  of  the 
solid  wall.  The  stretched  Field  of  view  is  5  mm  (stream-wise)  by  2.25  mm  (wall  normal). 


Figure  14:  Three  images  of  the  near-wall  polymer  concentration  fields  taken  in  the  transition  region  between  highly 
filamented  and  uniform  regions  (5.3  x  10'8  <  K  <  2.4  x  10'7)  (WSR-301,  Q/Qs  =  2,  C,  =  4000  wppm ,  (J  =  6.65  ms1). 
This  flow  regime  is  characterized  by  periodic  filaments  in  a  somewhat  uniform  background  polymer  concentration.  The 
images  selected  show  various  examples  of  background  polymer  concentration  with  lofted  polymer  Filaments.  The  flow 
is  from  left  to  right,  and  the  dashed  line  denotes  the  location  of  the  solid  wall.  The  stretched  field  of  view  is  5  mm 

(stream-wise)  by  2.25  mm  (wall  normal). 
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Figure  15:  %DR/cM  vs.  cM  for  the  three  Mw examined:  U  =  6.65  (O),  13.2  (•),  and  19.9  (□)  ms  '.  Curves  fitted  to  Eq. 
7.4  are  plotted  for  each  speed  (dashed  lines)  and  for  all  speeds  together  (solid  line).  Values  of  intrinsic  drag  reduction 
and  intrinsic  concentration  derived  from  the  fitted  curves  are  provided  in  Table  6. 


Appendix  C: 


BUBBLE-INDUCED  SKIN-FRICTION  DRAG  REDUCTION  AND 
THE  ABRUPT 

TRANSITION  TO  AIR-LAYER  DRAG  REDUCTION 

BRIAN  R.  ELBING,  ERIC  S.  WINKEL,  KEARY  A.  LAY,  STEVEN  L.  CECCIO, 
DAVID  R.  DOWLING,  AND  MARC  PERLIN 

Department  of  Mechanical  Engineering  University  of  Michigan,  Ann  Arbor,  MI  48109,  USA 
(Submitted  to  the  Journal  of  Fluid  Mechanics  -November  2007) 

To  investigate  the  phenomena  of  skin-friction  drag  reduction  in  a  turbulent 
boundary  layer  (TBL)  at  large  scales  and  high  Reynolds  numbers,  a  set  of  experiments 
has  been  conducted  at  the  US  Navy’s  William  B.  Morgan  Large  Cavitation  Channel 
(LCC).  Drag  reduction  was  achieved  by  injecting  gas  (air)  from  a  line  source  through  the 
wall  of  a  nearly  zero-pressure-gradient  TBL  of  a  flat  plate  test  model  that  was  either 
hydraulically  smooth  or  fully  rough.  Two  distinct  drag  reduction  phenomena  were 
investigated;  bubble  drag  reduction  (BDR)  and  air  layer  drag  reduction  (ALDR). 

The  stream-wise  distribution  of  skin-friction  drag  reduction  was  monitored  with 
six  skin-friction  balances  at  downstream  distance-based  Reynolds  numbers  to  220  million 
and  at  test  speeds  to  20.0  ms'1.  Near-wall  bulk  void  fraction  was  measured  at  twelve 
stream-wise  locations  with  impedance  probes,  and  near-wall  (0  <  Y  <  5  mm)  bubble 
populations  were  estimated  with  a  bubble  imaging  system.  The  instrumental  suite  was 
used  to  investigate  the  scaling  of  BDR  and  the  requirements  to  achieve  ALDR. 

Results  from  the  BDR  experiments  indicate  that:  significant  drag  reduction  (> 
25%)  is  limited  to  the  first  couple  meters  downstream  of  injection;  marginal  improvement 
was  possible  with  porous-plate  versus  slot  injector  design;  BDR  has  negligible  sensitivity 


to  surface  tension;  bubble  size  is  independent  of  surface  tension  downstream  of  injection; 
BDR  is  insensitive  to  boundary  layer  thickness  at  the  injection  location;  and  no 
synergetic  effect  is  observed  with  compound  injection.  Using  these  data  previous  BDR 
scaling  methods  are  investigated,  but  data  collapse  is  observed  only  with  the  “initial 
zone”  scaling,  which  provides  little  information  on  downstream  persistence. 

ALDR  was  investigated  with  a  series  of  experiments  that  included  a  slow  increase 
in  the  volumetric  flux  of  air  injected  at  free-stream  speeds  to  15.3  ms1.  These  results 
indicated  that  there  are  three  distinct  regions  associated  with  drag  reduction  with  air 
injection:  Region  I)  BDR;  Region  II)  transition  between  BDR  and  ALDR;  and  Region 
III)  ALDR.  In  addition,  once  ALDR  was  established:  friction  drag  reduction  in  excess  of 
80%  was  observed  over  the  entire  smooth  model  for  speeds  to  15.3  ms'1;  the  critical 
volumetric  flux  of  air  required  to  achieve  ALDR  was  observed  to  be  approximately 
proportional  to  the  square  of  the  free-stream  speed;  slightly  higher  injection  rates  were 
required  for  ALDR  if  the  surface  tension  was  decreased;  stable  air  layers  were  formed  to 
free-stream  speeds  of  12.5  ms'1  with  the  surface  fully  roughened  (though  approximately 
50%  greater  volumetric  air  flux  was  required);  and  ALDR  was  sensitive  to  the  inflow 
conditions.  The  sensitivity  to  the  inflow  conditions  can  be  mitigated  by  employing  a 
small  faired  step  (10  mm  height  in  the  experiment)  that  helps  to  create  a  clean  separation 


line. 


1.  Introduction 


As  skin-friction  drag  is  an  important  component  of  resistance  in  nearly  all 
transportation  systems  moving  in  a  fluid,  skin-friction  reduction  techniques  have  been 
investigated  for  several  decades.  Prior  investigations  have  shown  that  skin-friction  drag 
reduction  can  be  achieved  by  injecting  gas  (air)  into  the  near-wall  region  of  a  liquid 
(water)  turbulent  boundary  layer  (TBL).  Bubble  drag  reduction  (BDR),  has  been  shown 
to  produce  local  drag  reduction  in  excess  of  80%  in  laboratory  scale  experiments,  and  has 
therefore  attracted  the  interest  of  numerous  investigators.  Unfortunately,  the  majority  of 
BDR  experiments  to  date  have  been  conducted  at  low  Reynolds  numbers  (on  the  order  of 
107,  based  on  downstream  distance)  and  small  scales  (typically  one  meter  or  less),  and  the 
proper  scaling  laws  remain  unclear.  For  practical  application  of  BDR  to  full-scale  ships, 
experimental  results  at  large  scales  and  high-Reynolds  numbers  approaching  those  of 
full-scale  craft  are  needed. 

High  Reynolds  number  experiments  on  BDR  were  reported  by  Sanders  et  al. 
(2006)  for  flows  of  a  near-zero  pressure  gradient  TBL.  They  showed  that  significant 
levels  of  BDR  could  be  achieved  near  the  injection  location.  However,  they  also  reported 
that  there  was  limited  persistence  of  the  BDR,  which  lasted  only  a  meter  or  two  from  the 
location  of  injection.  For  BDR  to  be  feasible,  the  propulsive  power  saved  by  drag 
reduction  must  exceed  the  power  requirements  for  gas  injection.  And,  the  further 
downstream  drag  reduction  persists,  the  less  overall  gas  must  be  injected.  However,  with 
such  low  levels  of  persistence,  BDR  would  require  a  large  number  of  gas  injection 
locations  and  a  large  volume  of  gas.  Researchers  have  speculated  that  the  gas  injection 
process  can  be  optimized  to  minimize  pumping  requirements  and  improve  BDR 


persistence,  as  the  geometry  of  injector  will  influence  the  size  and  speed  distributions  of 
the  injected  gas  bubbles. 

Sanders  et  al.  (2006)  also  observed  that,  at  lower  flow  speeds  and  higher  gas 
injection  rates,  a  layer  of  gas  would  form  on  the  underside  of  the  flat-plate  and  persist 
along  its  entire  length.  Such  air  layers  would  lead  to  skin  friction  reduction  of  greater 
than  80%.  Air  layers  were  observed  for  only  a  limited  number  of  flow  conditions,  as 
examination  of  this  phenomenon  was  not  the  principal  goal  of  that  study.  However,  the 
prospect  that  stable  air  layers  can  be  formed,  even  at  higher  speeds,  has  led  to  increased 
interest  in  this  phenomenon,  which  we  term  air  layer  drag  reduction  (ALDR)  (see 
Winkel,  2006). 

We  present  here  a  brief  review  of  BDR  and  ALDR  research.  A  seminal  study  of 
BDR  was  conducted  35  years  ago  when,  using  electrolysis  to  generate  bubbles  near  the 
leading  edge  of  an  axisymmetric  body,  McCormick  and  Battacharyya  (1973)  showed 
promising  results  for  skin-friction  drag  reduction.  Net  drag  reduction  approaching  40% 
was  observed  on  a  body  0.91  m  in  length  at  flow  speeds  to  2.6  ms'1.  In  another  early 
study,  Bogdevich  and  Evseev  (1976)  observed  that  drag  reduction  peaks  immediately 
downstream  of  the  point  of  injection;  however  further  downstream  skin-friction  returned 
to  its  usual  level.  Since,  numerous  BDR  experiments  have  been  conducted  that 
investigated  the  parameters  influencing  drag  reduction,  some  of  which  are  described 
below.  Gas  is  injected  typically  through  a  slot,  porous-plate,  or  an  array  of  holes  at  the 
wall  of  a  TBL  on  an  axisymmetric  body  or  flat  plate.  Through  this  research  a  number  of 
parameters  have  been  identified  which  influence  the  level  of  drag  reduction  including: 
downstream  distance  from  the  injector,  gas  flow  rates,  free-stream  velocity,  plate 


orientation  (buoyancy),  injector  geometry,  and  surface  roughness.  A  review  of  skin- 
friction  drag-reduction  can  be  found  in  Merkle  and  Deutsch  (1992). 

Researchers  at  the  Pennsylvania  State  University  have  conducted  a  number  of 
BDR  experiments.  Here,  we  mention  some  highlights  of  their  efforts,  and  subsequently 
we  refer  to  their  work  as  necessary.  Madavan  et  al.  (1984a,  b)  measured  drag  reduction  in 
a  laboratory  scale  flat-plate  TBL.  From  these  experiments  they  observed  that:  drag 
reduction  was  relatively  independent  of  injector  pore  size  (from  0.5  to  100  fim),  bubble 
size  distributions  are  determined  primarily  by  the  free-stream  velocity  and  gas  flow  rate 
(increasing  with  increasing  injection  rate  and  decreasing  free-stream  speed),  and  drag 
reduction  improves  when  the  bubbly  layer  is  beneath  the  solid  surface  of  a  plate  rather 
than  above  it.  Clark  and  Deutsch  (1991)  investigated  drag  reduction  on  an  axisymmetric 
body  with  imposed  pressure  gradients.  Their  results  showed  that  a  mild  favorable 
pressure  gradient  reduces  the  efficiency  of  BDR  and  that  in  adverse  pressure  gradients, 
gas  injection  leads  to  separation  and  high-levels  of  drag  reduction.  (Similar  results  were 
found  more  recently  by  Kawakita  and  Takano  (2000)  where  pressure  gradients  were 
applied  to  a  flat-plate  TBL  with  gas  injection.  Skin-friction  measurements  showed 
improvements  in  drag  reduction  with  the  application  of  adverse  pressure  gradients  near 
the  point  of  injection.)  Finally,  Fontaine  and  Deutsch  (1992)  used  five  different  gases 
with  varying  density  and  solubility  and  showed  that  the  type  of  gas  has  little  or  no  effect 
on  the  observed  drag  reduction. 

The  influence  of  bubble  size  on  BDR  processes  has  been  examined  by  several 
researchers.  McCormick  and  Battacharyya  (1973)  reported  drag  reduction  of  10%  to  30% 
with  small  bubbles  generated  by  electrolysis,  despite  estimated  mean  void  fractions  in  the 


boundary  layer  on  the  order  of  only  1%.  Kawamura  et  al.  (2003)  produced  bubbles  with 
two  different  methods:  the  first  involved  injecting  air  through  a  slot,  producing  bubbles 
from  0.5  to  2  mm  in  diameter.  The  second  enlisted  bubble  formation  through  out-gassing 
to  produce  bubbles  l/10'h  to  l/100th  the  size  of  the  injected  bubbles.  (Here  the  “plate-on- 
bottom”  configuration  was  employed  and  hence  their  finding  may  be  a  result  of  the 
increase  in  rise-velocity  of  the  larger  bubbles,  causing  them  to  migrate  from  the  plate 
surface  much  more  quickly.)  They  report  that  the  smaller  bubbles  are  nearly  twice  as 
effective  at  reducing  skin-friction  at  equivalent  void  fractions.  In  contrast,  Takahashi  et 
al.  (2001)  examined  BDR  on  a  towed  body  and  used  surfactants  to  reduce  the  bubbles 
size.  They  reported  no  apparent  difference  in  the  resulting  drag  reduction.  The  majority 
of  BDR  experiments  have  been  conducted  in  fresh  water.  If  bubble  size  is  an  important 
factor  in  drag  reduction,  equivalent  flows  of  salt  water  would  be  expected  to  have  a 
significant  influence  on  the  bubble  sizes  produced  and  the  resulting  levels  of  drag 
reduction.  Winkel  et  al.  (2004)  showed  that  bubble  formation  is  modified  and  the 
resulting  size  distributions  are  reduced  by  approximately  a  factor  of  four  (in  diameter)  in 
saltwater  relative  to  freshwater.  However,  direct  skin-friction  measurements  made  by 
Shen  et  al.  (2006)  with  gas  injection  into  both  fresh  and  saltwater  produced  no 
measurable  difference  in  drag  reduction.  In  most  studies  of  BDR,  the  bubble  sizes  are 
much  larger  than  the  viscous  length  of  the  TBL,  and  it  is  possible  that  the  injection  of 
very  small  bubbles  with  dimensions  on  the  order  of  the  viscous  length  may  lead  to 
enhanced  bubble  turbulent  interactions  and  drag  reduction.  However,  Shen  et  al.  (2006) 
injected  very  small  (on  the  order  of  10/„,  i.e.  10  viscous  wall  units)  lipid-stabilized 


bubbles  into  the  boundary  layer  and  saw  no  measurable  improvement  in  drag  reduction 
compared  to  BDR  flows  with  larger  bubbles  and  equivalent  void  fractions. 

A  group  of  Japanese  researchers  have  performed  experiments  at  large  scales  to 
directly  address  the  issue  of  BDR  scaling.  Long,  slender  flat-bottom  ship  models,  ranging 
from  12  to  50  m  in  length,  have  been  towed  at  speeds  to  7  ms'1  with  gas  injection 
(Watanabe  et  al.,  1998  and  Kodama  et  al.,  1999).  The  test  models  were  very  slender,  with 
a  beam/width  of  0.6  or  1.0  m.  Gas  was  injected  over  the  center  50%  of  the  model  span. 
Drag  reduction  was  observed  to  decay  with  downstream  distance,  and  by  the  end  of  the 
50  m  model  only  slight  drag  reduction  existed,  on  the  order  of  a  few  percent.  Overall 
skin-friction  reduction  approaching  20%  on  the  50  m  model  was  reported.  However,  it 
should  be  noted  that  there  were  no  “skegs”  or  “strakes”  to  contain  the  gas  underneath  the 
models,  thus  gas  may  have  escaped  from  the  sides  of  the  vessel.  Kodama  et  al.  (2002) 
reported  that  drag  reduction  measurements  made  in  the  span-wise  direction  showed  poor 
uniformity,  decreasing  from  a  maximum  at  the  model  centerline  to  nearly  zero  at  the 
outer  edges.  This  observation  agrees  with  the  concept  of  gas  escaping  from  underneath 
the  test  models. 

The  same  Japanese  research  group  has  attempted  at  least  two  BDR  experiments 
on  full  scale  vessels  exceeding  100  m  in  length.  Gas  was  injected  near  the  bow  of  a  ship 
through  horizontal  arrays  of  holes  at  flow  rates  to  43  m3s‘.  During  the  first  at-sea 
demonstration  (Kodama  et  al.,  2000)  no  net  power  savings  was  observed.  The  researchers 
reported  that,  although  skin-friction  reductions  were  thought  to  have  been  achieved,  an 
offsetting  reduction  in  propulsive  efficiency  occurred  because  of  gas  entrainment  in  the 
propeller.  In  more  recent  at-sea  demonstrations  (Nagamatsu  et  al.,  2002;  and  Kodama  et 


al.,  2006),  net  friction  reductions  and  power  savings  (after  gas  injection  costs)  of  a  few 
percent  were  demonstrated.  Due  to  the  scale  and  complexity  (surface  curvature,  surface 
roughness,  and  sloped  surfaces)  of  these  at-sea  demonstrations,  almost  no  details  of  the 
distribution  (stream-wise  and  span-wise)  of  the  gas  and  corresponding  drag  reduction 
were  measured. 

The  underlying  physical  mechanisms  of  BDR  are  not  yet  fully  understood.  It  is 
generally  agreed  that  for  bubbles  to  reduce  skin-friction  they  must  in  some  way  reduce 
the  turbulent  momentum  exchange  in  the  buffer  region  of  the  boundary  layer.  A  number 
of  different  mechanisms  have  been  proposed,  and  results  suggest  that  multiple 
mechanisms  work  simultaneously.  A  reduction  in  the  bulk  density  may  play  a  role  in  the 
shear-stress  reduction,  by  reducing  the  Reynolds  stresses,-  p(u'v')  in  the  near-wall 

region.  Moreover,  if  the  drag  reduction  effect  is  due  to  a  simple  density  reduction,  the 
observed  level  of  drag  reduction  should  scale  with  the  near-wall  void  fraction.  Lumley 
(1973,  1977)  hypothesized  that  the  presence  of  bubbles  locally  increased  the  viscosity  in 
the  buffer  region.  This  increase  in  local  viscosity  would  then  reduce  turbulent  fluctuation 
levels  and  increases  the  thickness  of  the  sub-layer  and  buffer  region  of  the  TBL. 
Similarly,  Pal  et  al.  (1989)  observed  that  not  only  did  mean  skin-friction  levels  decrease 
during  gas  injection,  but  fluctuations  in  skin-friction  were  suppressed.  In  contrast  to  the 
proposed  reduction  in  turbulence  levels,  using  PIV  Nagaya  et  al.  (2001)  observed  an 
increase  in  turbulent  fluctuations  with  gas  injection.  Although  increases  in  both  UO  and 
vQ  fluctuation  components  were  observed  during  drag  reduction,  they  de-correlated, 
effectively  reducing  the  Reynolds  stresses.  Finally,  Meng  and  Uhlman  (1998)  propose 
“bubble-splitting”  as  a  mechanism  for  BDR.  Their  argument  is  that  energy  is  extracted  as 


a  larger  bubble  is  split  into  two  or  more  bubbles;  the  surface  area,  As,  increases  resulting 
in  a  higher  surface  energy,  oAs,  where  o  is  the  surface  tension. 

Numerical  simulations  of  bubbly  turbulent  flows  have  been  conducted  that  also 
yield  insight  to  the  mechanism(s)  of  BDR  at  relatively  low  Reynolds  numbers.  Using 
direct  numerical  simulation  (DNS)  of  a  bubbly  shear-layer,  Druzhinin  and  Elghobashi 
(1998)  showed  that  bubbles  migrate  to  the  vortex  cores  in  a  turbulent  shear  flow  and 
reduce  turbulent  fluctuations.  Lu  et  al.  (2005)  examined  bubbly  channel  flow  and 
demonstrated  the  importance  of  bubble  deformation  whereby  entrainment  of  bubbles  into 
near-wall  vortices  suppressed  the  formation  of  stream-wise  vorticity  and  drag;  this  effect 
was  not  observed  for  non-deforming  bubbles.  Similarly,  van  den  Berg  et  al.  (2005) 
examined  friction  drag  reduction  in  a  Taylor-Couette  flow  and  also  found  a  significant 
influence  of  bubble  size  and  deformation  on  the  bubble-vortex  interactions  and  friction 
modification.  These  studies  provide  valuable  insight  into  the  scaling  and  mechanisms  of 
BDR. 

Though  BDR  has  been  studied  extensively  for  more  than  30  years,  a  universal 
scaling  law  that  adequately  collapses  the  experimental  results  has  yet  to  be  developed. 
The  drag  reduction  is  presented  as  a  percentage  drag-reduction,  %DR,  where  r w  is  the 
shear  stress  without  gas  injection  and  zw  is  the  (reduced)  skin  friction  with  gas  injection: 

%DR  =  1-—  (1.1) 

^ wo 

We  would  like  to  determine  how  %DR  scales  with  the  other  independent  parameters  of 
the  flow,  such  as  the  flow  speed,  (/„,  and  the  injected  gas  flux,  Qa>  and  other  flow 
parameters.  Madavan  et  al.  (1985)  report  %DR  as  a  function  of  the  mean  void  fraction  of 
air  in  the  TBL,  a  ,  defined  by  (1.1) 


(1.2) 
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where  Qw  is  the  volumetric  flux  of  the  liquid  through  the  unmodified  momentum 
boundary  layer  defined  by  (1.3)  where  <599  is  the  boundary  layer  thickness,  <5*  is  the 
displacement  thickness,  and  B  is  injector  span, 

Qw=U„B(8„-8').  (1.3) 


While  this  scaling  can  collapse  the  data  for  a  singe  set  of  specific  experiments,  it  does  not 
adequately  collapse  experimental  data  across  different  flows,  as  it  neglects  several 
parameters  which  influence  bubble  drag  reduction. 

More  recently,  Deutsch  et  al.  (2003)  and  Sanders  et  al.  (2006)  suggest  other 
scaling  relationships  for  BDR  on  rough  and  smooth  plates,  respectively.  They  developed 
the  scaling  parameter 


Q, 

Q„+u_Be0 


* 


(1.4) 


for  both  rough  and  smooth  plates,  where  6o  is  the  unmodified  local  momentum  thickness, 
and  ur  and  u*  are  the  friction  velocities  of  smooth  and  rough  surfaces,  respectively. 
Sanders  et  al.  (2006)  proposed  a  similar  scaling  law  for  BDR  on  smooth  surfaces, 


Qa 

Qa+U„B(d0mj-00)' 


(1.5) 


that  produces  a  better  collapse  of  the  data.  However,  it  should  be  noted  that  the 
improvement  in  the  collapse  is  due  likely  to  the  reduction  in  the  length  scale  [60  in]-  00]  < 

e0. 

As  opposed  to  BDR,  ALDR  occurs  when  a  continuous  or  nearly-continuous  layer 
of  air  is  formed  between  a  solid  surface  and  the  outer  liquid  flow,  as  observed  for  a 


limited  number  of  conditions  by  Sanders  et  al.  (2006).  Although  little  or  no  experimental 
investigations  have  been  conducted  to  investigate  specifically  the  phenomenon  of  ALDR, 
BDR  researchers  have  inadvertently  created  and  reported  some  of  the  characteristics  of 
ALDR.  Madavan  et  al.  (1985)  reported  a  maximum  level  of  drag  reduction,  typically  in 
excess  of  80-90%,  where  the  increase  of  gas  injection  yielded  little  or  no  improvement  in 
skin-friction  reduction,  likely  associated  with  the  formation  of  an  air-layer.  Kodama  et  al. 
(2002)  described  an  attempt  to  establish  an  air-layer  along  a  flat-plate  test  model  at  a 
free-stream  speed  of  10  ms'1.  Although  the  exact  gas  injection  rates  were  not  described, 
no  stable  persistent  air-layer  was  established.  They  described  an  air-film  that  persisted  20 
cm  downstream  of  injection,  and  subsequently  separated  into  slugs  of  air  and  eventually 
bubbles  further  downstream. 

In  the  present  work,  we  continue  the  study  of  BDR  and  ALDR  using  the  test 
model  employed  by  Sanders  et  al.  (2006).  We  report  on  our  continuing  effort  to 
understand  the  mechanisms  underlying  the  limited  persistence  of  BDR.  Furthermore,  we 
significantly  extend  our  examination  of  ALDR.  And,  finally,  we  examine  the  influence 
of  surfactants  and  surface  roughness  on  the  observed  levels  of  drag  reduction.  The  paper 
is  organized  in  the  following  sections:  we  discuss  the  experimental  methods  used  in  the 
collection  of  the  data  in  §  2  as  well  as  the  experimental  program;  the  results  of  our 
experiments  are  presented  in  §  3  along  with  a  discussion  regarding  the  scaling  of  BDR 
and  ALDR;  and  in  §  4  the  paper  is  summarized  and  conclusions  are  drawn. 

2.  Experimental  methods  and  the  experimental  program 

Three  separate  tests  were  conducted  in  the  LCC  on  the  same  test  model,  and  table 
1  provides  a  summary  of  the  experimental  setups.  (For  a  complete  description  of  the  LCC 


facility  and  its  use,  see  Etter  et  al.,  2005.)  The  schematic  of  the  test  model  working 
surface  for  the  primary  experiment,  Test  1,  is  shown  in  figure  1.  The  two  additional 
ALDR  experiments  investigated  surface  roughness  (Test  2)  and  inlet  condition  sensitivity 
(Test  3).  Test  2  used  the  same  layout  shown  in  figure  1  without  the  near-wall  imagers  and 
electrical  impedance  probes.  Test  3  included  a  1  cm  step  at  the  injector  that  was  faired 
upstream  to  the  model  nose;  the  shear  stress  was  integrated  over  sections  of  the  body 
using  six  load  sensors  spanning  nearly  the  entire  model  length.  As  in  all  the  experiments, 
a  coordinate  system  is  used  with  the  origin  located  in  the  plane  of  the  plate’s  leading  edge 
and  X  defined  downstream,  Y  defined  plate  normal  from  the  working  surface  of  the 
model,  and  Z  the  span-wise  coordinate.  The  suite  of  measurement  devices  included  skin- 
friction  force  balances,  in  situ  remote-controlled  laser  illumination  for  the  onboard 
photographic  imagers,  electrical  impedance  probes,  and  flow  metering.  Bubble 
photographs  and  electrical  impedance  measurements  were  made  in  the  near-wall  region 
of  the  test  model. 

2. 1.  Instrumentation 

2.1.1.  Gas  Injectors  and  Flow  Metering 

Gas  was  injected  through  the  lower  surface  of  the  test  model  at  two  locations,  X  = 
XinJ  =  1.38  and  3.73  m  where  Xinj  is  the  injector  downstream  location.  Two  different  types 
of  injectors  (slot  and  porous-plate),  spanning  the  center  2.65  m  of  the  test  model  (~  87% 
of  the  model  span)  were  tested;  cross-sectional  schematics  of  the  porous-plate  injector 
can  be  found  in  Sanders  et  al.  (2006)  while  the  slot  injector  used  in  Test  1,  Slot  A,  is 
shown  in  figure  2.  The  porous-plate  injector  consisted  of  a  slot  inclined  at  a  mean  angle 


of  25°  from  the  working  surface  of  the  test  model  and  contracting  at  a  full  angle  of  10°. 
The  slot  was  capped  by  a  layer  of  porous  (40  pun  mean  pore  diameter)  sintered  stainless 
steel  (Mott  Corporation),  2.0  mm  thick  extending  25  mm  in  the  stream-wise  direction. 
The  layer  of  sintered  metal  was  flush  with  the  working  surface  of  the  model  and  fastened 
to  the  injector  to  ensure  that  all  of  the  injected  air  passed  through  the  porous  material 
(Sanders  et  al.,  2006).  Slot  A  was  inclined  12°  from  the  working  surface  with  a  constant 
throat  gap  of  5.75  mm  that  produced  an  opening  in  the  stream-wise  direction  of  28  mm. 
Tests  2  and  3  used  the  Slot  B  design,  also  shown  in  figure  2,  which  had  the  slot  inclined 
at  a  shallow  5.7°  from  the  model  surface  and  contracting  at  a  full  angle  of  6.1°.  The 
downstream  edge  was  broken  giving  it  a  convex  surface  that  produced  a  1.5  cm  opening 
on  the  test  model  surface.  During  Test  2,  a  porous  material  was  used  in  the  contracting 
slot  to  produce  pressure  drop  close  the  surface  opening  (required  for  a  separate 
experiment  not  discussed  here).  Test  3  used  a  modified  injector  from  test  2  with  a  1.3  cm 
step  affixed  to  the  upstream  edge  of  the  injector  opening.  The  injector  opening  was  1.4  m 
downstream  of  the  leading  edge  of  the  plate.  The  step  was  constructed  of  PVC  and 
continued  forward  to  the  nose  of  the  model  where  it  was  faired  with  fiberglass  reinforced 
epoxy.  The  injector  from  Test  2  was  modified  slightly  with  the  installation  of  nine 
equally  spaced  vanes  in  the  injection  slot  to  maintain  its  shape.  Two  of  these  vanes,  along 
with  the  previously  installed  baffles  and  screens,  divided  the  injector  into  three  sections. 
Each  of  these  sections  had  its  own  manifold  as  well  as  airflow  regulation  (as  described 
below). 

The  gas  injection  rate  was  monitored  using  two  insertion  thermal  mass-flow 
meters  (640S  Steel-Mass,  Sierra  Instruments)  mounted  at  the  center  of  a  6.3  cm  inner 


diameter  straight  steel  pipe.  To  ensure  that  the  flow  was  fully  developed  at  the 
monitoring  locations,  the  flow  meters  were  located  30  inner  diameters  (1.90  m ) 
downstream  and  10  inner  diameters  (0.63  m)  upstream  of  any  junctions  in  the  line.  The 
flow  meters  were  factory  calibrated  for  the  range  of  0  -  820  SCFM  or  equivalently  0  - 
23.2  SCMM (Standard  m3  min')  over  an  operating  range  of  0-345  kPa  and  10-54  °C.  An 
analog  voltage  signal  from  each  flow  meter  was  recorded  simultaneously  with  the  skin- 
friction  data.  Gas  was  delivered  to  the  injectors  via  40  ports  evenly  spaced  along  a 
manifold  that  spanned  the  rear  of  each  injector.  Inside  the  manifold,  three  layers  of 
screens  and  baffles  were  employed  to  generate  a  pressure  drop  to  ensure  the  even 
distribution  and  uniformity  of  gas  flux  along  the  injector  span. 

2.1.2.  Skin-friction  Force  Balances 

Local  skin-friction  measurements  were  made  at  six  stream-wise  locations  using 
floating-plate-type  drag  balances;  their  locations  are  shown  in  figure  1 .  The  floating  plate 
is  15.2  cm  in  diameter,  0.79  cm  thick,  and  made  of  17-4PH  stainless-steel.  The  plate  was 
fixed  rigidly  to  a  beryllium  copper  flexure  that  was  instrumented  with  a  full  wheat-stone 
bridge  of  semiconductor  strain  gages.  The  drag  balance  and  its  housing  were  flush 
mounted  using  an  eight  point  leveling  system.  The  annular  gap  between  the  housing  and 
the  floating  plate  was  60  ±  20  fj.m.  The  semi-conductor  strain  gages  were  excited  using  a 
Vishay  signal-conditioning  amplifier  (Model  2310,  Vishay  Measurement  Group).  The 
sensor  outputs  were  amplified  and  low-pass  filtered  at  10  Hz  with  the  same  Vishay  unit. 
The  output  signal  was  recorded  at  50  Hz  with  a  National  Instruments  NI-DAQ  data 
acquisition  board  and  a  LabView  virtual  instrument. 


The  skin-friction  balances  were  calibrated  in  situ  with  loads  ranging  from  zero 
through  8.9  N.  The  loads  were  applied  to  the  floating  plate  via  a  suction  cup  affixed  to 
the  plate  and  a  tensioning  cable  attached  to  a  precision  load  cell  (Model  LCEB-5,  Omega 
Engineering).  Tension  was  applied  to  the  cable  by  moving  the  precision  load  cell  with  a 
linear  traverse.  Such  a  setup  was  necessary  to  eliminate  biases  that  were  found  if  the 
weights  were  simply  suspended  with  cables  and  pulleys.  The  tensioning  load  cell  was 
calibrated  in  the  vertical  position  using  precision  weights  prior  to  its  use.  Multiple 
calibrations  were  performed  on  the  skin-friction  sensors  to  confirm  their  repeatability  and 
to  assess  the  uncertainty  of  the  measurement. 

Test  3  used  a  different  setup  that  measured  the  integrated  skin-friction  over  six 
sections  (approximately  1.6  m2  each)  spanning  the  entire  model  length.  To  accomplish 
this,  precision  rails  and  linear  bearings  were  used  to  affix  the  model  skin  to  the  body  of 
the  model.  Each  of  the  six  sections  was  instrumented  with  force  transducers/load  cells  to 
measure  the  total  frictional  force  on  the  respective  section. 

2.1.3.  Bubble  Imaging  Systems 

Images  of  the  near-wall  gas-laden  boundary  layer  were  acquired  using  three 
monochrome  V. i-inch  format  CCD  cameras  (Model  CV-M10-SX,  JAI,  659  x  494  pixel 
resolution)  with  50  mm  focal  length  lenses  (HF-50HA-1B,  Fujinon)  and  40  mm  extension 
tubes.  To  minimize  blurring,  the  camera  shutter  speeds  were  set  to  1.25  fis,  creating  a 
maximum  image  shift  of  less  than  two  pixels  (<  20  fim  in  the  stream-wise  direction)  at 
the  highest  flow  speeds.  Images  were  acquired  at  25  Hz  with  three  analog  frame  grabbers 
(Model  IMAQ-1409,  National  Instruments)  through  a  LabView  interface.  As  much  as  4 
W of  blue-green  light  from  three  different  argon-ion  lasers  (Model  5300,  Lexel;  Innova  I- 


90,  Coherent;  and  Innova  70C,  Coherent)  were  used  to  illuminate  the  bubble  images.  The 
laser  light  was  delivered  to  the  bubble  images  using  multi-mode  fiber  optics  at  each  of 
the  three  stream-wise  locations. 

A  special  technique  was  necessary  to  view  the  near-wall  region.  Using  a 
periscope-like  assembly,  the  image  path  utilized  a  “bump-out”  prism  protruding  5  mm 
from  the  surface  of  the  test  model.  A  cross-sectional  schematic  of  the  near-wall  bubble 
imaging  system  is  shown  in  figure  3.  The  image  plane  was  normal  to  the  lower  model 
surface  and  aligned  with  the  flow.  The  imaging  regions  and  scales  were  slightly  different 
for  each  of  the  three  near-wall  cameras,  but  all  were  approximately  8  mm  (stream-wise) 
by  6  mm  (vertical).  The  imagers  were  focused  at  the  surface  of  the  prism  and  the  depth  of 
field  was  measured  to  be  2  mm  in  air  using  a  target  inclined  at  45°  and  noting  where  its 
image  was  in-focus.  It  should  be  noted  that  recorded  bubble  populations  may  be 
contaminated  somewhat  by  the  flow  disturbance  caused  by  the  proximity  to  the 
protruding  prism. 

2.1.4.  Electrical  Impedance  Probes 

Electrical  impedance  probes  were  mounted  on  the  model  surface  to  measure  the 
bulk  void  fraction  in  the  near-wall  region  of  the  flow.  (These  measurements  along  with 
those  from  the  near-wall  cameras  provided  corroborative  data  regarding  void  fraction. 
Due  to  the  difficulties  involved  with  measurements  of  this  type,  redundant  information 
was  deemed  necessary  not  just  advantageous.)  Each  probe  consisted  of  two  brass 
electrodes  soldered  to  signal  conductors.  The  brass  electrodes  that  were  to  be  in  contact 
with  the  flow  were  machined  flush  and  mounted  in  an  1 1.4  cm  diameter  non-conducting, 
flat  PVC  disk.  Each  disk  had  a  large  (3.2  mm  diameter  electrodes  with  6.4  mm  cross- 


stream  separation)  and  a  small  (1.6  mm  diameter  electrodes  with  3.2  mm  cross-stream 
separation)  pair  embedded.  The  purpose  of  the  two  sizes  was  to  detect  qualitatively  void 
fraction  gradients  normal  to  the  surface  as  illustrated  in  Cho  et  al.  (2005).  The  probes 
were  positioned  1.91  cm  upstream  and  1.91  cm  downstream  of  the  disk  center  for  the 
small  and  large  probes,  respectively.  The  disks  were  positioned  slightly  off  the  cross¬ 
stream  centerline  and  centered  at  12  stream-wise  locations  (X  =  1.07,  1.96,  2.59,  3.41, 
5.09,  5.94,  6.61,  7.43,  9.23,  10.05,  10.68,  11.50  m).  Each  probe’s  interior  surface  was 
encased  in  an  epoxy,  which  electrically  isolated  each  electrode  prior  to  contact  with  the 
flow. 

The  basic  circuit  for  the  electrical  impedance  probes  is  shown  in  figure  4  and 
consists  of  an  AC  excitation  voltage  ( Vs ),  reference  resistor  (Rref),  probe  electrical 
impedance  associated  with  the  flow  (Ze/),  stray  capacitance  from  the  input  and  output 
wires  (Zc/  and  Zc2,  respectively),  and  stray  capacitance  between  the  input  and  output  wires 
(Zc3).  The  stray  capacitance  between  the  input  and  output  leads  was  minimized  by 
shielding  each  separately,  but  is  included  due  to  its  significant  influence  during 
measurement  of  higher  void  fractions.  The  input  and  output  stray  capacitance  could  not 
be  neglected  due  to  the  length  of  the  wires  required  for  the  experimental  setup.  This  stray 
capacitance  was  measured  in  situ  by  short  circuiting  the  probe  with  a  known  resistor  and 
measuring  the  voltage  across  Rrep  Vref.  The  stray  capacitance  was  measured  repeatedly 
throughout  testing  and  showed  negligible  variation.  The  circuit  AC  excitation  voltage 
was  produced  with  a  signal  generator  (8904A  multifunction  synthesizer,  HP),  and  the 
reference  resistor  had  a  nominal  resistance  of  75  kQ,  which  maximized  the  output  signal 


and  sensitivity.  Applying  Kirchhoff’s  current  law  to  the  circuit  provides  the  relationship 
between  the  electrode  impedance  and  the  known  parameters: 

2  fc-VK/Wrf  (2.1) 
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Four  impedance  measurement  systems  were  fabricated  to  facilitate  simultaneous 
measurements.  During  each  experiment,  measurements  were  acquired  from  each  pair, 
necessitating  a  system  for  rapid  probe  multiplexing,  performed  using  a  128-channel 
multiplexer  (PXI-2530,  National  Instrument)  that  was  capable  of  switching  between  each 
set  of  four  probes  at  more  than  200  Hz.  The  sampling  time  for  each  probe  was 
approximately  three  seconds,  the  data  were  acquired  at  1000  Hz.  Data  acquisition  was 
initiated  once  the  skin-friction  sensors  had  achieved  steady  state. 

Increased  probe  sensitivity  was  accomplished  by  recording  the  deviation  of  the 
measured  impedance  from  the  baseline  impedance  at  zero  void  fraction.  A  lock-in 
amplifier  (LIA)  (SR830  DSP  LIA,  Stanford  Research  Systems)  and  a  signal  generator 
(8904A  multifunction  synthesizer,  HP)  were  used  to  excite  and  balance  the  bridge  as  well 
as  to  demodulate  the  resulting  signal.  The  typical  output  from  the  LIA  was  the  amplitude 
and  phase  of  Vref-  Vw  with  a  gain  of  20,  averaged  over  3  ms,  and  band-pass  filtered  (0.01 
Hz  about  the  center  frequency).  Each  circuit  had  a  separate  excitation  frequency  (5,  7,  9, 
and  1 1  kHz).  A  single  computer  controlled  all  the  probe  multiplexing,  Vw  amplitude  and 
phase  variation,  and  data  acquisition  via  a  LabView  virtual  instrument  and  data 
acquisition  card  (PCI-6040E,  National  Instruments). 


2. 2.  Data  A  nalysis  Techniques 


2.2.1.  Analysis  of  Skin-friction  Drag  Reduction 

The  drag  reduction  is  presented  as  a  percentage  drag-reduction,  %DR,  defined  by 
(1.1),  where  is  the  shear  stress  without  gas  injection  and  rw  is  the  (reduced)  skin 
friction  with  gas  injection.  As  gas  was  injected  into  the  confined  volume  of  the  LCC, 
even  with  improved  pressure  control,  the  static  pressure  in  the  test  section  increased. 
And,  these  pressure  excursions  led  to  baseline  drift  in  the  measured  shear  stress  as  well  as 
changes  in  the  injected  volume  flux  of  gas.  To  account  for  this  bias,  a  time-record  of  the 
static  pressure  on  the  upper  wall  of  the  test-section  was  recorded  simultaneously  with  the 
skin-friction  and  gas  injection  rate  data.  As  BDR  scales  with  the  volume  of  gas  injected 
into  the  boundary  layer  (Shen  et  al.,  2006),  the  mass-flow  rate  recorded  by  the  flow 
meters  (as  standard  volumetric  flow-rate)  was  converted  to  a  true  volumetric  flow-rate  at 
the  test  model  surface  using  the  ideal  gas  law.  In  addition,  the  increase  in  test-section 
pressure  resulted  in  a  decrease  in  the  mass  flow-rate.  Thus,  over  the  course  of  a  gas 
injection  test  (typically  30-60  s )  the  volumetric  injection  rates  decreased  by  as  much  as 
30%  (dependent  on  the  free-stream  velocity  and  gas  injection  rate),  which  can  have  a 
significant  effect  on  the  level  of  drag  reduction  achieved.  To  account  for  this,  the  time 
traces  of  shear  stress  and  volumetric  gas  injection  rate  were  segmented  into  0.5  s 
intervals.  The  drag  reduction  observed  during  each  interval  was  distributed  among  ten 
bins  based  on  the  volumetric  injection  rate. 


2.2.2.  Near-wall  Void  Fraction  Analysis 

The  region  of  influence  of  the  electrode  pair  can  be  defined  as  the  radial  distance, 
R,  from  the  model  surface  through  which  a  nominal  percentage  of  the  current  passes  (Cho 
et  al.,  2005).  Here  the  “influence  volume”  is  bounded  by  two  arbitrary  values  of  R/L 
(with  L  the  electrode  spacing),  through  which  50%  (R/L  =  0.85)  and  75%  (R/L  =  1.93)  of 
the  electrical  current  flows.  A  value  for  R/L  of  1.4  was  selected,  the  average  of  the  50% 
and  75%  values,  which  yield  influence  radii  of  4.4  and  8.8  mm  for  the  3.2  and  6.4  mm 
electrode  pair  spacing,  respectively. 

The  measured  impedance  of  (2.1)  can  be  related  to  the  volume  averaged  void 
fraction,  a,  using  Maxwell’s  mixture  model,  (2.2)  (Ceccio  and  George,  1996).  This 
model  provides  an  estimate  only  since  the  analysis  of  Hewitt  (1978)  following  Maxwell 
(1881)  assumed  a  uniformly  disperse  bubbly  mixture,  which  is  not  expected  in  the 
present  experiments. 
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Zm  and  Zw  are  the  impedance  of  the  mixture  and  water,  respectively.  Gg  (~  0)  and  ow  are 
the  conductivity  of  the  gas  and  water,  respectively.  For  the  current  study  the  complex  part 
of  the  mixture  impedance  is  assumed  to  be  negligible  since  crw  »  2  n  feLe0,  where  aw 
is  the  conductivity  of  water  (~3.5  £iS  cm'),  f  is  the  input  frequency  (~10  kHz),  eL  is  the 
non-dimensionalized  permittivity  of  the  water  (~80),  and  e0  is  the  permittivity  of  a 
vacuum  (8.85xl012  F  m‘)  (George  et  al.,  2000).  The  real  part  of  the  mixture  impedance 
(i.e.  the  conductivity)  dominates  when  aw  is  greater  than  0.44  fiS  cm'. 


2.2.3.  Bubble  Population  Measurements 

The  resolution  in  the  near-wall  video  images  (Y  ~  0  to  5  mm)  was  determined 
from  images  of  a  scale  recorded  periodically  throughout  the  test.  The  image  resolutions 
from  the  three  near-wall  cameras  (X-Xinj  =  0.58,  4.56,  and  9.30  m)  were  determined  to  be 
10.9  ±  0.2,  12.0  ±  0.2,  and  12.5  ±  0.2  / im  pixel'1,  respectively.  The  stated  uncertainty  in 
the  resolutions  is  the  95%  confidence  interval  determined  from  four  different  image 
calibrations  and  the  Student’s  r-distribution 

The  images  used  for  the  bubble  populations  were  selected  at  random  and  analyzed 
using  National  Instruments  Vision  Assistant  program.  The  “measurement  tool”  of  this 
data  analysis  software  was  used  to  draw  line  segments  across  the  (approximate)  diameter 
of  the  near  spherical  bubbles,  in  the  stream-wise  direction.  Only  in-focus  bubbles  were 
tabulated  in  a  method  described  in  Sander  et  al.  (2006). 

Approximately  1000  bubbles  were  counted  from  as  many  as  10  images  for  each 
flow  condition.  This  is  a  sufficient  number  of  bubbles  to  provide  an  accurate 
representation  of  the  test  data  as  the  histograms  are  consistent  on  an  image-to-image 
basis.  An  effort  was  made  to  count  all  the  in-focus  bubbles  in  each  image  to  eliminate 
any  random  error;  the  worst  possible  bias  error  is  estimated  at  two  pixels  (~25  pm)  for  the 
detection  of  each  edge  of  the  bubbles.  This  provides  a  bound  on  the  smallest  detectable 
bubbles  in  the  images,  about  25  pm,  and  while  bubbles  smaller  than  one  pixel  might 
appear  as  specks  in  the  images,  this  was  not  observed. 

The  physical  size  of  the  bubbles  observed  in  the  image  analysis  ranged  from 
approximately  25  to  1,000  pm.  The  resultant  bubble  diameters  were  placed  in  25  pm  bins. 


The  bubble  size  distributions,  N(d),  were  normalized  such  that  N(dl )  =  1 ,  where  Ld 
is  the  number  of  bins  and  i  is  the  bin  index. 

2.3.  Experimental  Program 

2.3.1.  Experimental  Test  Matrix 

In  Test  1  the  free-stream  liquid  velocity,  gas  injection  rate,  injection  location, 
background  surface  tension  and  injector  type  were  varied.  Free-stream  speeds  ranging 
from  6.7  to  20.0  ms'1  and  gas  mass  flow-rates  of  100  to  800  SCFM  (2.83  to  22.7  SCMM) 
with  a  maximum  volumetric  flux  of  0.14  m2  s1  per  unit  injector  span  were  tested.  These 
experiments  investigated  BDR  and  ALDR  from  a  single  injection  location,  possible 
synergistic  effects  produced  from  compound  injection  from  two  slot  injectors  2.35  m 
apart  in  the  stream-wise  direction,  and  the  influence  of  background  surface  tension  on 
BDR  and  ALDR.  To  test  this  latter  effect,  the  LCC  background  water  surface  tension  was 
reduced  from  about  70  to  50  dyne  cm1  by  adding  15  wppm  of  a  common  soluble 
surfactant,  Triton-X-100.  Research  conducted  by  Winkel  et  al.  (2004)  showed  that  the 
addition  of  surfactant  to  the  tunnel  volume  can  reduce  mean  bubble  diameters  by  more 
than  a  factor  of  two  immediately  downstream  of  injection.  However,  this  phenomenon 
was  not  observed  in  the  current  results,  which  indicates  that  it  is  limited  to  near  the 
injector. 

Test  2  investigated  ALDR  on  a  roughened  surface  in  two  parts:  Test  2a  repeated  a 
subset  of  conditions  from  Test  lb  to  determine  sensitivity  to  injector  design.  Then,  in 
Test  2b  the  entire  model  surface  was  roughened  using  epoxy  paint  (High  Build  Semi- 
Gloss  97-130,  Aquapon)  with  glass  bead  grit.  The  particles  were  tightly  packed 


producing  a  sand  grain  type  roughness.  Based  on  the  skin-friction  measurements  and 
assuming  a  fully  roughened  surface  the  first  75%  of  the  model  was  very  uniform  yielding 
an  average  roughness  height,  k,  between  400  and  600  fim.  The  remaining  25%  of  the 
model  was  rougher  with  k  ranging  from  800  to  1100  /rm.  Free-stream  speeds  ranging 
from  6.7  to  15.8  ms1  and  volumetric  fluxes  from  0.03  to  0.14  mV  per  unit  injector  span 
were  tested. 

Test  3  investigated  ALDR  via  an  alternate  drag  measurement  method,  and  the 
effect  of  injecting  air  behind  a  small  step  to  mitigate  the  effect  of  large  turbulent 
structures  at  the  injection  point.  Six  large  (1.6  m2)  floating  plates  were  used  to  measure 
integral  drag  on  the  plate,  as  opposed  to  the  six  small  floating  balances  used  previously. 
Free-stream  speeds  ranged  5.9  to  11.8  ms'1  and  gas  injection  rate  per  unit  injector  span 
varied  from  0.049  to  0. 16  mV. 

For  all  drag  reduction  experiments,  prior  to  each  injection  a  baseline  (non¬ 
injection)  skin-friction  measurement  was  made  for  a  period  of  about  15  seconds.  For  an 
injection  run,  shear  stress  data  acquisition  was  started  five  seconds  prior  to  gas  injection 
and  stopped  five  seconds  post  injection.  In  addition  to  the  shear  stress  measurements,  the 
gas  mass  flow-rate  to  one  or  both  injectors  and  the  static  pressure  of  the  LCC  test  section 
was  recorded  simultaneously.  A  control  valve  was  preset  to  provide  a  desired  gas 
injection  rate  according  to  the  pre-determined  test  matrix.  Two  manual  10.2  cm  diameter 
vents  were  added  to  the  top  of  the  LCC  test  section  to  help  regulate  pressure  rise  in  the 
test  section  during  a  gas  injection  run.  While  the  vents  improved  the  pressure  control 
there  were  still  variations  and  some  pressure  rise  did  occur  during  gas  injection. 


2.3.2.  Comparison  with  the  Setup  and  Methods  of  Sanders  et 
al.  (2006) 

The  results  presented  aim  to  extend  and  improve  the  data  set  collected  from  a  first 
phase  of  BDR  tests,  see  Sanders  (2004)  and  Sanders  et  al.  (2006).  All  phases  of  the 
present  experiment  are  meant  to  extend  the  current  BDR  research  to  larger  scales,  higher 
Reynolds  numbers,  and  higher  speeds  in  a  well-controlled  test  environment.  The  previous 
phase  of  experimentation  investigated  drag  reduction  using  injection  from  a  single 
porous-plate  injector  located  either  1.3  or  9.8  m  downstream  of  the  leading  edge.  The 
primary  experiment  investigated:  i)  the  influence  of  injector  geometry;  ii)  the  influence  of 
compound  injection  via  simultaneous  injection  from  two  slot  injectors  (X  =  1.38  and  3.73 
m)  for  a  variety  of  injection  conditions;  and  Hi)  the  influence  on  BDR  due  to  a  reduction 
in  surface  tension  from  the  addition  of  surfactant.  Both  experiments  investigated  similar 
ranges  of  flow  speeds,  from  6.7  to  20.0  ms'1,  and  gas  injection  rates  to  22.7  SCMM. 

In  addition  to  the  different  set  of  test  parameters,  the  present  experiment  provided 
higher  fidelity  measurements  than  those  of  Sanders  et  al.  (2006).  The  main  improvements 
are  in  the  control  of  tunnel  pressure  through  the  addition  of  test  section  vents  and 
simultaneous  measurement  of  that  pressure  and  the  gas  injection  rate.  As  gas  was  injected 
from  the  test  model,  the  vents  were  controlled  manually  to  allow  air  and  water  to  be 
purged  from  the  tunnel  to  accommodate  the  injected  volume  of  air.  Secondly,  the 
accuracy  of  air  injection  rate  and  tunnel  pressure  were  improved  by  recording  the  time 
traces  with  the  skin-friction  balances,  whereas  for  Sanders  et  al.  (2006),  static  pressure 
and  air  flow-rate  were  monitored  visually  from  a  digital  display.  The  improved 


measurements  of  the  present  experiments  allow  for  real  time  comparison  of  drag 
reduction  with  volumetric  gas  injection  rate. 

3.  Results 

The  presentation  of  results  is  divided  into  three  subsections:  Non-injection 
(baseline)  results,  BDR  results,  and  ALDR  results.  The  BDR  and  ALDR  are  separated  as 
they  are  fundamentally  differing  flows,  though  air  layers  can  arise  as  a  consequence  of 
bubble  injection.  These  subsections  are  further  divided  according  to  our  principal 
findings,  including  the  effects  of  various  parameters  and  the  examination  of  potential 
scaling  methods. 

3.1  Baseline  Results 

3.1.1.  Smooth  Model  Baseline  Flow 

Figure  5  shows  the  non-injection  skin-friction  coefficients  obtained  from  Tests  1 
and  2a  on  the  smooth  model,  where  Cfo  =  Two  /  \  pUj .  Also  plotted  with  the  current 

baseline  results  is  the  best-fit  curve  from  Sanders  et  al.  (2006),  and  the  Schultz-Grunow 
(1941)  curve: 


Cf0  =  0.370  log"2  584  ( Re* )  (3.1) 

where  Re*  =U„X/v,  X  is  the  distance  measured  form  the  leading  edge  of  the  test  model, 
and  v  is  the  kinematic  viscosity.  It  should  be  noted  that  Sanders  et  al.  (2006)  suspected 
their  absolute  skin-friction  measurements  were  biased  low  due  to  both  calibration  bias 
from  pulley  friction  and  possible  flow  through  the  sensor  housing.  The  current  data  were 


acquired  after  measures  were  taken  to  eliminate  these  effects,  as  described  above.  The 
current  results  indicate  that  Sanders  et  al.  results  were  indeed  biased  low  and  that  the 
present  results  from  the  same  model  are  in  good  agreement  with  the  Schultz-Grunow 
(1941)  friction  line. 

3.1.2.  Rough  Model  Baseline  Results 

Figure  6  shows  the  non-injection  results  obtained  from  Test  2b  on  the  roughened 
model  surface.  The  particles  embedded  in  the  epoxy  coating  were  tightly  packed  giving 
an  average  roughness  height,  k  ~  460  ±  120  /urn.  This  k  value  was  uniform  across  the  first 
nine  meters  of  the  plate,  however  along  the  last  four  meters,  the  coating  was  not  applied 
as  uniformly  resulting  in  approximately  double  the  k  value.  These  skin-friction 
coefficient  curves  were  used  also  to  estimate  the  average  roughness  height  by  assuming 
the  surface  is  fully  rough,  which  is  supported  by  the  Reynolds  number  independence  of 
the  skin  friction  coefficient.  The  surface  roughness  was  inferred  from  the  measured  drag 
after  examining  White’s  and  Schlichting’s  friction  curves  for  fully  rough  flow  over  a  flat 
plate  (White,  1991).  This  method  yielded  average  roughness  heights  of  400,  550,  580, 
1100,  and  830  /urn  at  X  =  1.96,  5.94,  7.43,  9.23,  and  10.7  m,  respectively.  These  results 
are  also  consistent  with  a  visual  examination  of  the  surface  embedded  with  the  particles. 

3.2.  Bubble  Drag  Reduction 

3.2.1.  Down-stream  Persistence  of  Bubble  Drag  Reduction 

Skin-friction  drag  reduction,  %DR,  is  shown  versus  downstream  distance  from  the 
injector  (X-Xinj)  in  figures  7,  8,  and  9  for  the  three  primary  test  speeds  (£/„  =  6.7,  13.3, 
and  20.0  ms'1).  Both  the  porous-plate  and  slot-style  injectors  are  showm  with  injection 


from  the  first  injector  (X  =  1.38  m ).  Provided  in  the  figures  for  each  condition  are  the 
volumetric  injection  rates  corrected  for  the  static  pressure  in  the  test  section.  Also 
presented  in  the  legend  is  whether  that  condition  corresponds  to  BDR,  ALDR,  or 
transition  between  BDR  and  ALDR.  These  regions  are  defined  subsequently,  but  are 
provided  here  for  clarity.  For  the  current  discussion  of  BDR  only  the  three  lowest 
injection  rates  (2.55,  4.81,  and  9.06  m3 min')  at  13.3  ms'1  and  all  20.0  ms1  cases  are  used. 
Thus  for  the  BDR  conditions  near  the  injector,  high  levels  of  drag  reduction  are  observed 
and  then  decay  rapidly  with  downstream  distance.  The  decreasing  drag  reduction  with 
downstream  distance  is  associated  with  the  migration  of  bubbles  from  the  near-wall 
region  where  they  are  effective  at  reducing  drag.  Shear  induced  lift  forces  from  the  mean 
velocity  gradients  in  the  boundary  layer  are  suspected  of  driving  the  bubbles  from  the 
solid  surface  resulting  in  poor  persistence  of  drag  reduction.  Simple  models  that  map  the 
trajectory  of  a  single  bubble  in  an  unmodified  TBL  have  been  presented  in  Meng  and 
Uhlman  (1989)  and  Sanders  (2004)  demonstrate  this  effect.  As  a  result,  at  higher  speeds 
nearly  the  entirety  of  drag  reduction  was  lost  after  the  first  two  meters. 

3.2.2.  Near-wall  Void  Fraction  Measurements 
Prior  research  has  shown  that  gas  distribution  (void  fraction)  in  the  boundary 
layer  during  BDR  has  a  peak  occurring  some  distance  (typically  a  few  hundred  wall 
units)  from  the  wall  (Merkle  and  Deutsch,  1990).  Furthermore,  Pal  et  al.  (1989)  have 
reported  that  bubbles  do  not  reside  in  the  viscous  sub-layer  region,  indicating  that  the 
void  fraction  asymptotes  to  zero  at  the  wall.  This  near  bubble-free  zone  in  the  near-wall 
region  has  been  termed  the  “liquid  layer”.  The  existence  of  a  liquid  layer  is  evidenced  by 
the  near-wall  bubble  imaging  cameras;  shown  in  figure  10  is  a  wall-normal  image  of  the 


bubbly  layer  at  X  =  5.94  m  and  LL  =  13.3  ms'  (where  the  measured  drag  reduction  was 
about  10%).  The  image  shows  a  dark  region  immediately  adjacent  to  the  plate  surface; 
the  dark  area  indicates  that  there  are  few  air-water  interfaces  ( i.e .  bubbles)  to  scatter  the 
laser  light. 

Electrical  impedance  measurements  also  indicate  that  void  fraction  increases  with 
increasing  distance  from  the  wall.  Only  the  electrical  impedance  measurements  on  the 
plate  surface  at  20.0  ms'  are  included  here.  Void  fraction  estimates  versus  X-Xinj  are 
shown  for  the  two  electrode  pair  spacings  (3.2  and  6.4  mm,  i.e.  small  and  large, 
respectively)  in  figure  11.  At  20.0  ms'1  the  void  fraction  measurements  beyond  the  first 
two  meters  were  on  the  order  of  a  few  percent.  Also  of  interest  is  that  the  void  fraction 
measurements  made  by  the  electrodes  with  the  larger  spacing  consistently  show  higher 
void  fraction.  The  larger  electrode  spacing  resulted  in  a  larger  “influence  volume”  and 
therefore  measured  the  impedance  of  a  larger  portion  of  the  boundary  layer.  From  this, 
we  qualitatively  infer  that  the  void  fraction  increased  with  distance  from  the  wall,  as 
discussed  by  Cho  et  al.  (2005).  This  result  agrees  with  the  results  reported  by  Nagaya  et 
al.  (2001).  In  figure  12  the  void  fraction  versus  X-Xinj  is  shown  for  the  four  nominal 
injection  rates  to  compare  the  two-injector  geometries  and  the  addition  of  surfactant,  all 
at 

20  ms'1  free-stream  speed.  Similar  results  were  obtained  with  the  porous-plate  injector 
with  and  without  surfactant  in  the  background,  and  the  near-wall  void  fractions  were 
consistently  higher  with  the  porous-plate  injector  than  with  the  slot  injector.  This  increase 
in  void  fraction  is  consistent  with  the  mild  improvements  observed  in  the  drag  reduction 
using  the  porous-plate  injector.  The  skin-friction  drag  reduction  is  shown  in  figure  13  as  a 


function  of  the  estimated  near-wall  mean  void  fraction  measured  by  the  impedance 
probes  at  20.0  ms'1  for  both  electrode  pair  spacings.  Both  are  reasonably  well  represented 
by  a  straight  line  with  the  smaller  electrode  spacing  having  approximately  double  the 
slope  of  the  large  spacing.  This  reinforces  that  bubble  stratification  away  from  the  wall 
occurs  near  the  wall. 


3.2.3.  Scaling  of  Bubble  Drag  Reduction 
The  scaling  laws  developed  from  previous  BDR  research  (described  in  §  1)  are 
evaluated  with  the  present  data  set.  First,  the  scaling  law  developed  by  Madavan  et  al. 
(1985),  which  assumes  that  drag  reduction  is  proportional  to  the  mean  void  fraction,  a  , 
across  the  entire  boundary  layer  is  applied  and  the  results  shown  in  figure  14  for  the  two 
injector  styles.  This  scaling  produces  a  rather  poor  collapse  of  the  experimental  results 
obtained  in  the  present  study.  Second,  the  more  recent  scaling  law  developed  by  Deutsch 
et  al.  (2003)  is  implemented  and  the  outcome  is  shown  in  figure  15,  where  the  length 
scale  employed  is  the  momentum  thickness,  da  rather  than  the  boundary  layer  thickness. 
Once  more  this  scaling  produces  a  rather  poor  collapse  of  the  current  data  set. 

To  improve  BDR  scaling,  we  posit  a  new  “initial  zone”  scaling.  Following 
Madavan  et  al.  (1984a),  we  begin  with  drag  reduction  as  a  function  of  the  dimensionless 
quantity, 


Qa 

U^S’ 


(3.1) 


where  S  is  the  area  through  which  the  gas  is  injected.  In  their  study,  gas  was  not  injected 
from  a  line-source,  but  through  the  surface  area  of  the  plate  spanning  the  tunnel  width 
and  23  cm  in  the  stream-wise  direction.  This  scaling  is  modified  to  represent  a  nominal 


air-layer  thickness,  ta,  given  by  (3.2)  by  assuming  that  the  injected  gas  is  moving  at  the 
free-stream  velocity,  where  B  is  now  the  injector  span  rather  than  the  area  of  gas 
injection: 


(3.2) 


Drag  reduction  measurements  taken  from  the  first  shear-stress  balance  ( X-Xinj  =  0.58  m) 
are  presented  in  figure  16  as  function  of  ta.  The  present  results  show  a  good  collapse  of 
the  data  and  agree  well  with  the  findings  of  Madavan  et  al.  (1984a),  whose  results  are 
also  shown  in  figure  16.  Although  a  reasonable  collapse  is  evident,  unfortunately  this 
scaling  yields  no  information  about  the  downstream  persistence  of  BDR. 


3.2.4.  Bubble  Drag  Reduction  with  Compound  Injection 
To  investigate  whether  any  advantage  is  obtained  with  the  combination  of  nearby 
injection  locations,  experiments  were  conducted  where  gas  was  injected  into  the 
boundary  layer  simultaneously  from  two  slot-type  injectors  spaced  2.35  m  in  the  stream- 
wise  direction.  Drag  reduction  was  measured  for  a  variety  of  nominal  gas  injection  rates 
at  the  three  nominal  test  speeds.  Aggregate  gas  injection  rates  from  both  injectors  were 
limited  to  22.7  SCMM ,  the  maximum  flow  rate  of  a  single  injector.  To  analyze  the  effect 
of  compound  injection  on  drag  reduction,  the  sum  of  drag  reduction  achieved  with  slot  1 
and  slot  2  independently  is  compared  with  the  drag  reduction  obtained  with  the  same 
flow  rates  using  compound  injection.  In  many  cases  it  was  found  that  there  was  little 
difference  between  the  two,  such  as  shown  in  the  upper  graph  of  figure  17  for  13.3  ms1. 
However,  the  lower  graph  of  figure  17  shows  a  unique  condition  at  13.3  ms1  where  for 
the  same  aggregate  injection  rate,  a  single  injector  yielded  much  larger  drag  reduction 


than  simultaneous  injection  from  two  slots.  Here,  injecting  the  full  gas  flux  from  slot  1, 
produced  near  100%  drag  reduction  over  the  entire  plate.  This  observation  is  a  result  of 
the  transition  from  BDR  to  ALDR,  which  will  be  discussed  later  in  the  paper.  In  contrast, 
no  stable  air  film  was  formed  when  the  gas  flow-rate  was  split  evenly  between  the  two 
injectors  -  the  drag  reduction  decayed  drastically  -  even  if  the  total  volume  was  larger 
than  the  single  injector  case.  These  results  indicate  that  there  are  no  synergistic  effects  of 
compound  injection  that  improve  efficiency  of  BDR  over  the  range  of  parameters 
investigated;  rather,  splitting  the  gas  injection  between  two  injectors  can  actually  reduce 
the  total  drag  reduction  by  preventing  a  transition  from  BDR  to  ALDR. 

3.2.5.  Influence  of  Injector  Geometry 

The  skin-friction  drag  reduction  results  presented  in  figures  7  through  9  indicate 
that  changes  in  the  injector  geometry  as  tested  has  only  a  weak  effect  on  either  the 
effectiveness  or  the  downstream  persistence  of  BDR.  However,  rather  consistently  at  all 
speeds  and  injection  rates,  slightly  higher  levels  of  drag  reduction  were  observed  with  the 
porous-plate  injector  than  with  the  slot  injector,  although  the  differences  lie  within  the 
uncertainty  of  the  measurements.  This  result  qualitatively  agrees  with  the  findings  of 
Madavan  et  al.  (1985),  who  showed  that  gas  injection  from  porous-plate  injectors  with 
two  different  pore  diameters  (0.5  and  100  /xm )  showed  measurable  yet  minor  differences 
in  drag  reduction. 

3.2.6.  Influence  of  Surface  Tension  on  BDR 

To  investigate  the  influence  of  bubble  size  on  drag  reduction,  a  surfactant  was 
employed  to  reduce  the  surface  tension  and,  in  theory,  the  bubble  size.  A  common 


soluble  surfactant  (Triton-X-100)  was  added  to  the  background  volume  of  the  test  facility 
at  a  concentration  of  15  wppm.  Winkel  et  al.  (2004)  showed  that  this  can  reduce  bubble 
size  by  a  factor  of  two,  at  least  near  the  injector.  The  surface  tension  of  samples  from  the 
tunnel  volume  was  measured  with  an  in  situ  tensiometry  device,  which  is  described  in 
Lapham  et  al.  (1999).  The  surface  tension  of  the  tunnel  volume  was  measured  to  be  70  ± 
1  dyne  cm1  prior  to  the  addition  of  the  surfactant  and  about  50  ±  2  dyne  cm1  during 
testing.  The  surfactant  employed,  Triton-X-100,  has  been  shown  to  have  time  dependent 
properties  by  Lapham  et  al.  (2001),  but  repeated  measurements  during  the  surfactant  test 
period  revealed  that  less  than  3  dyne  cm1  change  occurred.  Drag  reduction  at  13.3  and 
20.0  ms1  for  the  porous-plate  injector,  with  and  without  surfactant  in  the  tunnel  volume, 
are  compared  in  figure  18,  exhibiting  little  significant  difference  in  drag  reduction. 
Although  there  was  a  measurable  difference  in  surface  tension,  the  bubble  cameras 
showed  little  or  no  measurable  reduction  in  bubble  sizes  (table  2)  in  the  near-wall  region, 
in  contrast  to  the  findings  reported  by  Winkel  et  al.  (2004)  where  measurements  were 
made  immediately  downstream  of  injection.  However,  the  observation  that  drag  reduction 
performance  was  unaffected  by  the  addition  of  surfactant  agrees  with  the  findings  of 
Takahashi  et  al.  (2001).  In  addition,  histograms  of  the  bubble  populations  in  the  near-wall 
region  at  20.0  ms'1  are  shown  in  figure  19  for  the  5.66  and  22.7  SCMM  injection  rates  for 
the  porous-plate,  slot,  and  surfactant  (porous)  injection  tests.  This  shows  that  bubble  sizes 
downstream  from  the  injector  are  independent  of  injector  design  and  background  surface 
tension,  and  that  the  mean  bubble  diameters  increase  with  increasing  gas  injection  rate. 


3.2.7. 


Influence  of  Boundary  Layer  Thickness 


The  influence  of  boundary  layer  thickness  at  the  line  of  injection  on  drag 
reduction  is  shown  in  figure  20,  comparing  the  drag  reduction  results  for  a  single 
injection  rate  from  either  slot  1  (X  =  1.38  m )  or  slot  2  (X  =  3.73  m)  at  13.3  and  20.0  ms'. 
Although  the  data  points  are  rather  sparse,  the  data  appear  to  collapse,  indicating  that 
boundary  layer  thickness  at  the  line  of  injection  may  have  minimal  influence  on  the 
magnitude  or  persistence  of  BDR.  The  ratio  (slot  2  /  slot  1)  of  the  boundary  layer 
thickness  and  friction  velocities  are  about  two  and  one,  respectively.  Therefore,  while  the 
boundary  layer  thickness  did  change  for  the  two  cases,  the  resulting  drag  reduction  was 
not  strongly  influenced.  This  brings  into  doubt  scaling  based  solely  on  outer  variables, 
such  as  that  proposed  by  Madavan  et  al.  (1985)  and  Sanders  et  al.  (2006)  where  TBL 
thickness  is  a  length  used  to  scale  %DR. 

3.3.  AIR  LAYER  DRAG  REDUCTION 
3.3.1.  The  Abrupt  Transition  to  Air  Layer  Drag  Reduction 
Under  certain  conditions,  drag  reduction  approaching  100%  was  achieved  over 
the  entire  downstream  length  of  the  test  model  (approximately  10  m  in  length).  This  high 
level  of  drag  reduction  resulted  from  the  formation  of  an  air-layer  between  the  model  and 
the  liquid  flow.  The  formation  of  an  air-layer  is  evidenced  by  dark  images  from  the  near¬ 
wall  cameras  as  none  of  the  laser  light  was  scattered  from  air-water  interfaces,  by  the 
saturation  of  the  electrical  impedance  probes  at  values  approaching  the  impedance 
measured  in  air,  and  by  the  qualitative  images  from  an  oblique  view  of  the  lower  plate 
surface  from  outside  the  tunnel  during  the  second  phase  of  the  present  experiments. 

Sanders  et  al.  (2006)  achieved  air  layer  drag  reduction  for  only  a  limited  set  of 
flow  conditions.  Reported  here  are  the  results  from  a  series  of  experiments  meant  to 


investigate  the  conditions  necessary  to  achieve  air  layers  and  hence  ALDR.  Experiments 
were  conducted  only  for  the  porous-plate  injector  (with  and  without  surfactant  in  the 
tunnel  background)  from  the  upstream  location  (X  =  1.38  m)  at  four  free-stream  speeds, 
nominally  6.7,  8.9,  11.1,  and  13.3 

ms'1.  ALDR  was  not  observed  at  speeds  higher  than  13.3  ms1  in  this  first  set  of 
experiments  due  to  the  limitations  in  the  capacity  of  the  gas  delivery  and  metering 
systems.  (In  §  3.2.3  additional  experiments  with  surface  roughness  added  are  discussed 
and  speeds  to  15.3  ms'1  are  reported.)  For  each  free-stream  speed,  the  gas  injection  was 
varied  slowly  (approximately  0.142  to  0.283  SCMM  s  ')  from  zero  using  a  manual 
controlled  valve,  until  a  maximum  level  of  drag  reduction  was  observed  over  the  entire 
length  of  the  test  model.  Once  a  maximum  level  of  drag  reduction  was  observed,  the  flow 
rate  was  decreased  slowly  to  zero  injection.  No  appreciable  differences  were  observed  in 
drag  reduction  values  when  increasing  or  decreasing  gas  flow  rates  were  employed,  i.e. 
no  hysteresis  was  observed. 

From  the  varying  injection  rate  experiments,  %DR  vs.  volumetric  gas  injection 
rate  per  unit  span,  q,  curves  were  generated  at  a  given  flow-speed  and  downstream 
location.  %DR  versus  q  is  presented  in  figure  21  for  four  sensor  locations  and  a  flow 
speed  of  11.1  ms1.  As  observed  in  BDR,  over  the  lower-range  of  gas  injection  rates, 
%DR  decreases  with  downstream  distance.  However,  above  a  critical  gas  injection  rate, 
qcri,  (see  figure  23),  high  levels  of  drag  reduction  were  observed  with  no  apparent  decay 
with  downstream  distance  over  the  entire  downstream  test  surface. 

To  investigate  further  the  transition  from  BDR  to  ALDR,  we  examine  more 
carefully  another  location,  X-Xinj  =  6.05  m,  not  shown  in  figure  21.  Three  distinct  regimes 


are  apparent  and  identified  in  figure  22  with  %DR  versus  q  for  1 1.1  ms  ':  Region  I  -  BDR 
zone  where  %DR  increases  linearly  with  gas  injection  rate;  Region  II  -  A  transition  zone 
between  BDR  and  ALDR  where  the  %DR  versus  q  slope  increases  dramatically  relative 
to  region  I,  but  still  exhibits  linear  change;  and  Region  III  -  ALDR  zone  where  the  level 
of  drag  reduction  has  reached  a  maximum 

(~  100%)  and  increasing  the  gas  injection-rate  shows  little  or  no  improvement  in  drag 
reduction  due  to  the  near  complete  elimination  of  friction  drag.  Two  gas  injection  rate 
thresholds  are  defined;  a  transition  threshold,  q„ms,  and  a  critical  gas  injection  rate  for 
ALDR,  qcrit.  These  gas  injection  rates  are  located  at  the  break  point  (abrupt  slope  change) 
between  regions  I  and  II  (q,rans)  and  between  regions  II  and  III  ( qcril )  as  seen  in  figure  22. 
Shown  in  figure  23  are  the  results  for  qcrit  and  qnans  from  the  primary  experiment  without 
surfactant  in  the  background  that  define  the  transition  between  drag  reduction  regimes. 
3.3.2  Scaling  of  Air  Layer  Drag  Reduction 
The  gas  injection  rate  thresholds  are  observed  to  be  influenced  strongly  by  the 
free-stream  velocity,  as  observed  in  figure  23  for  the  porous-plate  injector  without 
surfactant.  The  critical  gas  injection  rates  per  unit  span  with  and  without  surfactant  in  the 
tunnel  background  are  shown  in  figure  24  as  a  function  of  free-stream  velocity,  U„.  The 
lines  shown  in  figure  24  are  least  mean-square-error  quadratic  fits  to  the  experimental 
data.  Interestingly,  the  gas  requirement  for  ALDR  is  about  5  to  10%  higher  with 
surfactant  in  the  tunnel  background.  This  implies  that  the  interfacial  tension  plays  a  role 
in  the  formation  and/or  the  stability  of  the  air  layer.  Though  we  do  not  fully  understand 
this,  we  conjecture  that  the  reduced  surface  tension  (with  surfactant  present)  allows  the 
liquid-gas  interface  to  deform  more  readily  (decreased  pressure  with  the  same  curvature); 


hence  the  liquid-gas  interface  allows  larger  fluctuations  to  occur,  and  a  larger  flux  of  gas 
is  required  to  maintain  the  air  layer.  In  addition,  the  small  difference  seen  in  figure  18 
that  shows  less  %DR  with  surfactant  than  without  surfactant  may  manifest  itself  here. 

Additionally,  if  a  true  air-layer  is  formed  between  the  model  surface  and  the 
outer-flow,  the  nominal  thickness  of  the  air  layer,  ta,  can  be  estimated  using  (3.2), 
ta  =Qj4/BU00  ,  the  same  parameter  used  successfully  for  the  “initial  zone”  scaling  of 

BDR.  Although  no  measurements  of  the  velocity  of  the  gas-phase  were  recorded,  it  is 
assumed  presently  that  they  are  on  the  order  of  and  proportional  to  the  liquid  free-stream 
velocity.  Again,  when  the  critical  gas  injection  rate  is  plotted  with  the  corresponding 
nominal  thickness,  ta,  a  linear  relationship  is  observed  as  shown  in  figure  25.  And  as  can 
be  seen  in  the  figure,  values  of  ta  from  approximately  4  to  8  mm  were  required  to 
maintain  an  air-layer  over  the  entire  test  model. 

To  provide  perspective  for  our  ALDR  results,  we  resort  to  other  large-scale 
experiments.  The  large-scale  model  tests  on  flat  bottom  models,  to  50  m  in  length,  by  the 
Japanese  group  at  the  National  Maritime  Research  Institute  (i.e.  Kodama  and  co-workers) 
have  exceeded  the  gas  injection  rates  used  here.  However,  a  persistent  air-layer  with 
dramatic  drag-reduction  was  not  observed  over  the  length  of  their  models.  At  a  speed  of  7 
ms1 ,  with  ta  about  4  mm,  drag  reduction  exceeding  40%  was  not  observed  beyond  the 
first  few  meters  of  the  50  m  test  model  of  Kodama  et  al.  (2002).  The  difference  is 
possibly  the  result  of  span-wise  gas  diffusion  on  their  test  models.  The  test  models  used 
were  at  most  one  meter  in  width,  with  gas  injected  over  only  the  center  50%.  Kodama  et 
al.  (2002)  report  that  the  span-wise  uniformity  of  drag-reduction  was  quite  poor;  it 
decreased  linearly  from  the  centerline  to  almost  zero  at  the  edge  of  the  model.  In 


addition,  it  is  possible  that  gas  could  have  escaped  from  beneath  the  test  model  as  the 
models  had  no  “skegs”  to  help  capture  the  air  beneath  the  ship.  Contrastingly,  in  the 
present  experiments,  gas  was  injected  over  nearly  90%  of  the  test  model  and  it  was 
“trapped”  underneath  the  model  by  the  tunnel  walls  (i.e.  the  model  spanned  the  entire 
width  of  the  test  section)  perhaps  explaining  the  difference  in  results. 

3.3.3  Roughness  Effects  on  ALDR 
The  first  supplementary  experiment  to  better  understand  ALDR  was  that  of  the 
roughened  model  surface.  As  this  configuration  used  a  different  injector  design  (Slot  B 
rather  than  the  porous-plate),  data  were  collected  once  again,  as  described  above,  initially 
on  the  smooth  model  for  free-stream  speeds  ranging  from  6.7  to  15.3  ms'.  At  the  highest 
free-stream  speed,  15.3  ms'1,  ALDR  was  observed  on  the  skin-friction  balances,  but  the 
required  air  flux  exceeded  the  flow  meter  calibration  range.  At  the  conclusion  of  these 
additional  smooth  surface  experiments,  the  model  was  roughened  and  data  were  collected 
for  free-stream  speeds  ranging  from  6.8  to  12.5  ms  '.  Figure  26  shows  qcril  for  the  two 
smooth  model  tests  and  the  rough  model  experiment  with  outline  data  points 
corresponding  to  results  that  exceeded  the  calibration  range  of  the  flow  meter.  As  is 
evident  in  the  figure,  while  ALDR  was  observed  for  a  given  speed,  it  required  about  a  30 
to  50%  larger  volumetric  gas  flux  to  transition  than  did  the  smooth  plate.  Note  further 
that  the  lower  two  curves  in  the  figure  exhibit  reasonable  repeatability  -  these  first  and 
second  experiments  with  the  hydraulically  smooth  plate  were  conducted  about  a  year 
apart. 

Video  imaging  of  the  flow  captured  the  transition  to  and  from  BDR  and  ALDR  as 
the  flow  evolved.  The  imager  was  mounted  below  the  test  surface  and  had  a  view  of  a 


fixed  location  about  five  meters  downstream  of  the  leading  edge.  Figure  27  displays 
images  from  the  video  at  three  times  during  the  injection  process  with  a  free-stream  speed 
of  6.8  ms'1.  Figure  27  (left)  was  recorded  prior  to  air  injection  and  shows  the  roughened 
model  surface;  figure  27  (center)  was  recorded  prior  to  achieving  the  critical  volumetric 
flux  to  transition  to  ALDR  (/.  e.  during  BDR);  and  figure  27  (right)  was  recorded  following 
the  transition  to  ALDR.  During  BDR  (middle),  the  layer  adjacent  to  the  plate  is  opaque 
and  the  hatch  and  sampling  port  (small  rectangle  in  upper  right  used  for  a  separate 
experiment)  that  are  visible  in  the  left  image  are  no  longer  evident.  Once  transitioned  to 
ALDR  (right),  the  layer  becomes  near-transparent  again  as  evidenced  by  the  reappearance 
of  the  sampling  port.  These  images  support  the  concept  that  transition  from  BDR  to  ALDR 
corresponds  to  a  change  in  flow  regimes. 

3.3.4  The  Importance  of  Inflow  Conditions:  Use  of  a  Faired 

Step 

In  a  third  experiment  to  investigate  air  layers  (as  well  as  air  cavities),  we 
discovered  that  uniformity  and  lack  of  large  (with  respect  to  the  LCC)  turbulent  eddies 
(on  the  order  of  1  cm)  in  the  inflow  TBL  to  the  injection  location  was  required  to  obtain 
ALDR  with  flush  injection.  This  has  important  ramifications  for  implementation  of 
ALDR.  In  the  LCC,  the  disturbances  were  caused  by  an  attempt  at  in  situ  repairs  to  the 
test  model  surface  upstream  of  the  injector  that  resulted  in  a  highly  disturbed  TBL 
upstream  of  the  injector.  After  observing  that  we  could  no  longer  achieve  ALDR  at 
equivalent  gas  fluxes,  a  different  injector  configuration  was  implemented.  This  setup 
included  a  1  -cm  step,  faired  upstream  into  the  nose  of  the  plate.  Once  installed,  the  step 
provided  a  clean  separation  line  for  the  flow  between  the  gas  and  the  liquid.  This  allowed 


for  the  formation  of  air  layers,  although  at  increased  air-flux.  This  result  indicates  that 
care  must  be  taken  to  insure  that  the  TBL  upstream  of  air-layer  formation  should  be  free 
of  large-scale  non-uniformity,  and  that  the  use  of  a  small  step  will  help  insure  a  clean 
flow  separation  for  the  air  layer,  even  if  the  incoming  TBL  is  not  highly  uniform. 

4.  Summary  and  Conclusions 

From  the  high-Reynolds-number  drag  reduction  experiments  conducted  on  the 
test  model,  the  following  conclusions  can  be  drawn  with  regard  to  bubble  drag  reduction 
(BDR).  1)  Prior  scaling  laws  developed  by  Madavan  et  al.  (1984a),  Deutsch  et  al.  (2003), 
and  Sanders  et  al.  (2006)  do  not  adequately  collapse  BDR  and  ALDR  data  for  the  range 
of  parameters  investigated  here.  2)  The  present  results  collapse  using  an  “initial  zone” 
scaling  and  agree  reasonably  well  with  the  results  of  Madavan  et  al.  for  the  sensor 
immediately  (0.6  m)  downstream  of  the  gas  injector.  3)  Drag  reduction  decays  rapidly 
with  downstream  distance  at  higher  free-stream  speeds  and  drag  reductions  exceeding  10- 
20%  were  not  observed  beyond  about  two  meters  at  20.0  ms1.  4)  Injector  geometry  has 
only  a  marginal  effect  on  the  magnitude  or  persistence  of  drag  reduction,  at  least  for  the 
two  different  geometries  here.  5)  Reducing  the  surface  tension  by  employing  a  surfactant 
has  little  effect  on  BDR  away  from  the  injector.  6)  Downstream  of  the  injector  decreasing 
the  surface  tension  has  negligible  effects  on  bubble  size.  7)  Changes  of  TBL  thickness  by 
a  factor  of  about  two  at  the  point  of  injection  has  little  or  no  influence  on  BDR.  And,  8) 
Compound  injection  (gas  injection  through  two  injectors  spaced  in  the  stream-wise 
direction)  showed  no  “synergistic”  effects  ( i.e .  no  improvement  over  injecting  the  entire 
flux  through  a  single  injector),  and  under  certain  conditions  compound  injection  can 
reduce  drag  reduction  by  preventing  the  transition  to  ALDR. 


The  following  conclusions  can  be  gleaned  from  the  air  layer  drag  reduction 
(ALDR)  experiments  conducted  on  the  test  model.  1)  Drag  reduction  can  be  divided  into 
three  distinct  regions:  Region  I)  a  BDR  zone,  where  DR  grows  linearly  with  gas  injection 
rate;  Region  II)  a  “transition  zone”  where  drag  reduction  increases  linearly  with  a  much 
steeper  slope  than  in  Region  I;  and  Region  III)  an  ALDR  zone  characterized  by  a  thin  air 
film  between  the  test  model  surface  and  the  liquid  free-stream  where  90%  ±  10%  drag 
reduction  is  observed  and  hence  increasing  injection  rate  yields  little  or  no  improvement 
in  drag  reduction.  2)  The  “critical”  gas  injection  rate  required  to  form  a  persistent  air 
layer  is  approximately  proportional  to  the  square  of  the  free-stream  liquid  velocity  (or 
linearly  proportional  to  the  nominal  air  layer  thickness,  ta  =  Qa/BU „  ).  3)  The  “critical” 
gas  injection  rate  was  observed  to  increase  slightly  with  surfactant  solution  in  the  tunnel 
volume.  A  drop  in  surface  tension  from  70  to  ~  50  dynes  cm  ',  increases  the  gas 
requirements  by  ~10%,  implying  that  surface  tension  plays  a  role  in  the  formation  and 
stability  of  an  air  layer.  4)  ALDR  was  observed  to  persist  over  the  entire  length  of  model 
(X-Xinj  ~  10  m)  at  speeds  to  15.3  ms1  with  no  decay  in  the  level  of  drag  reduction. 
Therefore,  it  is  conjectured  that  ALDR  may  have  persistence  lengths  much  greater  than 
the  length  of  the  current  test  model.  5)  Significantly  increasing  the  surface  roughness 
results  in  an  increase  in  the  critical  volumetric  flux  to  transition  to  ALDR,  but  a  stable  air 
layer  still  forms  and  was  observed  to  free-stream  speeds  of  12.5  ms1.  6)  Inflow 
conditions  significantly  affect  the  stability  and  formation  of  the  air  layer,  but  this 
deficiency  can  be  mitigated  through  the  use  of  clean  separation  lines  via  a  faired  step. 
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Tables  and  Figures 


Test  Type 

Model 

Injector 

Skin-Friction 

Measurement 

Test  la 

BDR 

Smooth 

Porous-plate  &  Slot  A 

Local 

Test  lb 

ALDR 

Smooth 

Porous-plate 

Local 

Test  2a 

ALDR 

Smooth 

Slot  B 

Local 

Test  2b 

ALDR 

Rough 

Slot  B 

Local 

Test  3 

ALDR 

(Inlet  Condition) 

Smooth 

Slot  B 
(Modified) 

Integrated 

Table  1:  Summary  of  tests  conducted  and  the  experimental  setup  employed  for  the  BDR 
or  ALDR  experiments.  The  model  surface  condition,  injector  type,  and  the  method  of 
skin-friction  measurement  are  listed. 


Upstream  location,  X  =1.96  m 


Injection  Rate 

2.55  m3minl 

4.81  namin'1 

9.06  m3min1 

15.3  m3minI 

Sintered  Metal 

123  ptm 

146  f, im 

130  pun 

138  pun 

Surfactant 

135  pun 

105  pim 

N/A 

N/A 

Slot  1 

106  pun 

133  y um 

124  pun 

178  pun 

Downstream  location,  X  =  10.68  m 


Injection  Rate 

2.55  namin'1 

4.81  namin' 

9.06  m3minl 

15.3  m3min  l 

Sintered  Metal 

80  pun 

92  pun 

90  pun 

181  pun 

Surfactant 

76  pun 

88  pun 

94  pun 

149  pun 

Slot  1 

N/A 

99  pun 

99  pun 

209  pun 

Table  2:  The  mean  bubble  diameters  measured  with  the  near-wall  bubble  cameras  (Y  < 
5  mm)  at  the  upstream  and  downstream  locations  with  a  free-stream  speed  of  20.0  ms1. 
Bubble  diameters  are  tabulated  for  the  porous-plate  (sintered  metal)  injection,  porous- 
plate  with  surfactant  injection,  and  slot  A  injection. 


SS  -  Skat-Friction  Sensor 

NWC  -  Near-Wall  Conte ra 

SEP  -  Surfbce  Electrical-Impedance  Probes 


Figure  1.  Schematic  of  the  test  model’s  working  surface  showing  the  injectors  and 
instrument  layout  for  Test  1.  Gas  injectors  are  located  1.38  and  3.73  m  from  the  model’s 
leading  edge.  Six  skin-friction  force  balances  and  twelve  surface  electrical  impedance 
probes  are  positioned  in  the  stream-wise  direction.  Near-wall  imagers  are  positioned  at  X 
=  1.96,  5.94,  and  10.68  m. 


Figure  2.  Cross-sectional  schematic  of  (upper)  Slot  A  used  during  Test  1  and  (lower) 
Slot  B  used  during  Test  2.  Test  3  used  the  lower  injector  without  the  porous  material  in 
the  slot  and  a  slightly  modified  air  delivery  system.  Gas  was  distributed  to  the  injector 
through  40  -  12.8  mm  diameter  ports  evenly  spaced  along  the  span  of  the  injector 
manifold.  Three  layers  of  baffles  and  screens  served  to  create  a  pressure  drop  and  evenly 
distribute  the  gas  along  the  injector  span.  See  Sanders  et  al.  (2006)  for  the  screen  and 
baffle  specifications. 


Tunnel  wall 


Figure  3.  An  elevation-view  schematic  of  the  near-wall  imaging  system  used  in  Test  1. 
A  periscope  type  assembly  of  two  prisms  was  used  to  deflect  the  camera  line  of  sight  to 
the  near-wall  region.  The  image  area  was  illuminated  from  the  side  using  laser  light 
delivered  with  a  fiber  optic  cable. 


Figure  4:  Circuit  diagram  of  the  electrical  impedance  probes  used  in  Test  1.  Shown  is 
the  voltage  source,  Vs,  the  impedance  of  the  bubbly  flow,  Ze„  reference  resistor,  Rreft  and 
the  stray  capacitance  from  the  lead  wires  Zc„  Zc2,  and  Zc3. 


Figure  5:  Baseline  results  obtained  from  Test  1  and  2a  on  the  smooth  model.  Also  shown 
is  the  best  fit  curve  obtained  by  Sanders  et  al.  (2006)  on  the  same  model  and  the  Schultz- 
Grunow  (1941)  flat  plate  skin-friction  curve. 
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Figure  6:  Baseline  results  from  Test  2b  on  the  roughened  model  surface.  Data  from  each 
skin-friction  balance  are  plotted  separately  (X  =  1.96  to  10.7  m)  to  show  the  Reynolds 
number  independence.  Also  shown  is  the  Schultz-Grunow  (1941)  curve  that 
approximates  the  smooth  model  baseline  results. 
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Figure  7:  %DR  versus  X-XinJ  at  a  free-stream  speed  of  6.7  ms'  during  Test  1.  A 
comparison  of  the  Slot  A  (solid  symbols)  and  porous-plate  (open  symbols)  injectors  at 
the  four  gas  injection  rates  is  presented.  Given  in  the  legend  is  the  volumetric  gas 
injection  rate  corrected  for  test  section  static  pressure,  and  in  parentheses  is  whether  this 
injection  rate  corresponds  to  BDR,  ALDR,  or  the  transition  between  BDR  and  ALDR. 
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Figure  8:  %DR  versus  X-Xinj  at  a  free-stream  speed  of  13.3  ms'1  during  Test  1.  A 
comparison  of  the  Slot  A  (solid  symbols)  and  porous-plate  (open  symbols)  injectors  at 
the  four  gas  injection  rates  is  presented.  Given  in  the  legend  is  the  volumetric  gas 
injection  rate  corrected  for  test  section  static  pressure,  and  in  parentheses  is  whether  this 
injection  rate  corresponds  to  BDR,  ALDR,  or  the  transition  between  BDR  and  ALDR. 
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Figure  9:  %DR  versus  X-Xinj  at  20.0  ms‘  free-stream  speed  during  Test  1.  A  comparison 
of  the  Slot  A  (solid  symbols)  and  porous-plate  (open  symbols)  injectors  at  the  four  gas 
injection  rates  is  presented.  Given  in  the  legend  is  the  volumetric  gas  injection  rate 
corrected  for  test  section  static  pressure,  and  in  parentheses  is  whether  this  injection  rate 
corresponds  to  BDR,  ALDR,  or  the  transition  between  BDR  and  ALDR. 


Figure  10:  Photographic  image  of  the  near-wall  bubbly  flow  at  13.3  ms'  ( X  =  5.94  m), 
5.10  mJmin  1  injection  (~10%  %DR).  The  image  is  approximately  5x8  mm,  and  the  flow 
is  from  right  to  left.  The  dark  region  adjacent  to  the  wall  (that  is  completely  black  and  has 
been  identified  by  the  superposed  white  dashed  line)  exhibits  the  liquid-layer  that  is 
almost  entirely  depleted  of  bubbles. 


Figure  11:  Near-wall  void  fraction,  a,  versus  downstream  distance  from  the  injector  at 
20.0  ms1  for  the  porous-plate  injector  (Test  1)  measured  by  the  surface  electrical 
impedance  probes  for  different  gas  flow  rates  (rimin').  Two  electrode  pair  spacings  are 
shown:  3.2  mm  (small)  and  6.4  mm  (large). 


Figure  12:  The  near-wall  void  fraction,  a,  measured  with  electrical  impedance  probes 
versus  downstream  distance  from  the  injector  at  the  four  nominal  injection  rates  is  shown 
for  the  porous-plate  (with  and  without  surfactant  in  the  background)  and  upstream 
injection  from  Slot  A  without  surfactant  background.  These  data  from  Test  1  all 
correspond  to  free-stream  speeds  of  20.0  ms'1  and  the  large,  6.4  mm,  electrode  spacing. 
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Figure  13:  %DR  versus  a,  the  near-wall  mean  void  fraction  in  the  boundary  layer 
measured  by  the  surface  impedance  probes  for  the  porous-plate  injector  at  20.0  ms 1  from 
Test  1. 
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Figure  14:  %DR  versus  a  ,  the  mean  void  fraction  across  the  entire  boundary  layer, 
shown  for  the  two  injector  styles  and  the  three  free-stream  speeds  from  Test  1. 
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Figure  15:  %DR  versus  the  scaling  of  Deutsch  et  al.  (2003),  (1.4),  for  the  smooth  surface 
where  uT  =  u*.  Results  are  from  Test  1  with  porous-plate  and  Slot  A  injectors. 
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Figure  17:  %DR  obtained  from  Test  1  versus  downstream  distance  from  injector  2  ( XinJ2 
=  3.73  m)  at  13.3  ms'  with  gas  injection  from  slot  1  only,  slot  2  only,  and  simultaneously 
from  both  slots  at  (upper)  4.81  namin  '  and  (lower)  9.06  m3miri‘.  For  comparison  with 
compound  injection,  solid  lines  that  represent  the  simple  sum  of  drag  reduction  from  slots 
1  and  2  only  (i.e.  from  the  filled  squares  and  circles)  are  shown;  dashed  lines  represent 
compound  injection.  Also  shown  is  (upper)  9.06  m3miri 1  and  (lower)  15.3  m3min  ‘ 
injected  only  from  slot  1. 
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Figure  16:  The  first  skin-friction  balance  (X-XmJ=  0.58  m)  %DR  versus  an  “initial  zone” 
scaling  nominal  air-layer  thickness,  ta.  Data  are  presented  for  the  porous-plate  and  slot 
injectors  (Test  1)  along  with  the  results  of  Madavan  et  al.  (1984a),  where  X-Xinj  =  0.254  m 
measured  from  the  injector  to  the  downstream  edge  of  the  drag  balance. 
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Figure  18:  %DR  vs.  X  -Xinj  using  the  porous-plate  injector  with  and  without  surfactant  in 
the  background.  Results  from  Test  1  indicate  negligible  differences  in  drag  reduction 
with  the  reduction  in  surface  tension  from  70  to  about  50  dynes  cm1. 


Figure  19:  Bubble  size  distributions  at  20.0  ms1  free-stream  speed  with  the  porous-plate 
injector  with  and  without  surfactant  and  the  upstream  slot  injector  at  ( X-Xinj  =  9.3  m)  for 
(upper)  4.81  m3miril  injection  and  (lower)  15.3  m3miri'  injection.  Data  were  collected 
from  the  near-wall  bubble  cameras  during  Test  la. 


Figure  20:  %DR  versus  X-  Xinj  (m)  from  Test  la  from  Slot  1  (Xinj  =  1.38  m)  or  Slot  2  (Xinj 
=  3.73  m)  at  free-stream  speeds  of  13.3  and  20.0  ms'.  Open  symbols  represent  slot  1 
injection  and  closed  symbols  represent  slot  2  injection.  Results  indicate  that  boundary 
layer  thickness  at  the  point  of  injection  has  little  effect  on  drag  reduction. 


FIGURE  21:  %DR  versus  q,  the  volumetric  gas  injection  rate  per  unit  span  (mV')  from 
Test  lb  for  11.1  ms'1  at  four  downstream  locations.  The  solid  horizontal  line  at  %DR  = 
80%  is  the  threshold  used  for  defining  ALDR. 


Figure  22:  %DR  versus  q  ( m2s ')  from  Test  lb  11.1  ms'1  at  X-Xinj  =  6.05  m,  which 
illustrates  the  three  regimes  of  BDR  and  ALDR:  I)  a  BDR  regime  where  drag  reduction  is 
nearly  linear  with  gas  injection;  II)  a  transitional  region  with  a  much  steeper  slope;  and 
III)  an  ALDR  regime  where  a  maximum  level  of  drag  reduction  is  achieved. 


Figure  23:  Transition  gas  injection  rate,  q ^  and  critical  gas  injection  rate,  qcri„  ( m2s ') 
versus  free-stream  speed  determined  from  Test  lb  data.  These  two  curves  define  the 
boundaries  for  the  three  drag  reduction  regions:  I)  BDR,  II)  transition  region  between 
BDR  and  ALDR,  and  III)  ALDR. 


Figure  24:  Critical  gas  injection  rate,  qcri,  ( m2s '),  for  ALDR  for  the  porous-plate  injector 
without  (0)  and  with  surfactant  ()  in  the  background  tunnel  volume.  The  lines  represent  a 
least-squares  quadratic  fit  to  the  experimental  data  collected  during  Test  lb. 
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Figure  25:  Critical  air  layer  thickness,  ta  =  QjBU„  ,  as  a  function  of  free-stream  speed, 
C/„ ,  for  ALDR  with  injection  from  the  porous-plate  injector  with  and  without  surfactant 
in  the  tunnel  background.  These  data  were  from  Test  lb. 


Figure  26:  qCru  versus  U„  for  ALDR.  Outlined  symbols  correspond  to  data  points  that 

exceeded  the  flow  meter  calibration  range  and  are  only  estimates.  Smooth  model  ALDR 
experiments,  Tests  lb  and  2a  were  conducted  about  a  year  apart  with  different  injector 
geometry.  The  roughened  model,  Test  2b,  was  tested  immediately  following  Test  2a  with 
the  same  injector. 


Figure  27:  Images  recorded  during  Test  2b  of  the  roughened  model  surface  at  6.8  ms'1: 
(left)  no  injection;  (center)  bubble  drag  reduction;  and  (right)  air  layer  drag  reduction. 


